Contents 


Certificate i 

Acknowledgements iii 

Synopsis v 

List of Figures xxi 

List of Tables xxx 

List of Symbols and Abbreviations xxxi 

1 Introduction 1 

1.1 Preamble 1 

1.2 Technological Significance of MeV Ion Implantation 2 

1.3 MeV Implantation : Energy Deposition and Range 4 

1.4 MeV Implant Induced Damage and Defects : A Brief Review 6 

1.4.1 Point defects 9 

1.4.1. 1 Divacancy 11 

1.4. 1.2 Other vacancy associated defects 14 

1.4.1. 3 Self-interstitials 15 

1.4.2 Defect Migration and Clustering 17 

1.4.3 Extended defects 19 

1.4.4 Defects and amorphization by irradiation 21 

1.5 Motivation and Focus of The Present Work 23 


XVI 



CAPACITANCE BASED STUDIES OF DEFECTS 
IN DEEP BURIED LAYERS PRODUCED BY MeV 
HEAVY IONS IN SILICON 


by 

PRAVAT KUMAR GIR! 



DEPARTMENT OF PHYSICS 

Indian Institute of Technology Kanpur 

OCTOBEe, 1997 



CAPACITANCE BASED STUDIES OF DEFECTS 
IN DEEP BURIED LAYERS PRODUCED BY MeV 
HEAVY IONS IN SILICON 


A Thesis Submitted 

in Partial Fulfilment of the Requirements 
for the Degree of 
DOCTOR OF PHILOSOPHY 


by 

PRAVAT KUMAR GIRI 


to the 

DEPARTMENT OF PHYSICS 

Indian Institute of Technology Kanpur 


OCTOBER, 1997 



XVll 


2 Principles of Experiments and Techniques of Data Analysis 2 1 

2.1 Introduction 

2.2 Defect Study by C-V Technique 2^ 

2.2.1 Defect profile in presence of uniformly distributed traps 29 

2.2.2 Defect profile in presence of nonuniform distribution of traps .... 32 

2.2.3 Factors affecting C-V profiling technique 33 

2. 2. 3.1 Series resistance 33 

2. 2. 3. 2 Excess leakage current 36 

2. 2. 3. 3 Minority carrier and interface trap 36 

2.3 Current- Voltage (TV) Characteristics 37 

2.4 Thermally Stimulated Capacitance (TSCAP) 38 

2.5 Capacitance Transient 39 

2.5.1 Exponential transient 40 

2.5.2 Nonexponential transient 40 

2. 5. 2.1 Series resistance effect 41 

2. 5. 2. 2 Leakage current 43 

2.6 Deep Level Transient Spectroscopy (DLTS) 43 

2.6.1 DLTS lineshape analysis 45 

2.6.2 Effect of series resistance and leakage current 46 

2.6.3 Trap parameter extraction 47 

2.7 Time Analyzed Transient Spectroscopy (TATS) 48 

2.7.1 Analysis of nonexponentiality : A case study 52 

2.7.1. 1 Typical current transients and PITS 53 

2.7.1. 2 TATS of positive peak 57 

2.7.1. 3 TATS of negative peak 59 



2 0 JUL 

l.(DKAKI 
} i. T., KANPUa 

128-5g4 

TH 


XXVll 


6.3 First order TATS spectra at 211.6K for different filling times (tf) showing 
progressive changes in occupancy for low dose (5 x 10^^ cm~'^) implanted 
sample. Peak Pi is fitted (dotted line) to two different centers PI A and 
PlB. Peak Dl can be fitted to a Gaussian broadened peak with a FWHM 

of 6 meV (not shown) 152 

6.4 Occupancy of individual peak as a function of filling time for samples with 

two different doses of Ar+ ions : (a) 5 x 10^^ cm'^ and (b) 1 x cm-^. 
Arrows indicate the measured emission time constants for the correspond- 
ing center 

6.5 Simulated occupancy behaviour of three trap levels by solving rate equa- 
tions under the condition n, (a) assuming three independent trap 

levels, (b) assuming that fastest and slowest sates are two configurations 

of the same defect. Parameters of simulation are given in box 166 

6.6 CC-TATS spectra for different filling time (t/) at a fixed temperature show- 

ing evolution of a new peak whose time constant increases with increasing 
filling times 171 

6.7 Change of emission time constant (r) with filling time (ty) for the evolved 

peak of Fig. 6.6. This shows stabilization of the relaxed state energy (E cc 
ln(r)) for longer filling times 172 

6.8 CC-TATS spectra at different temperature for a fixed filling time (t/=ls) 

showing temperature independence of the evolved peak 173 

6.9 Power law dependence of r and t/ for spectra shown in Fig. 6.6. Solid line 

refers to a best fit to relation : r ~ (l+t//te//)“ 174 

6.10 CC-TATS spectra with varying filling time at a fixed temperature for Au'^ 
irradiated (5 x 10® cm~'^) n-Si. Multiple peak structure in evolved peak for 


longer filling times can be noted 177 

7. 1 Temperature dependence of C-V characteristics of Ar'*’ irradiated p-Si Schot- 

tky diode 180 








CERTIFICATE ' 






It is certified that the work contained in the thesis entitled Capacitance based 
studies of defects in deep buried layers produced by MeV heavy ions in silicon 
by Mr. Pravat Kumar Giri, has been carried out under my supervision and has not 
been submitted elsewhere for a degree. 


October, 1997 




Dr. Y.N. Mohapatra 
Department of Physics 
Indian Instutute of Technolog}' 
Kanpur 



Introduction 


5 


in 

c 

111 


Matl. modification 
regime 
(Heavy Ions Ar,Kr, etc. 
of few KeV to few MeV) 
Ion implantation 
Ion beam mixing 
Ion beam assisted, 
deposition 


Matl. analysis regime 

of few MeV) 

Rutherford back scattering 
Channeling 

Elastic recoil detection 









X 


\ 


\ 


X 


2 3 ' 10 20 30 

is the energy in T.F. units) 


Figure 1.1: Universal stopping power vs. ion energy plot in solid in Thomson-Fenni (T.F.) 
unit as predicted by LSS theory. S„ and Sg are the nuclear and electronic stopping powers 
respectively (after ref. [14]). 

between swift ions and the target nuclei is small. Therefore, stopping power goes through 
a maximum at intermediate energies. Typical ranges of application for different energy 
ranges are also mentioned in the figure. 

There are other energy loss mechanisms, which are not so significant in the energy and 
mass regime of conventional ion implantation, such as elastic collisions with free electrons, 
inelastic collision with nuclei resulting from Cerenkov radiation etc [14]. 

Knowledge of energy loss mechanisms help determine the path that an implanted ion 
takes before coming to a halt. The ion trajectories being a statistical process, the resulting 
ion profile is a distribution, often assumed to be a Gaussian for low implant energies of 
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in electron and neutron irradiated silicon. Lee et al. [39] argue that the divacancy , 
A-center and divacancy-oxygen [V 2 + O] are equally probable in electron irradiated high 
oxygen content Si below 350°C. These anneal above 350°C giving the (divacancy + 2 
oxygen) [V 2 + O 2 ] and (trivacancy -f- oxygen) [I3 + O] centers. These in turn anneals 
out at 400°C yielding the (trivacancy + 2 oxygen) [V3 + 0\ and (trivacancy + 3 oxygen) 
[V3 + O] centers. Lee et al. have provided a comprehensive list of all such defects [39]. 

1.4. 1.3 Self-interstitials 

Self-interstitial is one of the most interesting defects in silicon because of its nearly ather- 
mal migration and ability to form pairs with impurities. Experiments to study directly the 
self-interstitial in silicon have not yielded much information on it’s configuration. Using 
EPR, Lee and Corbett found evidence for self-interstitial in neutron irradiated intrinsic Si 
at room temperature [40]. It appeared to be stable in a<(100)' di-interstitial configuration. 
Interstitial silicon atoms have been found to be highly mobile at 4.2K under electron irra- 
diation in p-Si [41]. In n-Si, its presence has been detected below 175K, though migration 
properties are less clear. In the absence of direct experimental observation of isolated in- 
terstitials, most of its properties are based on theoretical calculation. It is predicted to be 
a negative U center state binding two electrons deeply with first ionization state at Ec-0.8 
eV and second one at Ec-0.4 eV. The self-interstitials readily form pairs with impurities 
such as B, Al, Ga, C, P and 0 [42]. A number of EPR centers has been attributed to 
interstitial related defects (P6, A5, B3 and 02) [43] and their annealing characteristics 
are roughly established. Figure 1.4 shows schematic of thermal stability of several such 
interstitial related defect in Si as identified with EPR technique [44]. 

There has been resurgence of interest in silicon interstitial in order to account for defect 
processes quantitatively under irradiation conditions. We review some of these phenomena 
related to interstitial clusters later. In this section, in the context of the work described 
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of species. Gold ions are used as comparatively more massive particle than argon, and 
processing temperatures are kept well below that required for electrical activation. The 
issues that depletion layer capacitance based methods can illuminate best are also the 
ones least studied in such samples. We seek to study those aspects of defect phenomena 
which would complement simulational and structural studies through electrical signatures 
of defects created due to heavy ions. The questions and issues that these studies can 
potentially address are grouped below. 

• Are the electrically active defects created due to heavy ions in silicon, same as well- 
characterized point defects observed with light particle irradiation apart from being 
in larger concentration ? Or are distinctly new ones created ? Does agglomeration 
lead to change in the spectrum of density of states ? Do electrically active defect 
behave like point defects in such cases ? Which of these control compensation ? 

• Electrically active defects are sensitive to inhomogeneities in their environment. 
Does the damage in the environment of the defects lead to their energy broadening? 
Are they also sensitive to relaxation of the damaged layer brought about by say 
thermal treatment ? 

\ 

• The physics of the defects are best studied by monitoring trapping / detrapping 
kinetics, which also provide information on relaxation mechanisms. It is only to be 
expected that the disorder created by MeV heavy ions can sustain metastable states 
of defects. There are several such not so well understood cases of metastabilities 
reported in the literature in implanted silicon as we point out in more detail in 
chapter 6. Transient based spectroscopic methods are highly suitable to explore 
the nature of defect relaxation and metastabilities. We planned this work with 
the expectation that high energy, high fluence implants would be ideal laboratories 
for such studies. This was hoped to be specifically relevant to charge relaxation 
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curves are simulated for Q=0, 0.88, 1.16, 1.46 and 1.73 where Q^loTsCs. Note that for 
Qjl, the effect of series resistance is negligible. This effect becomes severe with Q>1. 
Consequently, wrong carrier concentration profile will be obtained from these curves. 
When series resistance is not negligible, then there is additional phase shift (j) between the 
ac voltage applied and ac current flowing through the device for capacitance measurement. 

Assuming that the diode being profiled can be represented by a simple RC equivalent 
circuit, if one measures the phase angle (p during profiling, true doping profile N is given 
by [88] 


with 


u = 


eA sin^ (f) 
Cm 


(2.5) 


If there is also shunt resistance in the equivalent circuit, the analysis becomes more com- 
plicated. Although shunt resistance does not affect capacitance measurement, it cause 
changes in the phase angle which could invalidate Eqns. 2.4 and 2.5. In practice, shunt 
resistance is easily detectable by monitoring the dc leakage current during profiling. In 
presence of high density of traps, the analysis is more complicated since both will have 
similar effect in C-V profiling. 

There are several possible ways to reduce the effect of series resistance since its presence 
is felt through the product {2Tr frsC)^. One could reduce the capacitance C by using test 
diode of smaller area. Because of importance of edge-effect, stray capacitance correction 
and difl&culty in measuring small area accurately, it is not practical to use diode areas 
below 10~^ cm^. The measurement frequency could be reduced to reduce the product, but 
the signal to noise ratio becomes worse at lower frequencies. In any case the frequency 
must be kept high compared to the relaxation times of the any traps present in the 
material. Series resistance can be altogether avoided by using a feedback circuit to keep 
the capacitance constant and measure the corresponding voltage. We have used this 



SYNOPSIS 


.Heavy ions accelerated to MeV energies are finding increasing use in semiconductor 
device technology.- Their use is likely to emerge as a standard industrial process for their 
ability to introduce, at distances larger than several microns from the surface, controlled 
modifications in doping, conductivity, lifetime, composition and degree of amorphization. 
- The success in the use of MeV heavy ions is closely linked to development of coherent 
understanding of defects and defect related phenomena accompanying these processes.-’ 

Though study of defects produced by energetic particle irradiation in semiconductors 
has been a field of intense research for many decades, there has been a recent upsurge 
in interest specific to heavy ions both from fundamental and technological point of view. 
The primary electrical defects produced by irradiation of light particles, such as electron, 
have been well characterized with various magnetic resonance spectroscopy and electrical 
techniques such as deep level transient spectroscopy (DLTS). On the other hand, defects 
and lattice damage induced by MeV heavy ions with high dose are comparatively much 
less understood. Our present knowledge of heavy ion induced damage is mostly limited 
to studies related to structural tools. For example, transmission electron microscopy has 
been extensively used to study evolution of structural defects during post-implantation 
annealing, and often complimentary information is obtained by Rutherford backscattering 
spectrometry studies. 

The use of conventional electrical characterization in heavy ion damaged semiconduc- 
tors, however, have been mostly limited to resistivity studies. Some recent studies have 
utilized DLTS technique to characterize defects in case of low dose implantation or in sam- 
ples subjected to post-implantation annealing at high temperatures. This thesis attempts 
to study electrically active defects in deep buried layers produced by MeV heavy ions in Si 
prior to removal of damage by high temperature annealing.- k variety of capacitance based 
techniques including transient spectroscopies such as DLTS, novel methods of isothermal 
spectroscopic techniques such as time analyzed transient spectroscopy (TATS), and their 
suitable extension for meaningful extraction of defect parameters, have been used to study 
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2.6.1 DLTS lineshape analysis 

DLTS lineshape is dependent on both temperature and energy of the trap for a chosen 
value of r. For an exponential transient, DLTS signal is given by 

S = Co[e-""‘‘ - 

where e„ = AT’^expl—Er/kT). S goes though a peak at temperature If Ti and T 2 
are temperatures corresponding to the half maxima positions of the DLTS peak then 

p-ent _ p-2ent _ ^ 

^ “8 

for a choice oi r = ^ = 2. 

Solving the quadratic equation for we get 

ln(e„i) - ln{en^) = 2.496 

2 ln&) - -^) = 2.496 (2.21) 

k ^Ti T 2 

Thus, for DLTS signal the linewidth can be defined by Eqn. 2.21 which depends on Et, 
Ti, T 2 and r factor. A prior knowledge of activation energy is required to verify the 
deviation from non-exponentiality from the observed spectra using Eqn. 2.21. This is a 
major limitation of DLTS in applying it for lineshape analysis. In deriving this equation, 
we have assumed that Co is not changing with temperature. For temperature dependent 
prefactor Co in Eqn. 2.9 , lineshape is not defined by Eqn. 2.21. Hence, lineshape analysis 
is not straightforward in case of DLTS. 

Usually, the first term in Eqn. (2.23) is negligible for exponential transients. In that 
case from the knowledge of Ti and T 2 , calculation of activation energy is straightforward 
from a single spectra. However, lineshape analysis is quite easy using isothermal spec- 
troscopic techniques such as time analyzed transient spectroscopy which is discussed in a 
later section. 



deep buried layers. 

In. order to separate out effects due to heavy ion damage alone from electrical activity 
of implanted species, we have used noble gas Ar'^ ions, with energies around 1.45 MeV 
and doses cm~^ to produce buried layers in both n and p-type epitaxial silicon. 

Samples implanted with heavier Au"^ ions at 4.5 MeV and dose 5x10^ cm but not 
annealed to avoid electrical activation, have also been studied for comparison. Irradiation 
doses are chosen to be below amorphization threshold but comparatively higher than used 
by most other workers. Conventional Schottky barrier diodes are fabricated for charac- 
terization using capacitance based techniques. Both as-implanted and low temperature 
annealed samples have been studied to monitor progressive changes in electrically active 
defects. Use of a variety of electrical spectroscopic techniques yield interesting physics 
of defect phenomena related to migration, clustering and carrier kinetics at deep levels. 
Specifically, defects primarily responsible for carrier compensation have been identified 
in both n-Si and p-Si. Origin of many interesting phenomena have been traced to the 
trapping kinetics of these defects within the depletion layer. 

The work has been presented in eight chapters. The first three Chapters deal with 
review and motivation, principle and critique of techniques employed, and experimental 
details respectively. Chapter 4, 5, and 6 present results on n-type samples, and Chapter 
7 is devoted to p-n^pe samples for comparison. Chapter 8 summarizes salient conclusions 
of this work. In what follows, we summarize contents of each chapter. 

Chapter 1 provides a brief introduction to the field of irradiation induced defects, 
defect processes and damage in Si with specific focus on MeV implantation. With this 
perspective, the motivation of the present study and its regimes of interest have been 
defined. 

Chapter 2 discusses principles behind use of various conventional experimental tech- 
niques and their extensions to analyze complex situations relevant to our present work. 
Limitations of standard use of capacitance-voltage (C-V) techniques and various factors 
affecting capacitance measurements are discussed. Common origin of nonexponentialities 
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Figure 2.7: T 3 ^icaJ photocurrent transients in SI-GaAs samples : (a) decaying current at 
174.9K, (b) rising current at 293.2K. 

believe that the current dipping to a minimum before rising is an essential part of the 
whole transient and this major deviation from exponentiality needs to be explained. Any 
analysis of this phenomena which focuses only on the later rising part, as has been done 
in the literature, can be misleading. 

The most popular way of characterizing traps in these materials has been PITS. The 
decaying transients give rise to positive peaks in PITS spectra. Transients such as given 
in curve (b) of Fig. 2.7 results in negative peak feature in such spectra and has been 
widely discussed in hterature. Figure 2.8 shows conventional PITS spectra in the higher 
temperature regime where the negative feature is clearly observed. This figure also shows 



in. capacitance transients and their effect on trap parameter evaluation from conventional 
deep level spectroscopies are assessed. 

The critical role of spectroscopic lineshape analysis of transient signals in the study of 
defect phenomena is emphasized and suitability of different spectroscopic techniques from 
this point of view is analyzed. Specifically, use of time analyzed transient spectroscopy 
(TATS), an isothermal spectroscopic technique, is shown to provide cleaner interpreta- 
tions and better insights into nonexponential transients which indeed manifest in \nrious 
forms in our experiments. A brief case study involving analysis of experimental data of 
photo-induced current transients in semi-insulating gallium arsenide has been presented 
to demonstrate superiority of TATS in analyzing nonexponential transients. Brief account 
of principles of other well known techniques which have been used for the present work is 
also included. 

Chapter 3 gives details of sample preparation, experimental setup and measurement 
procedures. Various steps in sample preparation including Schottky diode fabrication, 
MeV ion implantation and thermal annealing are described. Design of cryostat and ex- 
perimental setups developed and automated in-house have been described. 

In Chapter 4 we present results on steady state capacitance-voltage (C-V) measure- 
ments on Schottky diodes of n-Si containing damaged buried layers. Occurrence of a 
flat voltage-independent region in 1/C^ — V plot, for low reverse bias, is common to all 
irradiated devices. C-V measurements also show hysteresis which is dependent on voltage 
sweep rate. This in addition to temperature dependence of C-V curves is strongly indica- 
tive of trap controlled nature of observed phenomena. In as-implanted device, zero bias 
depletion width is seen to be very large compared to the control device. It is observed to 
be several microns larger than the range of ions predicted from Monte Carlo simulations 
using TRIM, and it is also larger for higher irradiation doses. The experimental features 
in C-V curves have been simulated using model charge profiles taking into account cross- 
ing of Fermi level with a midgap trap within the depletion layer. The simulations suggest 
that the irradiation damage seem to produce two regions : (i) a compensated region witi^ 
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ity of our model to mimic these essential features of the data opens up new possibilities 
as regards study of underlying mechanisms. We have assumed that e„ is the dominant 
thermally controlled parameter. However, it is possible that even the rate of conversion 
r is temperature dependent and hence the energy measured through Arrhenius plot is an 
effective energy associated with this combined process. The existence of metastable states 
with involvement of slow relaxation processes has often been invoked in relation to deep 
levels, as in the case of EL2 in GaAs, DX centers in AlGaAs and M center in InP. The 
occurrence of complex transient forms can be taken as possible indicator of involvement 
of metastable states even at such relatively high temperatures in these materials. 

In summary, we have demonstrated the efficacy of TATS in comparison to PITS, 
a temperature scanning spectroscopic technique, in analyzing nonexponential transients 
in SI-GaAs. Different kinds of nonexponentiality playing role at different temperature 
range are estimated using TATS. A phenomenological model to account for the observed 
nonexponentiality at higher temperature rising transients has been proposed. 

We will use TATS operated in constant capacitance mode (CC-TATS) to study a range 
of defect phenomena observed in heavily ion damaged layers in Si. CC-TATS will be an 
obvious choice to eliminate nonexponentiality due to series resistance effect and high trap 
concentration. Though we will use TATS and DLTS both to analyze data, CC-TATS will 
be used to study nonexponentiality due to other reasons than those discussed above. 



most dopants being deactivated, and (ii) a sharp negatively charged defect profile with 
large concentration close to the depletion edge which is much larger than expected range 
of ions. Our experimental observations lend support to recent results of molecular dy- 
namic simulations according to which primary defects such as vacancies and interstitials 
migrate and form clusters and complexes. 

A clear proof of occurrence of a defect-controlled near intrinsic region comes from 
measurements of forward current-voltage (I-V) characteristics of the irradiated de\'ice in 
which crossover from ohmic conduction to square law space charge limited conduction 
regime is observed. A nearly vertical rise in current at a voltage corresponding to trap 
filled limit is an indicator that the associated defect is a discrete level. 

C-V characteristics of samples after vacuum annealing of irradiated samples are also 
presented. They show progressive reduction in zero bias depletion width with isochronal 
(30 minutes) annealing upto 600°C. Reduction in depletion width on annealing is accom- 
panied with dopant activation in the damaged region. The nature of variations in charge 
profile with annealing has been shown using simulations of progressive changes in C-V 
characteristics. 

Chapter 5 deals with deep level spectroscopic investigation of defects in the buried 
layers. Some of the interesting and unusual features observed in transient response are first 
described. This includes occurrence of a ’spike’ in capacitance response during capture of 
carriers, small difference between reverse-bias and zero-bias depletion width following cap- 
ture, large magnitude emission transient and its premature termination. These features 
are shown to be consequences of trap filling giving rise to large negative charge density 
within the depletion layer and are discussed in the light of crossing between quasi-Fermi 
level and the major trap controlling these phenomena. DLTS, in spite of its inadequacy 
in this situation to provide reliable defect parameters through standard analysis, has been 
useful as a survey technique, and does provide helpful guide to further experimentation. 
Both DLTS and thermally stimulated capacitance (TSCAP) spectra in the range 90K- 
320K reveal occurrence of mainly two discrete trap levels of which one is the well known 




Figure 3.5: Layout of the home-made dipstic cryostat used for electrical measurement in the 
temperature range 80K-350K. 




Chapter 4 

Steady State Capacitance and 
Current Characteristics 


In this chapter, our central concern would be ways of obtaining useful information from the 
study of capacitance- voltage (C-V) characteristics of Schottky barriers containing damage 
due to implantation in the depletion width. In chapter 2, we have given a brief account of 
usefulness of conventional C-V measurements, also indicating major factors that influence 
such characteristics. Here we first present unusual and interesting features of results 
obtained from C-V measurements, and then describe their dependence on temperature 
and sweep rate of measurement, and different stages of annealing. We also present current 
- voltage (I-V) characteristics data to the extent needed as support for interpretation of 
C-V data. In the last section we present a discussion of the view that emerges from these 
results regarding the nature and role of defects in the damage layer, band bending and 
charge profile. The discussion is supplemented with results of model simulations using 
simple physically consistent assumptions to reproduce all the major observed features in 
the data. 


4.1 C-V Characteristics 

4.1.1 Typical C-V Characteristics 

A number of sets of Schottky diodes were made on unimplanted and implanted wafers 
using procedures described in Chapter 3. All the results presented in this chapter deal 
with n-type Si implanted with 1.45 MeV Ar^ ions and 4.6 MeV Au^ ions. The samples 
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divacancy center (V 2 ^°, Ec-0.42 eV) and the other is a newly found midgap trap level. 
This new trap level is the dominant trap occurring in very large concentrations in all our 
MeV irradiated samples and is responsible for controlling hysteresis in C-V and space 
charge limited conduction in forward I-V characteristics. In order to study the nature 
of nonexponentiality and proper evaluation of defect parameters, we resort to isothermal 
transient spectroscopic analysis in time domain using TATS. To avoid known sources of 
nonexponentiality due to large trap concentrations and series resistance, transient mea- 
surements for time domain analysis are performed in constant capacitance mode. 

TATS lineshape corresponding to the major trap is seen to be broad for low tem- 
peratures and, narrow and skewed at high temperatures. The peak shape narrowness is 
linked to occurrence of premature termination of transients which in turn is explained 
on the basis of time evolution of quasi-Fermi level and its crossing with the trap level 
within the depletion layer. It is possible to estimate the true time constant and peak 
height associated with these skewed peaks from a careful fitting of the useful part of the 
lineshape. The signature of the dominant trap is obtained from such analysis for various 
samples. The emission energy’ of the dominant trap is observed to vary from Ec-0.49 eV 
to Ec-0.56 eV depending on processing conditions such as low temperature oven annealing 
and irradiation dose. This sensitivity of the defect parameters seem to be characteristic 
of pramorphized regions and is due to varying degree of relaxation of disorder in its envi- 
ronment. At low temperatures, in the absence of skewness in lineshape due to premature 
termination, the damage related peak for high dose samples shows a Gaussian broaden- 
ing in energy with FWHM of 25 meV. This broadening is disorder induced since it is 
observed to be only 6 meV for low dose samples. The capture cross-section values are 
typical of neutral trapping center. The possible origin of this defect has been discussed, 
most likely candidates being dangling bonds, higher order vacancy complexes or intersti- 
tial clusters. However, consistent with our observation of involvement of migrating species 
and in accordance with mounting evidence from recent studies on involvement of inter- 
stitials, di-interstitials or their clusters seem to be responsible for the dominant damage 



related trap. High temperature furnace annealing studies show near complete removal 
of the dominant trap on heat treatment at 600 C for 30 minutes. In annealed samples 
minority carrier related traps could also be observed possibly due to mild inversion in the 
damaged layer. 

Chapter 6 is devoted to the study of trapping kinetics by monitoring relative oc- 
cupancies of multiple traps after allowing filling of traps for different duration over 5-6 
orders of magnitude in time. Surprisingly we observe two distinct mechanisms which 
could be isolated at different temperature regimes in different group of samples. .-Vt lower 
temperatures, filling time dependence of TATS peaks corresponding to multiple traps 
show nonmonotonic change in trap occupancy. With filling time evolution it wfis possible 
to recognize that the peak attributed to divacancy alone, in reality, constitutes contribu- 
tions from two distinct emitting centers. The resolution of peak structure hiis been clearly 
demonstrated using higher order implementations of TATS. The progressive evolution of 
occupancies over many orders of magnitude in time reveals coupling in capture kinetics. 
We show that this coupled kinetics is due to charge redistribution among multiple traps, 
in which deeper states gain at the cost of shallower states with time. Broad experimental 
features in occupancy as a function of filling time could be reproduced by solving cou- 
pled rate equations with charge redistribution constraint. However, closer match with 
experiment is obtained when in addition a configurational interdependence of two traps 
is assumed. The mechanism of charge redistribution had been invoked earlier to explain 
charge relaxation in amorphous silicon, though it was considered controversial in that 
context due to lack of unambiguous experimental proof. This mechanism is expected to 
be operative whenever number of trapping sites are larger than availability of carriers for 
capture. 

Evidence of a novel defect relaxation mechanism was encountered during the course 
of studies on filling time dependence of occupancies for low dose implanted samples {Au+ 
at dose 5 x 10 cm and for partially relaxed Ar"^ irradiated samples - presence of 
a preamorphised region with comparatively lesser degree of disorder being common to 
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rate not sufficiently fast and partly due to saturation in trap filling. The occurrence of 
negative slopes in tliese curves is indicative of profile in the trap concentration. Consider, 
for example, the filling curve denoted by solid lines during voltage change from reverse 
toward forward bias. Around zero bias, an increase in forward bias leads to so much of 
increase in negative trapped charge that the depletion width is forced to widen at the 
same bias to uncover additional positive charge at the edge. Hence the slopes are due to 
combined effect of trap occupation in the depletion layer and dopant charges at the edge 
of the depletion layer. 

The observation of hysteresis with different degree of filling and various sweep rates 
establish that C-V is controlled by large trap concentration due to irradiation damage in 
these samples. We go on to show their manifestation in I-V measurements in the following 
section. 


4.2 I-V Characteristics 

I-V measurements are useful for understanding the transport properties across junction 
of different materials. In our irradiated devices, a major deviation from conventional 
Schottky barrier (our control diode) behaviour occurs due to the presence of damaged 
layer. Such behavior is seen in both Ar"*” and Au"^ irradiated device. 

4.2.1 Typical I-V Characteristics 

Figure 4.8 shows typical forward and reverse TV characteristics at room temperature for 
Ar"*" implanted Schottky diodes. For irradiated device, both the forward and reverse 
current are significantly smaller compared to corresponding currents in the unirradiated 
device. This also confirms the presence of a compensating region in the depletion region of 
the device. For bias voltages corresponding to low forward bias, the magnitude of current 
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dynarnics of defect occupation. In case of a depletion layer without deep traps, voltage 
changes acros.s the device are accommodated through changes at the depletion layer edge. 
However, in case of large trap concentration, changes of occupancy at the Er-Eir crossing 
would control the total depletion width. 

From the experimental results it is abundantly clear that the occupied dominant traps 
produce charges opposite in sign to the background doping i.e. they get negatively charged 
after capture of electron. Hence, when reverse bias is increased (from, say, zero bias) the 
extra voltage drop can be accommodated either by increasing the depletion width, which 
leads to uncovering of positive charges at the edge, or by emitting electrons from the traps 
thereby increasing net effective positive charge in the depletion layer. The former process 
would occur within the time constant of dielectric relaxation time while the later would 
be controlled by trap emission rate. At sufficiently high temperature, if the trap emission 
rate is faster than voltage sweep rate, detrapping at the crossing would itself accommodate 
the voltage change. The width of the depletion layer and hence capacitance would remain 
constant even when bias is being changed causing the feature of flatness. For larger biases 
the trap concentration may not be sufficient to do so and width would also increase in 
part. This would continue till Fernai level loses contact with trap level. Beyond this the 
C-V would reflect the background concentration giving rise to conventional linearity in 
1/C^ vs. V curve. In this sense, Et-Ef crossing is the primary determinant in studies 
involving high dose implanted samples, while it is only a minor irritant in conventional 
depletion capacitance studies. 

This explanation of major features in C-V characteristics is fully consistent with (or 
can be viewed as derived from) the occurrence of hysteresis and temperature dependence. 
Hysteresis arises from the fact that the traps get occupied by a fast capture process at the 
crossing point during decreasing bias cycle while on the increasing cycle the capacitance 
change is dependent on slower emission process from the traps. Similarly, when temper- 


these set of samples. The phenomenon, as described below, could most conveniently 
be studied for temperatures above 270K. When filling time is increased at a particular 
temperature, TATS emission spectra shows a distinct new peak whose time constant is 
strongly dependent on the duration of filling. The time constant varies in a typical case 
from 7ms to 14s (corresponding to energy deepening of 180 meV) for filling times in the 
range 20ms-60s, showing a power law dependence on filling time. This dependence of 
emission characteristics on history, in terms of waiting time, is a strong indicator of a 
non-Markovian process of relaxation occurring during capture. The height of the new 
evolving peak remains almost constant for all filling times. Further, free energy changes 
seem to involve an entropy-like term since the emission time constant is nearly temperature 
independent for a particular filling time. All these features suggest that the relaxation 
mechanism follows a hierarchically constrained dynamics as proposed for glassy systems. 
A closely similar defect relaxation mechanism has been invoked to explain capacitance 
transient data in hydrogenated amorphous silicon though it is more involved in that case 
due to presence of broad density of states. We discuss plausible scenario in which such a 
hierarchical mechanism may be operative in our case. 

In Chapter 7 we have discussed the case of p-type silicon under irradiation condi- 
tions (MeV Ar+ and Au"*" ions) similar to n-type. C-V characteristics showed similar 
features as n-type though the increase in depletion width due to damage was to a lesser 
extent. Conventional thermally stimulated capacitance (TSCAP) and constant capaci- 
tance TSCAP studies reveal two trap levels in irradiated samples : one is due to^well 
known CiOi (carbon interstitial-oxygen interstitial pair) complex related hole trap and 
other is damage related trap newly found by us. Conventional DLTS studies showed three 
majority carrier trap related peaks and one minority carrier related peak. Major peak in 
this case is due to damage related trap and the DLTS peak shape is broader than peak 
expected from exponential transients. Trap emission time constants, capture cross-section 
and trap concentration are determined using constant capacitance TATS analysis. The 
major trap is located at E„-f-0.52 eV having a large capture cross-section of 9.9x10“^^ 



cm2. 

Carrier capture rate at this major trap is seen to be surprisingly slow and thermally 
activated. From DLTS-like processing of capture transients as well as from variable filling 
pulse width TATS measurements, capture barrier was estimated to be 0.66 eV which is 
quite large compared to the emission barrier mentioned above. A part of this barrier is 
associated with a macroscopic barrier for holes to reach the trap levels at the damaged 
layer. With a high macroscopic barrier, a mild inversion is the most probable source of 
filling of minority carrier traps which show up as negative peak in CC-TATS spectra. The 
time constant of the major negative peak matches with that of damage related majority 
carrier peak in n-Si. Thus the dominant majority carrier trap in n-type and dominant 
minority carrier trap in p-type Si is of same origin and it acts as an efficient recombination 
center in heavy ion damaged Si. The majority carrier peaks in p-type also seem to have 
origin in interstitials and their clusters. 

High temperature furnace annealing of damaged p-Si upto 600°C showed appearance of 
new traps in DLTS and TSCAP spectra. In contrast to the case of n-Si, large concentration 
of secondary defects are found in p-Si on annealing and their characterization would 
require more systematic studies than envisaged in this work. 

In Chapter 8 summary and conclusions of this work are collated highlighting their 
implications for study of MeV heavy ion induced defects, and defect dominated semicon- 
ductors in general. At the end of this chapter, we suggest directions in which this work 
can be extended. 

Part of this work has been communicated to various journals and presented in confer- 


ences as listed next. 
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related question of how is one able to change the interface between trap-dominated and 
trap free region through irradiation dose, annealing and even long duration storage. 

An} residual channeling in spite of misorientation of the surface with respect to the 
beam, cannot lead to observation of huge concentration of traps so far away from the 
location of the nuclear damage. We have confirmed this in another set of experiments 
[135], not included in this thesis, in which n-Si samples were implanted with very heavy 
dose of 1.5 MeV He ions at an angle as large as 60° to the surface normal so as to keep 
range of these light ions to approximately 1-1.5/im. In these experiments also the total 
depletion width turned out to be larger than that predicted from TRIM simulations. 

We propose that the dominant compensating deep acceptor traps have their origin 
in species which can migrate over large distances from the source of their creation even 
at room temperature. As migration occurs, these species would tend to form clusters 
rendering them immobile. This process of migration and cluster /complex formation can 
indeed lead to a sharp profile moving deeper into the sample. 

W^e go on to show that the pieces of supporting evidence for such a mechanism is 
increasing in some recently reported experimental and theoretical studies. Recently Priv- 
itera et al [51] demonstrated room temperature migration of Si interstitials over several 
microns. In these experiments interstitials generated at the surface migrated deeper into 
the sample to interact with preexisting damage induced defects.^ Similarly, Larsen et al. 
[50], in study of migration of Si seif interstitials at room temperature using spreading 
resistance technique, detected dopant deactivation upto depth of several microns beyond 
the region directly modified by the ions. High fluence MeV implants, specifically with 
heavy ions, have usually been reported to result in very wide spatial distribution of defects 
(2-4;im) in contrast to TRIM predictions of FWHM of only a fraction of micron [51, 141]. 

Jaraiz et al. [52] have recently reported a set of significant results from molecular 

'■Note that these experiments ■were done with low fluences and low energy ions and clustering is not 
expected to limit the range of migration in this case. 
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Figure 5.3: Arrhenius plot for peak V 2 obtained from DLTS spectra of Ar+ and Au+ irradiated 
n-Si. 

has also been detected [20]. However, we do not see any detectable density of this level 
as production density is smaller for heavier mass [20]. As the density of V 2 is smaller, 
nonexponentiality corresponding to its emission in transients is insignificant. This enables 
standard use of Arrhenius plot for estimation of trap energy and capture cross-section. 
Arrhenius plot shown in Fig. 5.3 yields activation energy and capture cross-section of Ec- 
0.43 eV and 4 x 10 cm^ respectively for V 2 level. These results are in good agreement 
with the reported values in the literature [20, 21, 23]. 
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situations where faster processes progressively appear at later times arc highly unlikely to 
occur. Our theoretical simulation shows that it’s possible to get narrow peaks in TATS 
due to simultaneous trapping anc detrapping at the depletion edge with inclusion of a 
time dependent capture term in rate equation. However, it is difficult to visualize such a 
mechanism in our case where complication occurs due to large density of localized traps. 
Without proper knowledge of the nature of nonexponentiality, calculation of the activation 
energy from these apparent peak position may involve substantial error. 

Analysis of nonexponential transients, following an external perturbation, has been 
an exciting field of research in semiconductors due to it’s connection Avith fundamental 
processes occurring in the material. In chapter 2 we have discussed many possible cases 
of nonexponentiality observed in physical systems where they have given rise to peaks 
in TATS or DLTS broader than exponential case. However, we have a distinct case here 
where peaks are narrower than exponential case. 

In Fig. 5.8(a) we have plotted two constant capacitance voltage transient in logarithm 
of time scale and corresponding CC-TATS spectra are shown in Fig. 5.8(b). The sudden 
flatness of the transient towards the end of the transient can be noted from the figure. 
In analogy with DLTS analysis, in TATS signal construction when the time windows 
is chosen such that ti is on falling part of the transient and t 2 on the flat part of the 
transient, a peak skewed at right side will occur as shown in bottom figure. Considering 
the apparent peak as the true peak, an exponential peak is fitted with dotted line in 
Fig. 5.8(b) and the broken line corresponds to the exponential fit only considering the 
left side of the peak. The dotted line shows that the experimental peaks are considerably 
narrower than corresponding exponential transient, especially for 160‘’C annealed samples. 
As the apparent peak is an artefact of premature termination of the transient as already 
seen in DLTS spectrum (Fig. 5.2), it is only the later fitting by broken line which is 
meaningful. Note that for low temperature annealed sample, D1 peak height is more and 
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more defects are activated with time and defect seem to f)e rriaxint; to der.p,.r and stable 
states even at room temperature. The V 2 related peak ha.^ not changed ;t,s they are 
thermally stable upto 350°C. Prior to these measurements, the samples were kept at 
room temperature for few months before devices were made on the irradiated wafers. 
Hence, room temperature annealing may not be alone rc.spoii.sihle fui this behaviour and 
it could have been induced by some intended process cycling during measurements. In 


earlier studies with alpha irradiated n-Si, room temperature transformation of the defect 
are reported [27], though no specific reasons have been iissigned to it. In a recent .study of 
defect stability on ion damaged n-Si using TV measurement, it wtis r«»{>orted that typical 
decay time of defects lies within 2-3 days after implantation [TI8] . 

Effect of relatively low temperature oven annealing at i60‘’C for 30 mirmtes is shown 
as curve 3 in Fig. 5.13. In this case spectra from as-implanted sample anti annealed sample 
are different in many respects. Firstly, the V 2 related peak heigiit seem to have reduced 
with annealing and Dl peak is shifted toward higher temperattirc^s. 'riii.s is attributed to 
relaxation of the damaged layer where defects Dl find a more stable l•oufiguratil(n upon 
such low temperature annealing. The broad feature Hi do(?s not seem tt) b<‘ entirely due 
to minority carrier related trap. 

Effect of further annealing at 160°C for the irradiated sample is shown in curve 4 of 
the same figure. It shows increased peak (Dl) height with small shift in peak position. 
Note that this peak is substantially narrow compared to the as-implanted case (curve 
1 of Fig. 5.13). We have discussed the cause of such narrow peaks in DLTS or TATS 
spectra. Minority carrier related peak height is increased in longer duration annealed 
sample. Thus low temperature annealing of the damaged layer shows more electrically 
active defects, and the emission energy of the dominant defect goes deeper into the gap 
with annealing. Note that lower dose implanted samples have larger concentration and 
similar defect parameters as 160°C annealed samples. The gradual change in energy with 
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Chapter 6 

Trapping kinetics of defect states in 
n-Si 

6.1 Introduction 

In the previous chapter we have shown the presence of discrete trap levels of which diva- 
cancy and the major damage induced traps are dominant and that the latter occurs in 
large concentrations with energy approximately at the middle of the bandgap. 

In this chapter, we use variable pulse-width filling experiments to monitor occupancy 
changes <is is normally done to obtain a direct measure of capture cross-sections. However, 
due to the special characteristics of these materials, which we discuss later in the chap- 
ter, such measurements yield some interesting new physics of defects instead of providing 
a routine method of obtaining capture cross-sections. Significantly, these experiments 
shed light on a long debated question of how does charge relaxation occur in disorder 
semiconductors. This has been of paramount concern in amorphous and glassy systems. 
We avoid amorphous regime of silicon, by choosing ion doses smaller than the threshold 
dose, as has been mentioned earlier. Our choice of implantation parameters yield em- 
bedded damaged regions which have been referred to in the literature as ’preamorphized 
silicon’ h Therefore, the contents of this chapter can be viewed as an investigation into 
mechanism of charge relaxation in preamorphized silicon, which turns out to be a cleaner 

” ^The term ’preamorphized’ has been used in the literature in two distinct senses. Some use it to 
denote that a sample is amorphized prior to further implantations, or is ’amorphous to start with’. Other 
use it to mean that the induced damage is heavy but below amorphization limit. We have used this term 
in the second sense meaning ’damage below amorphization threshold . 
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where nn is the occupied states for i-th level with e, and its emission rate and capture 
rate constants and Nn its total concentration. We have numerically solved the above set 
of equations ith m — 3 using model parameters for N^i and Cj and experimental emission 
time constants. For three independent trap levels with roughly equal concentrations, the 
occupancy behaviour is shown in Fig. 6.5(a). 

A distinct feature of the independent model is that the final equilibrium concentration 
of the fastest trap is not as low as seen in experiments. We have checked that this feature 
is not merely dependent on the choice of capture and emission parameters, but is integral 
characteristics of these set of equations. This suggests that the fast peak (PlA) and the 
slowest peak (Dl) may not be independent. Hence we modify the above model to include 
the assumption that PlA and Dl are two different configurations of the same defect. This 
additional constraint modifies the equations above as follows ; 


dTij'i 

dt 


— —eiTixi + Cin(iVj’i + Nxz — titi — uj’s) 


dnT2 

dt 


—62^X2 + C2n{NT2 ~ ’^r2) 


J 

■ - — = —ezTiTz + Czn{NTi + A^rs — 'nri ~ ’^tz) 
dt 

with the constraint on number of electrons as. 


( 6 . 3 ) 

( 6 . 4 ) 

( 6 . 5 ) 


n = Nti + Nt2 ~b A^rs ~ ~ ^t2 • 


( 6 . 6 ) 


Figure 6.5(b) shows the occupancy beha\'iour in such case assuming the same emission 
and capture time constant as in the earlier case. Here, we have assumed that Nt^+Ntz = 
0A5Nj' and ~ O-SSiV^. 

This indeed reproduce features close to experimental occupancy behaviour, particu- 
larly for reduction of occupancy of PlA. Here there is substantiaUy large reduction of 
occupancy of PlA compared to the earlier case and this resembles more the experimental 
obseryation. Thus PlA peak can be considered as unstable configuration of Dl related 
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As niintioned in the beginning of this chapter, relaxation phenomena with similar 
characteristic features have been reported through a series of papers by Cohen and cowork- 
ers [159. 165] using junction capacitance measurements and modulated photocurrents in 
lighth n-t\pe and intrinsic amorphous silicon. Since a broad density of states is clearly 
invoh ed in the case of amorphous silicon, the corresponding capacitance transients are 
highly nonexponential [159] in contrast to our case. The defect spectra in preamorphized 
damage Si is dominated by discrete density of states making observation and analysis 
simpler. They report temperature independence of relaxation time for initial few orders 
of magnitude in filling time and usual activation once time constant ceases to change with 
further increase in t/. In our case, the filling time needed to go over to the activation 
regime is very large and experimental limitations prevents access to this regime at con- 
venient temperature of investigation. Su and Farmer [160] had sought to explain the 
experimental features in observations of Cohen and coworkers by charge redistribution 
among a broad density of states. They reproduced the features of progressive deepening 
of emission energj', nonexponential decay and temperature independence of r for initial 
filling time t/. However, they had to invoke a very specific form (exponential) of density 
of states to fit the data. In our case, we have got rid of broad density of states by choosing 
the fluence to be below amorphization threshold and further have observed near exponen- 
tial decay spectroscopically. Hence there is no room for interpretation involving charge 
redistribution in this case. The fact that charge redistribution would occur in presence of 
large concentration of multiple states vying for capture of small number of carriers has 
been demonstrated in separate temperature regime in our samples. 

Though the mechanism which drives this novel relaxation is not clear, a plausible 
physical scenario as visualized by Cohen and coworkers [159] can be described. With 
capture of a carrier a positive polarization cloud in the vicinity of the defect deepens the 
potential well and emission is monitored from this deepened state in absence of capture. 
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Figure 7.5: Arrhenius plot for peaks (PI, P2 and P3) obtained from DLTS spectra. For major 
peak (P2) data from CV-TATS is also shown for comparison. 

in this case, the change in capacitance is too large for conventional (constant voltage) 
capacitance transient to be exponential. It is known that in case there exists a source 
of high series resistance in the sample, DLTS peak shifts in temperature and displays 
narrower width along with reduction in peak height. Further, lattice strain in the dam- 
aged region may be responsible for broadening of energy level. As argued previously, we 
rely on proper lineshape analysis using TATS, with constant capacitance implementation 
overcoming nonexponentiality due to high trap density and high series resistance. 

The feature marked P4 in Kg. 7.4 shows a signal in the positive direction, indicating 
possible emission &om minority carrier traps. Its presence also contributes to line shape 
distortion to the majority carrier peaks (P2 and P3). Definite conclusions regarding traps 
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Figure 7.12: Arrhenius plot of the capture time constants of the major trap obtained from 
CC-TATS measurements at different temperatures. 

7.4 Minority Carrier Related Traps 

In Fig. 7.4 we have shown the presence of a negative peak P4 along with conventional 
majority carrier related positive peaks. Occurrence of the negative peak (P4) in the CC- 
TATS spectra cannot be ascribed to the series resistance of the damaged layer created 
bv /Ir"^ ions since a high dc resistance does not distort transients in the CC mode of 
operation. A careful evaluation of tune constants of this negative peak shows that this 
trap corresponds to major trap observed in n-type Si. Hence, the major damage related 
trap in both n-type and p-type Si seems to be of common origin acting as an efficient 

recombination center. 

It is somewhat unusual to observe minority carrier traps using Schottky diode. Under 
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annealed samples have been studied to monitor progressive changes in electrically active 
defects. 

In addition to conventional techniques of steady-state capacitance- voltage (C-V), deep 
level transient spectroscopy (DLTS) and thermally stimulated capacitance (TSCAP), 
novel methods of transient based isothermal spectroscopic techniques such as time an- 
alyzed transient spectroscopy (TATS) have been extensively utilized. Implementation of 
constant capacitance (CC) mode of operation in spectroscopic characterization overcomes 
interpretational difficulties. Complimentary methods of electrical of characterization such 
as impedance analysis and current-voltage (I-V) measurements have also been used in this 
study. 

The major contributions of the present work are summarized below. 

A. Methodological Contributions 

A proper appreciation of limitations of standard use of defect characterization tech- 
niques and various factors affecting capacitance measurements have been central to this 
work. Common origin of nonexponentialities in capacitance transients and their effect 
on trap parameter evaluation from conventional deep level spectroscopies have been as- 
sessed critically. Spectroscopic lineshape analysis of transient signals in the study of defect 
phenomena plays important role and determines suitability of spectroscopic techniques. 
Specifically, use of time analyzed transient spectroscopy (TATS), an isothermal spectro- 
scopic technique, is shown to provide cleaner interpretations, and better insights into 
nonexponential transients which indeed manifest in various forms in our experiments. A 
brief case study involving analysis of experimental data of photo-induced current tran- 
sients in semi-insulating gallium arsenide (SI-GaAs) has been presented to demonstrate 
superiority of TATS in analyzing nonexponential transients. Our work shows that TATS 
should be the preferred method in the study of physics of defects in defect dominated 
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Chapter 1 
Introduction 


1.1 Preamble 

Ion beams as instrument for modification and analysis of semiconductors have contributed 
in large measure to the remarkable progress during leist decades in semiconductors. The 
success of ion implantation as an industrial process are due to its many advantages which 
include precision, speed and reproducibility. This is achieved through accurate control of 
energ}’, dose and species of ion beams. It has become indispensable in standard device 
technology mainly as a method of introducing dopants at desired depths, at temperatures 
lower than that needed for diffusion and with less stringent requirements on purity of the 
source. The application of ion beams in silicon processing and characterization has been 
recently reviewed [1]. One of the primary disadvantages of the method, however, is that 
it creates radiation damage requiring high temperature annealing to repair the crystal and 
to activate the dopants by moving them to substitutional sites. The need to understand 
the nature of damage and defects created in the implantation process and their evolution 
under processing conditions have spawned a vast literature over many decades [2]. The 
remarkable progress in this field has not only helped in optimization of process, but also 
has enriched fundamental understanding of a whole range of solid state phenomena [3]. 
The technological benefits of this knowledge are now being reaped in increasingly large 
number of applications involving ion beams to an extent that the field has come to be 
recognized as ’defect engineering’. 

However, most of these developments have centered around applications of low to 
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medium energ\' ions, or light particles such as electron, protons and neutrons nniiiil} to 
understand the nature of radiation induced defects. It is only recent!} that therc^ has 
been an increasing interest in the use of heavy ions in the MeV eiierg}" range ciri\en 
primarily by the technological demands of ability to modify semiconductor materials at 
distances larger than several microns below the surface. The problems and issues arising 
out of use of high energy heavy ions are in many ways different from that of their ke\ 
counterparts. The present work is an attempt to contribute to this new emerging area 
of studies on defect phenomena due to MeV heavy ion implantation primarily through 
electrical studies. We postpone a more comprehensive statement of our objectives and 
scope to a later section in this chapter. We first give a brief review of issues in Me\ ion 
implantation with primary focus on defect related studies of relevance to this work. 

1.2 Technological Significance of MeV Ion Implan- 
tation 

Though there are several potential applications of MeV ion beams, they have mostly been 
used in three distinct areas : (i) deep implantation of dopants to fabricate retrograde wells 
or tubs, (ii) production of buried dielectric layers such as silicon dioxide in silicon or semi- 
insulating isolation in III-V compounds, (iii) introduction of controlled gettering sites to 
remove harmful impurities such as iron [1]. Of these, MeV technology is already in use for 
over a decade in the fabrication of retrograde wells for DRAM, CCD and CMOS applica- 
tions [4, 5] which would otherwise need lengthy high temperature dopant diffusion. This 
has led to the increasing importance of studies on defects and the interaction between 
multiple implants in new processing protocols. Production of buried dielectric layers has 
focused on implantation of oxygen in silicon and the efforts has been towards reduction of 
defect density levels to acceptable limits. Similar techniques of isolation is being adopted 
in III-V processing by implanting impurities or engineering defects which led to deep 
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compensated layers [6]. These applications typically need high energy and high dose 
implants. Effective gettering of detrimental impurities have been demonstrated by using 
intrinsic implant induced damage, nanocavities and high energy carbon implants. Com- 
mercial exploitation of these gettering techniques is imminent if the process is understood 
in sufficient detail [1]. 

Apart from these, there are several emerging applications involving mechanical and 
optical patterning using ion microbeams to form optical waveguides and cavities and selec- 
tive disordering of superlattices [7]. Ion beam induced synthesis of desirable compounds 
such as sillicides and germanides using heavy ion implantation has also shown in potential 
contact technologies [8]. MeV heavy ion implantation has been shown to be an effec- 
tive way of obtaining pure amorphous layers which is essential for technologies involving 
solid phase epitaxy. The growth of applications has been stimulated by availability of 
purpose-built MeV ion accelerators. 

There are many additional advantages that accrue with the use of MeV implantation 
technology in processing. Some of the important ones include : (i) increased yield and 
shorter design cycle time due to simplification of process steps, (ii) reduction in thermal 
budget, (iii) avoidance of warping in large diameter wafers owing to low temperature 
processing. In DRAM fabrication technology, improvements have been reported in several 
device characteristics such as lower offset error rates, latch-up immunity, increased RAM 
retention time, improved low voltage operation etc. 

With so many emerging applications, high energy implants are likely to become an 
integral part of semiconductor processing technologies, just as low energy implanters are 
standard industrial and laboratory tools today [9]. The key to the success of these 
applications and their widespread adoption hinges on our ability to predict process effects 
through a detailed understanding of defect related issues in each case. We can illustrate 
the importance of defect issues by citing another example apart from those mentioned 
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above. There has been a concerted effort over the recent years to introduce raie earth into 
silicon^ and other semiconductors using MeV ions to take advantage of optica! transitions 
in rare earth atoms within a semiconductor host [10, 11]. There has been a recent 
breakthrough in obtaining room temperature electroluminiscence at telecommunication 
wavelength of 1.5 jivci using MeV Er and O co-implantation in silicon [12]. It is believed 
that the mechanism of excitation is mediated by deep levels defects, though the defect 
issues arising out of these implantations are little understood at present [13]. 


1.3 MeV Implantation : Energy Deposition and Range 


The production of range and distribution of implanted ions is based on knowledge of 
processes through which the projectile looses energy in the medium. The two most im- 
portant mechanisms of energy loss are [14] : (i) electronic stopping, in which ion looses 
energy by inelastic collisions with bound electrons in solid, and (ii) nuclear stopping, in 
which part of the kinetic energy is transferred to a nuclei in near elastic collision. The 
relative importance of the two stopping mechanisms depends on the energy and mass of 
the implanted ions, and the atomic density of the solid. The cross-section of stopping S 
associated with each mechanism is represented as 


Se,n — 


'iV 


dE 


dx 


eV.cm ^ 


( 1 . 1 ) 


-1 e,n 


where ^ is the energy loss per unit distance for either electronic or nuclear stopping and 
N is the atomic density of solid. Figure 1.1 shows a schematic of these stopping power 
as a function of energy [14]. For low energy ions typically used for shallow junction 
formation, the nuclear stopping is the dominant mechanism. The nuclear stopping power 
is small for low energies since electronic screening reduces the effective atomic number 
of target nuclei, and it becomes again small at very high energies since interaction time 


^Rare earth doped semiconductors has become a field in itself at present. 
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Figure 1.2: Vacancy and ion profile for 1.45 MeV Ar"^ ions (dose 5x10^^ cm in Si using 
TRIM simulation. 

approximately the range of the ion divided by its velocity. The thermal spikes created by 
ionization and excitation along ion track decays away by ~ 10“^^ s after entry of the ions, 
and many of the vacancies also created along the ion track and the tracks of the cascading 
secondaries recombine in ~ 10“® secs. Light ions tend to leave tracks characterized by 
relatively small amount of damage. They slow down by electronic stopping process ini- 
tially with little displacement damage. Eventually nuclear stopping becomes predominant 
at the end of their range where they undergo collisions involving large scattering angles. 
Thus there is little lattice damage to the crystal except near the end of the ion range. 

Heavy ions by contrast may create damage clusters along their path, as shown in 
Fig. IfigS. These ions undergo a relatively higher degree of nuclear stopping than the light 
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Figure 1.3; Schematic of the ion track in solid and associated damage for a light ion (top) and 
a heavy ion (bottom) (after ref.[6]). 


ions, displacing target atoms right from the surface inwards. Large density of collision 
cascades give rise to considerable lattice damage within relatively small volume. For both 
light and heavy ion damage, the volume of the crystal in which the ion energy is deposited 
is usually larger than the volume in which lattice damage occurs. Much of the ion energy 
is transferred to the lattice as a localized thermal spike. 

There has been a large number of studies dealing with implantation of dopants at 
JvIeV energies. The focus of these studies were on the issue of electrical activation of these 
dopants and damage evolution upon annealing. Here again resistivity, secondary ion mass 
spectrometry (SIMS) and transmission electron microscopy (TEM) studies have been the 
main tools. The annealing studies get often complicated by the fact that the implanted 
ions are substitutional electrically active dopants. Also in-situ dynamic annealing have 
been found to influence defect introduction. Even a reverse temperature dependence of 

i 

annealing has been reported at high dose rates due to overlapping cascades in time [18]. 
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The role of damage induced purely by ballistic processes and species dependant effects of 
damage appear simultaneous in each studies. The effects are due only to damage induced 
by massive projectiles can be separated from species dependent effects using noble gas 
ions which are ultimately electrically inactive in silicon. However, studies in this direction 
has been rare. 

The nature of final disorder depends on a variety of factors such as ion mass, energy, 
dose, dose rate, temperature of implantation. For light particles and low fiuences only- 
point defects are produced while for high fiuences of heavy particles one encounters the 
other extreme of amorphous layers. These are the two extreme regimes of MeV implants 
that have been intensively studied in recent years though with different methods. The 
point defect regime has been studied mostly using electrically sensitive techniques and 
magnetic resonance methods, whereas heavy damage and their evolution has been studied 
using structural tools such as cross-sectional transmission electron microscopy (XTEM). 
It is the middle regime of the damage, below amorphization but much above point defect 
regime, that is least studied. In fact, a recent study christens this regime of damage as 
’missing link’ [19]. We review recent results in each of these regimes for MeV implantation. 

1.4.1 Point defects 

The regime of low damage, specially with light particles such as electrons, neutrons and 
protons, only point defects seem to be generated. These defects have been studied in 
great detail in the past and it would be impossible to review this vast area within the 
scope of this work. It is only recently that there have been studies on similar point defects 
generated by heavy ions in MeV range by keeping the fluence of the ion very small (< 10® 
cm“^) [20]. An impressive list of many defects have been studied using spectroscopic 
techniques such as infrared (IR), electron paramagnetic resonance (EPR) and transient 
based spectroscopies such as deep level transient spectroscopy (DLTS). In this regime of 
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defects, identification of electrical signature of the defects have played a central role in 
determining their nature and properties. The information is generally sought in terms of 
trapping and detrapping parameters such as energy, capture cross-section and relaxation 
mechanisms, and mechanisms of formation, diffusion, and interaction with other impuri- 
ties. The actual defects produced in any given experimental situation depends critically 
on a number of parameters such as implant temperatures, ion species and energy, dose, 
dose rate, conductivity type, impurity present, annealing conditions etc. Such variability 
has contributed to continuing interest in this field. We will limit ourselves only to those 
well known aspects which are relevant to this work. 

The hierarchy of stable defects generated from primary defects such as frenkel pairs, 
vacancies and interstitials have been most thoroughly studied for electron irradiation. 
These intrinsic defects pair among themselves or with commonly found impurities such 
as carbon, oxygen or with dopants. The most of these simple point defects pair anneal at 
well characterized temperature below 450°C. Some of the well known interstitial complexes 
are Cs — Ci, Ci — Oi, Bi — Bg, while vacancy complexes include divacancy (I'a), V-0, V- 
P. Many different interaction mechanisms between them have been proposed to explain 
relative concentration of these defects under different experimental situations. Though a 
wealth of data has been generated, a coherent picture of mechanisms is yet to emerge. 
Most commonly detected defects in ion irradiated n-Si are oxygen vacancy pair (V-0) 
or A center, two different charge state of the divacancy {V 2 ~ and Vf) and phosphorus 
vacancy (P-V) or E-center. A summary of energy and capture cross-section of reported 
defect levels in as-implanted n-Si is presented in Table 1.1. 

In p-type Si, the assignment of defect identification has been less straightforward. This 
difficulty may be attributed to a general lack of consensus in the literature regarding the 
identity of similar defect levels observed in electron irradiated p-Si. Hence it is difficult to 
tabulate defect levels for p-Si. Nevertheless, divacancy and interstitial-carbon-interstitial- 
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Table 1.1: Comparison of reported defect levels due to various ion in n-Si 



Ions 

o.v 

Vi- 

Vs/P.V 

Unknown(?] 

Wang [28] 

H+,He+,B+,P+ 

E(0.18) 

E(0.23) 

E(0.41) 


Troxell [21] 

~H+,B+,Si+ 

E(0.16) 

E(0.23) 

*3 X 10-^® 

■MU 


Hallen [22] 

J?+,iTe+ 

E(0.18) 

23 X 10-^5 

E(0.24) 

2.6 X 10-1® 

E(0.42) 

1.7 X 10-1® 

E(0.32), E(0.4J 
E(0.5) 

Svensson [20] 

H+,He^+,0^+ 
S'^'^ , Br^-^ , 


E(0.23) 

BSIM 

E(0.35) 

Palmetshofer [23] 

H+,D+ 

He+ 

E(0.17) 

8 X 10-^5 

■H 

HiB 

E(0.30), E(0.35 
E(0.39) 

Indusekhax [24] 


E(0.18) 

1.2 X 10-^® 


E(0.48) 

5.3 X 10-P 

E(0.28) 

E)(0.51) 

Krynicki [25] 

C 

... 


E(0.41) 

E(0.63) 

Mary [26] 

u 

E(0.17) 

2 X 10-1^ 

E(0.22) 

7 X 10-1® 

nil 


Zafar [27] 

He^+ 

E(0.17) 

2 X 10-^3 

E(0.20) 

2.2 X 10-1® 

E(0.43) 

3.4 X 10-1® 

E(0.26),E(0.3i: 
1.7 X lO-i'i 

Benton [19] 

Si+ 

E(0.18) 

E(0.23) 

E(0.40,0.44) 



* Capture cross-section in cm^ 


oxygen {CiOi) have been observed in most of the works [21, 28] with hole traps at 0.23 
eV and 0.36 eV respectively above valence band. Other interstitial related defect such as 
B,Cs at Eu -1-0.26 eV and clusters are sporadically noted in the literature [29]. 

1.4.1. 1 Divacancy 

Among radiation induced defects in silicon, divacancy (V 2 ) is the most ubiquitous and 
therefore is also one of the most well studied centers. Unlike monovacancy, the divacancy 
is stable and immobile at room temperature. As far as electrical activity is concerned, di- 
vacancy introduce deep levels corresponding to multiple charge states ranging from singly 
positive to doubly negative in heavily doped n-type Si. The level scheme for divacancy 
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Table 1.2: Charge state transition and energj' level scheme for divacancy trap in silicon 


Charge state 

Transition 

Encrg}' level 


emission of electron from neutral state 

E„+0.25 e\' 


electron emission to become neutral 

Ec-0.42 e\' 

v\-i- 

electron emission to become singly charged 

Ec-0.23 e\‘ 


that has emerged from several DLTS studies are given in Table 1.2 [20]. 

There have been controversies regarding whether Ec-0.23 eV level and Ec-0.42 e\' level 
are indeed the two charge states of the same defect ‘or not. The concentration of the two 
levels do not show one to one correspondence in many cases. For example, irradiations 
using h'leV proton and alpha particles in n-Si showed difference in dose dependence of 
the two levels [22]. Recently, Svensson et al. studied the two levels in n-Si irradiated 
with a series of high energy’’ particles starting from light particles to very heavy ions 
(electron, 'H+, "He2+, , 323 ^+ , and [20], Though the concentration 

of the two levels are same in electron irradiated material, their ratio was dependent on 
depth in case of ion bombardment. From a study of depth profile it was concluded that 
the concentration of level corresponding to Ec-0.42 eV is significantly larger than that of 
Ec-0.23 eV level in regions of large damage most probably due to strain favouring one 
state over the other. The Ec-0.42 eV peak in DLTS studies seemed to be broader for 
increasing projectile mass. Observations such as these and other properties have been 
sought to be explained on the basis of a model which invokes motional averaging of the 
two states at higher temperatures. The absence of Ec-0.23 eV level in MeV heavy ion 
implants is attributed to hindrance in reorientation among states at low temperature due 
to internal lattice strain associated with damage [20]. However, unresolved issues in 
terms of relative concentration of these two oft-observed levels persist and require more 
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studies by complimentary techniques. The divacancy is paramagnetic in the two charge 
state and Vj. There has been several thorough studies using EPR, IR and ENDOR 
spectroscopies by several workers [30, 31]. The defect has a point group symmetry Da,/ 
with four equivalent orientation in < 111 > directions. However, in its ’+’, ’0’ and 
’ charge states, presence of Jann-Teller distortion lowers the symmetry to C 2 / 1 , a fact 
that has been confirmed by EPR observations below 20K. For higher temperatures, the 
defect does not relax to lower symmetry configurations, rather a motionally averaged 
Dad state appears making V 2 a six-silicon atom center. This electronic bond switching 
was first observed by Watkins and Corrbett [32] and subsequently confirmed by Sieverts 
and coworkers [30, 31]. The symmetry of V 2 state should be Dad since no Jann-Teller 
distortion is expected. However this symmetry state is not observable since it is EPR 
inactive. 

The formation mechanism of divacancy has also been controversial. Pairing of single 
vacancies due to bombardment in neighboring sites has been considered improbable. A 
more likely mechanism is that two single vacancies, created directly through primary 
knock-ons or indirectly in recoiling atom collisions migrate and subsequently agglomerate 
to form the more stable V 2 state. These defects anneal in the temperature range 300-350‘’C 
in low oxygen content silicon. Involvement of oxygen has been suspected in annealing 
kinetics by comparison of results for float zone and Chrochralski silicon. 

The observed divacancy concentration in MeV implanted samples is much smaller than 
that would be predicted by a simulational package such as TRIM because of dynamic 
annealing and recombination. It has also been reported that divacancy concentration 
saturates for high dose implantations, the saturation value being one order of magnitude 
larger close to the surface compared to the deeper region (> l/xm for MeV implant) 
[33]. Interestingly, no divacancy centers are detected through EPR measurements in 
amorphous silicon [34]. The reasons for observed low concentrations of divacancy in high 
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fluence implantations needs further systematic studies. 

1.4. 1.2 Other vacancy associated defects 

Vacancies produced in silicon during a room temperature irradiation are mobile and have 
been shown to combine with dopant atoms [35]. This combination appears to be the' result 
of a Coulombic interaction. Formation of such a pair is more efficient in n-Si compared 
to p-Si. 

The vacancy-phosphorus (P.V) center or E-center is a well understood defect and has 
an acceptor level 0.44 eV below the conduction level [36] and have capture c:ros.s-section 
~ 10“^^ cw?. It has been observed to form even in p-type doubly doped silicon. The 
E-center has been shown to have a charge-state dependent annealing [37]. In DLTS 
spectra, E-center related peak is overlapped by the presence of V-f related peak fts they 
have similar trap parameters. The overlapping nature of the two peaks in DLTS has been 
a source of controversy as regards the relative concentration of the two levels. However, 
differences in annealing behavior is used to differentiate between E-center and divacancy 
peak. E-centers are annealed out at ~ 150®C whereas V-z centers are stable upto SOO^C. In 
moderately doped (phosphorus) Si, E-center concentration is reported to be only < 20% 
of the di vacancy concentration in low dose MeV ion implanted n-Si [18]. However, in 
GeV heavy ion implanted Si, presence of E-center in full concentration has been reported 
by Mary et al. [26]. 

The vacancy - oxygen center (A-center) is a fairly well understood defect [36, 38] with 
an acceptor level at Ec — 0.17 eV. It occurs even in low oxygen content silicon. This defect 
is formed when a \’acancy is trapped by interstitial oxygen atom. This trap level exhibits 
a very large capture cross-section for electrons, cr„ ~ 1 x 10“^“^ cm^ [36]. The oxygen - 
vacancy pair has been reported to be stable upto ~ 350°C [36]. 

A large number of oxygen plus multivacancy defects have been observed and studied 
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1.4.2 Defect Migration and Clustering 

Studies on diffusion of impurities in semiconductors has been of central importance for 
both fundamental and technological reasons. This has led to a wealth of data on high 
temperature diffusivities and various mechanisms of diffusion [48]. Without attempting 
to review these, we only focus on the recent interest in migration and agglomeration of 
intrinsic point defects such as silicon self-interstitials and vacancies. 

Since point defects in crystalline silicon generally occur in low equilibrium concentra- 
tions, measurements of their diffusive properties have been fraught with uncertainties. 
For example, high temperature diffusion data show very large discrepancies from different 
experiments. This has been illustrated by Gosele and Tan [49] by collating information 
from several authors. Extrapolation of defect diffusivities determined from high tem- 
perature measurement, to room temperature shows that vacancies should show no room 
temperature migration (with diffusivities < cm/s). There has also been along-held 
view in the literature that defect evolution in Si under irradiation is fully determined by 
recombination and clustering vidthin the damaged region. However, as mentioned ear- 
lier, long distance migration even at cryogenic temperature have been shown from EPR 
experiments and DLTS. From these experiments room temperature diffusion coefficients 
for vacancies and interstitial have been estimated to be D/ = 3.4 x 10“'^ cm^/s and 
D„=4.2xl0“® cm^/s, respectively. These large discrepancies could be attributed to ei- 
ther change in structures at high temperature or ionization enhanced migration under 
irradiation condition. However, there are several recent studies which show long range 
room temperature migration of these defects. Larsen et al. [50] have shown the impor- 
tance of long-range trap limited diffusion for self-interstitials using spreading resistance 
measurements on heavy ion irradiated silicon. Coffa et al. [51] have provided evidence 
of room temperature migration of self-interstitials over several microns through specially 
designed experiments. The interstitial and vacancy diffusivities being used for quantita- 
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tive prediction of phenomena in simulation at present are Dj — 0.01<’xp( O.D/A'T) and 
Dy = 0.001exp(-0.43/^T) [52]. 

The ion generated interstitials and vacancies, being mobile at room temp('ranir<?. are 
most likely to end up in clusters if they escape recombination. Though presence of such 
clusters have been invoked in many studies, direct correlation between their occurrence 
and electrical activity is still lacking. The vacancy clusters have been grown in the past by 
cooling down as-grown crystals and were reported to be electrical inactive [53. 54] . Early 
observations of interstitial type cluster had come from electron microscopy in relation to 
formation of dislocations. There has been so far no easy effective experimental way of 
bridging the gap between observations of simple point defects as in EPR and DLTS on 
the one hand, and electron microscopic investigations of observable extended defects on 
the other hand. It is for this reason that intermediate defect configuration between two 
such clusters have not been studied in the literature. 

There has been recently attempts to simulate the effects of migration and clustering 
and indirectly infer their properties from matching with experiments [55, 56]. The number 
of point defects produced due to ion bombardment depends on ion species, energ}' and vary 
over large orders of magnitude per ion. A survey of literature does suggest that clustering 
is predominant when the projectiles are heavy and doses are high most probably because 
long-lived excess interstitials are available under such conditions [52]. 

The binding energies of interstitial clusters have been estimated from molecular dy- 
namics (MD) simulations [57], and from smooth fit to the data are given by £'w(n) = 
2.1 - 1.45(V^-v'n - 1) where n is the number of interstitials. Similarly, the cluster bind- 
ing energies for vacancies have been estimated to be Ebv{n) = 3.65 — 5.15[n^^^ — (n— 1)^/^]. 
Jaraiz et al. [52] have shown from MD calculations that a fraction of interstitial pop>- 
ulation remain free to migrate, annihilate and agglomerate even after three hours of 
irradiation at room temperature for a dose of 5x10^^ cm“^ of Si. The depth distributions 
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of interstitials .were shown to peak at twice the implant range which was attributed to 
momentum transfer from the implant and profile distortion due to initial annihilation 
at the surface. In annealing studies, clusters are believed to undergo Ostwald ripening 
process eventually dissolving at high temperature at the surface. The simulation in par- 
ticular predict continual formation of interstitial atomic cluster with increasing size at 
room temperature. In our view, a careful study of this regime is essential at present in 
order to enable corrections among a variety of phenomena. 

The electrical activity corresponding to these clusters is a largely unexplored area. 
Benton et al. [19] have recently associated two hole levels at 0.29 eV and 0.48 eV with 
clusters in p-type Si. They also argue that DLTS peaks correspond to cluster smaller 
than 50A. However, it is surprising that one observes very well defined point defects like 
signature in the presence of expected distribution in sizes of these clusters. 

1.4.3 Extended defects 

When ion damaged materials in the dose regime > 10^^ cm“^ are annealed, extended 
defects have been observed to form. There have been fairly large number of studies 
on extended defect formation using cross-sectional TEM as the principal instrument. 
On annealing in the range of 670-800°C following high dose and high energy implants, 
the damage is known to evolve into rod-like interstitial agglomerates large enough to 
be observed by TEM and are known as {311} defects [56]. These intermediate defect 
configurations can either dissolve giving rise to self-interstitials or provide a nucleation site 
for the generation of full blown dislocation loops. The concentration of these agglomerates 
decreases while sizes increase on continued heat treatment as in Ostwald ripening process. 

The nature of evolution of damage in this range of extended defect formation has been 
systematically categorized by Jones et al. [58] on the basis of TEM studies of a very large 
number of implanted samples under varying conditions. There has also been an effort to 
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determine a threshold dose below which extended defect formation wonhi not orrur on 
annealing [59]. This dose is believed to be below the amorphization threshold, hut the 
parameters that determine it are not so well known [60]. 

The electrical signature of extended defects in DLTS studies, remains controv(‘rsiai 
in spite of several studies especially in plastically deformed silicon. Generally, a DLTS 
signature is attributed to an extended defect if it is a broad signal corresponding to a 
many-electron center. A characteristics of such many-electron center is believed to be 
slow increase in occupancy (logarithmic in time), due to changing Coulombic barrier dur- 
ing capture. However, such increase can be due to various other reasons such as local 
macroscopic barriers or coupled kinetics in multi-trap situation. Hence careful interpre- 
tation of conventional DLTS spectra are needed in such cases [61]. In high do.se ion 
implanted Si, a commonly observed level at 0.5 eV is generally associated with an ex- 
tended defect [29]. However, systematic correlations leading to definitive identification 
is still lacking. 

Injection of interstitials from extended defects at high temperatures and issues of 
interstitial migrations and clustering are being studied extensively in relation with a 
technologically important phenomena of transient enhanced diffusion (TED) of dopant in 
Si. TED is the phenomena of large nonequilibrium of diffusion of dopants in implanted 
samples at about 800®C. The non-equilibrium diffusion vanishes at large times. It is by 
now well established [46, 62] that the anomalous diffusion axises from excess silicon self- 
interstitials that are generated from extended agglomerates, such as {311} defect, though 
consensus regarding mechanisms of diffusion is yet to emerge. A recent article by Stolk 
et al. [46] discusses many aspects of these phenomena in great details and provide an 
agenda for researches in this field. Significantly many of the issues raised deal with defect 
migration and clustering and interstitial injection kinetics. 

In this work, in a bid to study as-implanted damage , we have avoided the regime of 
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extended defects by limiting ourselves to processing temperatures below SSO^C. 


1.4.4 Defects and amorphization by irradiation 


It is well known that ion-beam induced damage for sufficiently high dose, can give rise to 
amorphization. Though it is possible to make estimates of dose requirement from simple 
considerations, experimentally amorphization tends to depend on many parameters. The 
total number of atoms displaced by an incoming ion can be estimated as [63] 


Nd ~ 


3l 

2Ed 


where En is the total energy loss of a particle in primary and secondary nuclear colli- 
sions. and Ed is the displacement energy of a lattice atom. Typically, Ed is ~15 eV in 
semiconductors. The number of displaced atoms, for example, with a 1.4 MeV Ar ion 
in Si is ~3400. Thus at doses near 10^“^ cm~^ the displacement disorder concentration is 
approximately equal to the Si atom density, leading to amorphization. 

If one assumes all target atoms must be displaced to form an amorphous layer, then 
the critical dose for its formation is [14] 

2EdN 

{dE/dx)n 

This relation always underestimates the critical dose because of self-annealing or dynamic 
annealing of vacancies and interstitials. Some attempts have been made to take this into 
account by including defect diffusion and recombination processes. The critical dose for 
amorphization depends strongly on the temperature of the implantation [64]. With an 
increase in the dose rate, the critical dose usually decreases because the damage is able 
to accumulate before it self-anneals. For very high dose rates however, the sample can 
be heated by irradiation, and dynamic annealing of the damage occurs. In this case an 
amorphous layer may never form. The critical dose for creation of an amorphous layer in 
Si for room temperature Ar+ implant has been shown to be ~ 1.2 x 10^^ cm~^ [65]. The 
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annealing of point defects and locally amorphous zones occurs by defect diff\ision and 
recombination at ele^^ted temperatures. The energy deposition incretises with nuclear 
energy loss (g)„ which in turn increases with projectile mass as in modified Kinchin- 
Pease model. The earliest model of ion beam induced amorphization w;is by Morehead 
and Crowder [66] who attributed it to sufficient overlap of the amorphized cylindrical 
cascades. This is viewed as a heterogeneous mechanism of amorphization. 

However, there has been homogeneous amorphization models, proposed initially by 
Swanson et al. [67] and more recently by Holland et al. [68], according to which amor- 
phization is a phase transition induced by sufficient number of defects in crystalline silicon. 
Motooka et al. [55] have also proposed a model of homogeneous amorphization using 
molecular dynamic calculations in which accumulation of divacancy and di-interstitial 
pairs beyond a threshold value lead to amorphization. However, empirical evidence of 
identity of accumulating defect is still lacking. 

Apart from process of amorphization, implant induced pure amorphous silicon (a-Si), 
in contrast to hydrogenated a-Si, has also aroused interest and it has been studietl by 
a variety of methods [69, 70]. It has been shown that pure a-Si undergoes structural 
relaxation due to kinetics of point defects. Both structural relaxations (i.e. short range 
ordering) in a-Si and defect annihilation in crystalline silicon is found to obey bimolecular 
relaxation kinetics. The nature of short range ordering in pure a-Si, as obtained from ion 
implantation, turns out to be different from that obtained in usual hydrogenated a-Si, 
where creation of metastable defects and their annealing dominates. 

In spite of the realization that most phenomena in implant induced a-Si can be un- 
derstood within the framework of a defect mediated structural relaxation, the nature of 
participating defects is unclear. This is mainly due to difficulties in establishing clear 
connection between structural studies and defect spectroscopic techniques. In our view, 
such correlations can be best done using proper choice of regimes of controlled disorder. 
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The choice of dose for heavy ion implantation in this work is motivated by this consider- 
ation. The doses are kept below amorphization threshold but sufficiently high to create 
conditions of large disorder in crystalline Si with the hope that it would enable possible 
identification of dominant defects. We show in a later chapter that the study of charge 
relaxation mechanisms in samples so disordered can provide significant insights into re- 
laxation mechanisms in a-Si in general. However, we postpone a review of related issues 
to a more appropriate context in Chapter 6 where these results are presented. 

1.5 Motivation and Focus of The Present Work 

The foregoing review shows that there are many unresolved issues in our understanding of 
radiation induced damage in case of high fluence, MeV heavy ion implantation in silicon, 
though earlier studies with light particle irradiation provide a relatively firm basis. It is 
not clear as to which aspects and to what extent our understanding of light particles, low 
energy and low dose implantation can be carried over to the cases with combinations of 
parameters such as high energy, heavy ion and high fluence. The inherent differences are 
primarily due to nature of disorder created due to high density collisional cascades, and 
complexities arising out of movement and agglomeration of large concentration of defects. 
There is a lack of studies on electrical manifestation of phenomena such as defect migra- 
tion and agglomeration. The empirical support to test results of recent molecular dynamic 
simulational studies [52] are inadequate. The connection between damage as studied by 
structural tools such as cross-sectional transmission electron microscopy (XTEM) and de- 
fects controlling electrical characteristics such as compensation and trapping are missing. 
This is primarily because electrical characterization has been used for limited purposes 
in this regime of parameters. In fact, studies involving electrical characterization of point 
defects deliberately avoid complexities attributable to heavy damage in order to preserve 
validity of the techniques used and simplicity of interpretations. Hence the power of spec- 
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troscopic electrical defect characterization techniques has so far been grossly underutilized 
in the study of phenomena related to damage created by hea^*}' ions. 

The principal motive of this work is to study electrical manifestation of damage created 
in deep buried layers using MeV heavy ions and relate them to the nature of dominant 
electrical defects created. It was essential therefore to critically evaluate validity ol stan- 
dard electrical methods of characterization for these cases. The sources of difficulties in 
electrical studies are many. Firstly, the concentration of defects are largo inx'alidating 
straightforward use of popular characterization methods such as DLTS. Secondly, their 
spatial distribution is highly nonuniform leading to inhomogeneous modification of the 
host. Thirdly, the damage layer of interest is located deep within the sample. Though 
transport studies would certainly be intrinsically interesting in these samples, we delib- 
erately avoid them as methods of characterization since transport parameters typically 
constitute contribution in varying degree from many different mechanisms. Hence we 
limit the scope of this study to capacitance based techniques. Characterization tech- 
niques based on depletion layer capacitance provide a degree of control unavailable in 
transport measurement in terms of defect depth from the surface and spectroscopic char- 
acter. Capacitance based techniques, specially transient based spectroscopic methods, 
have been highly successful in the study of point defects in the well characterized regime 
of radiation induced damage (as for example, in the case of electron irradiation). A run- 
ning theme of this work has been to suitably extend and modify these techniques keeping 
in view anticipated difficulties such as various sources of non-exponentiality of transients. 

With this in focus, the damaged layers of interest are created using MeV heavy ions 
(Ar'*' and Au"*") within the depletion layers of Schottky devices, which serve as laboratory 
for all defect studies in this work. The doses are kept below amorphization threshold. The 
choice of Ar"^ ion has been motivated by the need to avoid issues of electrical activation 
of the implanted ion so as to focus only on damage induced by heavy ion irrespective 



Chapter 2 

Principles of Experiments and 
Techniques of Data Analysis 

2.1 Introduction 

Defects and impurities in crystalline semiconductors introduce states in their otherwise 
forbidden gap. Substitutional dopants give rise to shallow states ( states close to band 
edges) which control carrier concentration, while deeper states control carrier lifetime. 
Hence electrical characterization of these states has been both convenient and appropri- 
ate for most applications. Among electrical measurements the use of capacitance based 
techniques has several advantages over transport measurements. One of the principal 
advantages is the wide availability of robust instrumentation for carrying out sensitive ca- 
pacitance measurements. Further, in the context of semiconductor devices, space charge 
regions in a wide variety of structures provide convenient laboratories for implementing 
capacitance based techniques to study defects in environments most suitable from ap- 
plication point of view. The use of capacitance-voltage (C-V) measurements is one of 
the most popular methods of profiling shallow dopants. The experimental study of deep 
levels in semiconductors was revolutionized following introduction of capacitance based 
methods. In fact we owe much of our present understanding of deep levels and the huge 
amount of information collected over the last two decades to transient based capacitance 
techniques. 

The electrical characterization of deep level defects is normally carried out in terms of 
phenomenological parameters of these levels provided by Shokley-Read-Hall (SRH) model 
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of trapping and recombination processes. The basic information that one hocks from 
such measurements include thermal activation energy (Ej'). carrier capture rrciKS-sectioiis 
(cr), concentration (Nr) profile etc. of traps. With the advent of transitmr spectroscopic 
methods such as deep level transient spectroscopy (DLTS) [71], easy electrical rrrre-:itiori 
of defects through their emission signatures made a large number ot dynamical studies 
possible. There are number of excellent books, monographs and review.s dcn-med to the 
subject in the hterature [72, 73, 74, 75, 76]. We confine ourselves in this chapter to brief 
discussions of only those aspects which are relevant to our study of deep levels in ion 
induced buried layers. 

In this chapter, at the outset we discuss in presence of high trap density and .series 
resistance in a Schottky diode, validity of standard C-V profiling. The effect of .series 
resistance, leakage current, high defect concentration (Nr/Nd > 1, NV is traj) concentra- 
tion and Nd is shallow doping concentration), on capacitance transients are examined in 
this chapter. Their effect on trap parameter estimation is analyzed in detail. A variety 
of nonexponential transients are encountered in the study of defects in the ion-damaged 
materials. Critical comparison between temperature scanning techniques such as DLTS 
and isothermal techniques such as time analyzed transient spectroscopy (TATS) are made 
as regards their suitablity for quantitative evaluation of nonexponentiality and to under- 
stand underlying physics of defects. Nonexponentiality in transients are estimated from 
lineshape analysis of the DLTS and TATS spectra. A case study is presented with exper- 
imental evidence of different kinds of nonexponentilaity in different temperature regime 
in case of semi-insulating gallium arsenide (SI-GaAs) . Superiority of TATS over conven- 
tional DLTS kind of analysis in this case is demonstrated with this example. Limitations 
and advantages of various techniques and methodologies used in this work are discussed. 
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2.2 Defect Study by C-V Technique 

Techniques using depletion capacitance rely on the fact that the width of a reverse bi- 
ased space-charge region (scr) of a semiconductor junction device depends on the applied 
voltage. The C-V profiling method has been used with Schottky barrier diodes using 
metal and liquid electrolyte contacts, pn junctions, metal oxide semiconductor (MOS) 
capacitors, and metal oxide semiconductor field effect transistors (MOSFET). 

V e first recall the principles behind profiling semiconductors for the majority carrier 
distribution by the C-V technique. In the case of a Schottky barrier diode, a reverse bias 
produces a space charge region of width W in the semiconductor side of the junction. The 
capacitance (^) is determined by superimposing a small-amplitude ac voltage v on the 
dc voltage V. The ac voltage typically varies at a frequency of 1 MHz with an amplitude 
of 10-100 mV. The capacitance of the Schottky diode is 

C = A^f^^Nscr/iVa - V) (2.1) 

where Nsa- is the ionized impurity concentration in the space charge region. It is also 
assumed that potential drop due to leakage current is small and all the reverse biasing 
potential V drops across the depletion layer. It is assumed that diffusion capacitance is 
negligibly small compared to the depletion capacitance and no transverse non-uniformities 
exist. 

2.2.1 Defect profile in presence of uniformly distributed traps 

It is known that special analyses are required for the interpretation of junction capacitance 
measurements on semiconductors containing deep traps [77, 78, 79, 80, 81]. A deep 

trap is generally defined ^ as a trap whose energy is located below Fermi level in n-type 

^Though a srtict definition does not exist, magnitude of ionization energies and their pressure coeffi- 
cients can be used as distinguishing characteristics between deep and shallow levels. From a theoretical 
point of view, defect levels whose properties cannot be described by hydrogenic effective mass theory are 
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material or above the Fermi level in p-type material. Particularly, deep traps h<i\e been 
found to greatly influence apparent free carrier profiles determined by C-\' meitsuremenTs 


[79, 81, 82]. 


Under the assumption of one sided step junction, Fig. 2.1(a) and 2.i(b) shows the 
resulting band bending of the junction at a given reverse bias \ for the case of deep donors 
and deep acceptors, respectively. The unique feature of a structure with deep traps in 
the region (x-y) in which generation-recombination processes are active in determining 
the equilibrium occupation of the traps. The more general case is exemplified by the rest 
of the depletion region, OjWiy which is fully depleted of mobile carriers, and in which 
trap occupation is controlled only through emission processes. The existence of these 
two regions results in a space-charge nonuniformity. This nonuniformity is represented in 
Fig. 2.1 by a staircase function [83]. This transition region has a ’Debye tail’ of majority 
carrier which spills into the depletion region. This carrier concentration acts to give a 
finite equilibrium population probability of deep traps in that portion of the depleted 
region. In addition, this fuzziness of the depletion width boundary serves to limit the 
spatial resolution of C-V profiles to about a Debye length [84] even in materials without 
deep traps. To find majority carrier profile for uniformly doped substrate, Eqn. 2.1 can 
be written in the form 


Nscrix) 



( 2 . 2 ) 


with X = KgeoA/C. The trap profile can be extracted from the knowledge of background 
profile. 


deep levels. 
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Tran- 



Tran- 





(a) 


(b) 


Figure 2.1: Band diagram of a Schottky diode under reverse bias with one deep trap. Vs is the 
junction barrier height, Vr is the applied reverse bias, Ec is the energy position of the conduction 
band edge, Ed is the shallow donor level, Ep is the Fermi level, . jBt is the trap level, y is the 
point at which Ep and Er intersect and x is the edge of the depletion region. Deep donor traps 
are shown in (a) and deep acceptor traps are in (b). (after ref. [83]) 
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2.2.2 Defect profile in presence of nonuniform distribution of 
traps 

The shape of the profile in the presence of a nonuniform distribution of deep levels is 
now considered. A discussion of the C-V technique applied to nonuniforni distribution 
of trap have been considered by Kimerling [83] and by Schulz [85]. It has been shown 
that great care must be taken in capacitance measurement of free carrier distribution in 
semiconductors containing deep traps. Kimerling [83] has shown that spatial distribution 
derived from C-V can be anomalous unless measured conditions such tis temp^erature and 
frequency of the measurements are carefully determined based on the junction structure 
and the electronic nature of the traps. Schulz [85] has shown that when large concen- 
tration of deep donors/acceptors are present in a certain region of space charge layer, 
one measures an ’apparent’ profile which is shifted from the actual profile of the defect. 
In case of deep donors this profile appears broad and for deep acceptor it appears .sharp 
compared to the actual trap profile. This is due to crossing of the Fermi lev(>l with the 
trap level which causes occupation change deep inside the depletion region where defect 
are present. In such cases, a knowledge of the defect profile helps in determining the 
energy level of the deep traps by locating the cross-over point of Fermi le\'el and trap 
level. 

In case of hea\w damage in the semiconductor due to ion implantation, presence of 
an insulating layer may make a Schottky diode behave more like an Metal-insulator- 
semiconductor (MIS). In case of buried oxide layers, for example, C-V characterization is 
carried out modelhng the structure as M-S-I-S [86]. The doping profile can be obtained 
in such a case from appropriate modifications to standard MIS analysis. With combined 
measurements of capacitance at high and low frequency in regions of inversion, depletion 
and accumulation, a number of parameters such as thikness of the insulating layer and 
interface traps density can be inferred. 
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2.2.3 Factors affecting C-V profiling technique 

C-\' measurements can be affected by many factors which include (i) series resistance, (ii) 
diode edge and stray capacitance, (iii) excess leakage current, (iv) minority carrier and 
interface traps, (v) Debye length, and (vi) voltage breakdown etc. We discuss below some 
of these factors relevant to the present work. 

2. 2. 3.1 Series resistance 

A simple equivalent circuit for a pn junction or Schottky diode consists of a junction 
capacitance C, a junction conductance G, and a series resistance r^,. Conductance deter- 
mines the junction leakage current and can be varied by processing conditions. The series 
resistance depends on the bulk wafer resistivity and on the contact resistance. Capaci- 
tance meter assumes the device to be represented by the parallel circuit in Fig. 2.2, where 
Cm and Gm are measured capacitance and conductance. Comparing the two circuits 
allows Cm to be written as [87] 


O Q 



Figure 2.2: (a) Actual circuit with a series resistance Tj and (b) capacitance meter equivalent 
circuit for a pn junction or Schottky diode. 
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Voltage (V) 


Figure 2.3: Effect of series resistance on the C-V plot in a Schottky device. Simiilatioii param- 
eters : Doping Nd = 1 x 10^® cm~^, Diode area A=7.85xl0“^ cm“, measurement frequency 
f=l MHz. 

= (l+r.G)=f(2,r/r.C7)2 

For the measured capacitance to be approximately equal to the true capacitance, it is 
necessary that -C 1 and {27r frsC)^ <C 1. For junction with low leakage current, 
rsG <C 1 is satisfied. But, (27rfrsC)^ ^ 1 is not always satisfied. For example, if the 
true capacitance is 100 pF measured at f=l MHz then r, should be less than 160 ft 
for 1% accuracy in capacitance measurement. In this way, series resistance interferes 
with capacitance measurement. In Fig. 2.3 the effect of series resistance on the 1/C^ 
vs. V plot is shown for various values of series resistance in a Schottky diode. These 
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2.3 Current- Voltage (I-V) Characteristics 


I-\ studies are very important for complimentary information on carrier transport be- 
haviour in ion-damaged materials. Current transport in a Schottky device under forward 
bias constitutes four basic processes [90] namely, (i) thermionic emission, (ii) tunneling, 
(iii) carrier recombination in depletion region and (iv) minority carrier injection. When 
the diode is forward biased, the current raises exponentially with voltage and in presence 
of series resistance, the part of the applied voltage drops across it for larger current at 
high voltage. The series resistance and ideality factor of the diode can be estimated from 
an iterative fitting of I-V data. However, in presence of an insulating/semi-insulating 
layer between metal and semiconductor contact, the current conduction is determined 
by a number of processes namely, (i) Schottky emission, (ii) Poole-Frenkel emission, (iii) 
Tunnel or field emission, (iv) space charge limited, (v) Ohmic and (vi) ionic conduction. 
In particular, we are concerned with the case of space charge limited current in the case 
of p — i — n type device. Space charge limited current for a unipolar trap free case is 
proportional to the square of applied voltage (i.e. J oc V"^). At low voltage and high tem- 
perature, current is ohmic with exponential dependence on temperature. In presence of 
high concentration of trap, on the log-log plot of current versus voltage, there is a nearly 
vertical section in the characteristics. This is attributed to the filling of a discrete set of 
traps. The voltage at which this vertical transition occurs, is known as trap-filled-limit 
voltage {Vtfl) [91]- 

Trap filled limit voltage measures that fraction of the total concentration of trap that 
is empty in thermal equilibrium. From the knowledge of dimension of the insulator region, 
trap concentration is estimated as [91] 



where k is the relative static dielectric constant, L is in cm. For distributed trap states, 
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the distribution can be extracted from fitting of I-\ characteristics. Huwevei. in MIS 
structure, presence of metal-insulator and insulator-semiconductor junction can hav<> con- 
tribution to the total current making the analysis difficult. .A.part from estimation, 
one can determine trap energy (Et) from the relation [91] 

Ec-Et = kT \n{NJ9gNT) (2.7) 

where g ^ 2, is conduction band effective density of states and & is determined from the 
knowledge of k and g (mobility). In materials with high degree of structural disorder, even 
though chemically very pure, one may expect multiple trap levels distributed in energy. 
A distribution in energ;^’’ of trap in an insulator would reflect in the sharpne.s.s of current 
rise at Vtfl- In principle, energy distribution can be recovered from obscr\’ations of I-V 
curve in such cases. We show later that our samples do show trap limited space-charge 
conduction when forward biased. 


2.4 Thermally Stimulated Capacitance (TSCAP) 


In this method, a reverse biased junction or Schottky diode is cooled to sufficiently low 
temperatures so that thermal emission rates of carriers is negligible. A non-equilibrium 
occupancy of deep levels is then induced with zero bias pulse of duration long enough 
to fill the traps. As the sample temperature is increased, the deep levels relax to their 
equilibrium occupancy. The steady state capacitance is measured as a function of tem- 
perature. Capacitance steps are observed as traps emit their carriers. The midpoint of 
the TSCAP step in temperature (Tm) is related to the activation energy Ee by [92] 


Ec = kTm In 


/?(£. + 2fcr„) 


( 2 . 8 ) 


where P is the heating rate and u is attempt to escape frequency. Trap concentration is 
estimated from the step height of TSCAP curve. Activation energy is determined from an 
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Arrhenius plot of ln(r^//3) versus l/T^ for varying heating rate (/?). One can distinguish 
between majority and minority carrier trap by noting sign of the capacitance change. 

The technique works well for Nt > O.INd and Eg > 0.3 eV i.e. its sensitivity is 
much less than the spectroscopic techniques described earlier. TSCAP measurement is 
less influenced by diode quality. Baseline determination in TSCAP data without trap 
filling is necessary to know the TSCAP step height. Often, there is a large change in 
capacitance during heating cycle even without filling. This is mostly due to Fermi level 
movement with temperature which causes the effective change in depletion width. In 
terms of convenience of use and interpretation, TSCAP is less convenient than DLTS and 
TATS due to requirement of controlled thermal scan rate. 

2.5 Capacitance Transient 

The junction capacitance in a semiconductor diode is related to the trapped charge con- 
centration. If nonequilibrium occupancy of deep traps can be induced, the equilibrium can 
be monitored through measurements of capacitance consisting of capacitance transients 
[93]. Under reverse bias, the traps in the part of the depletion region are unoccupied 
since trap level is above the Fermi level. A non-equilibrium situation is created tj’pically 
by applying a forward bias or zero bias for short duration, during which the depletion 
width reduces. This duration is called filling time, tf. Under filling condition, a trap 
level (say, an electron trap) goes below the Fermi level and the unoccupied traps capture 
conduction band electrons. On reapplication of reverse bias, the carriers emitted from 
the trap subsequently are swept out of the depletion region. As already shown in Fig. 2.1, 
there can be a region within depletion width where the traps are occupied till the point 
trap level crosses the Fermi level. In case of high trap concentration, consideration of this 
region can play an important role. 



Principles of Experiments and Techniques of Data Analy-sis 


40 


2.5.1 Exponential transient 

The measured capacitance transient would be exponential if the rhange in rap.iriianre is 
small as compared to the steady state capacitance and can be expresseti as 


C{t) = Co 




( 2 . 9 ) 


where Co is steady state capacitance value at reverse bias, tit is occupied trap concentra- 
tion and Nd is shallow doping concentration, Tg is characteristic emussion time eon.stant 
of trap. The condition that change in capacitance should be small as I’i >iar’ar- d to Cg 
is same as the condition Nt/Ncl ^ 1 i.e. for small trap concentration. Use of constant 
capacitance voltage transient does not have any limitation on trap den.sity. This has been 
discussed in the next section. The emission time constant in the al)ove equation has 
an Arrhenius dependence on temperature given by 


— = an <vth> Ngexpi-ErlkT) ( 2 . 10 ) 

Te 

where Et is the activation energy of the trap and 0-„ is the capture {Tos.s-section of 
the trap. Other symbols have their usual meaning. Hence. moTiittjriiig of capacitance 
transients at various temperature allows determination of as a function of Kunperature 
and hence the trap parameter Ey and cr„. Experimental characterization is performed by 
use of Arrhenius plot to obtain the activation energy and attempt to escape frequency 
from the slope and intercept, respectively. 


2.5.2 Nonexponential transient 

The capacitance transient signal represented in Eqn. 2.9 is an exponential. However, 
nonexponential capacitance transient may occur when (i) trap density is of the same 
order as that of shallow dopants, (ii) several trap levels decay with similar emission rates, 
and (iii) emission rate is field dependent. There are large number of papers in which 
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attempts have been made to take into account non-exp onentiality due to each of these 
reasons in interpretation of data [94, 95]. There have been also varieties of experimental 
ways of avoiding nonexponentiality [96]. There are large variety of methods available for 
extracting the trap parameters of exponential signals [97]. 

The nonexponentiality due to high trap density can be avoided by performing the 
transient measurement in constant capacitance (CC) mode. In this method capacitance 
or the depletion width is held constant during the carrier emission measurement by dy- 
namically varying the applied voltage through a feedback circuit [96, 98, 99]. The voltage 
transient in that case is expressed as 

Vit) = - Nrexpi-t/re)] + (2.11) 

without any restriction on ^ ratio. Equation 2.11 shows V(t) response to be exponential 
in time. One of the limitation of the CC transient measurement is its slower circuit 
response due to feedback circuit. Initial part of the transient can be affected due to 
carrier capture even during the emission process [100]. Some of the sample dependent 
causes of nonexponential transient may include series resistance effect and leakage current 
effect in the device. 


2. 5. 2.1 Series resistance effect 


Distortion in capacitance transient due to series resistance is also well known [101, 102]. 
Carrier freezeout or trap compensation may lead to temperature dependent series resis- 
tance. For non-leaky devices with series resistance r^, measured capacitance transient 
Cm{t) is related to the sample capacitance C's(t) by 


l-ka;2r2C?(t) 


( 2 . 12 ) 


where u) is the test signal frequency of the capacitance meter. 
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It is clear from Eqn. 2.12 that C^{t) will go through a numnia when > I. 

At this point of inversion (time tm) of the nonmonotonic transient. t,n hs siu li tfiat the 
condition = 1. This is easily obtained by differentiating Eqn. 2.12 witii respect 

to time and setting it to zero. Hence, 

r = I = (2.13) 

* UjCsitm) 20jCm{tm) 

Evaluating for a particular temperature (in case of isothermal transient), we get from 
Eqn. 2.12, 


Cs = 2cJ 

for t < tm 

(2.M) 

1 + \/l 

Cs = 2Cm-^^^ 

for t > tm 

(2.15) 


where x = 2ursCm{t). The presence of series resistance in the transient is first detect (?d by 
noting peak shape in the transient and the actual transient is constructed using Eqn. 2.14 
and 2.15 [103]. 

The applicability of this method relies on the fact that series resistance should be ap- 
preciably high such that > 1 is satisfied. If is not so high, one can add e.xtonial 

resistance in series with the device to fulfill the condition stated above. However, when 
both majority carrier and minority carrier traps occur at the same time, nonmonotonic 
transient may also result. One may mistakenly attribute such occurrence to series resis- 
tance. It is advisable to cross-check the presence of series resistance from current-voltage 
characteristics of the device as well as from impedance and loss measurement. Therefore, 
independent estimate of series resistance and prior knowledge of types of traps (majority 
or minority carrier) are needed. We have already discussed ways of obtaining series re- 
sistance, in section 2.2.4.I. Particularly, we use impedance and loss measurement in our 
device to estimate the effect of series resistance in capacitance measurement. 

Series resistance effect in transient can be reduced by reducing the measurement fre- 
quency u! [104]. However, it severely restricts the range of usable time constants, also low 
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frequency measurement gives rise to noisy data. In presence of defect /damaged region in 
the device, at lower frequencies one normally observes increase in capacitance which in 
turn does not help in reducing the quantity The iise of constant capacitance (CC) 

mode in which Cg is kept fixed during transient overcome the series resistance problem 
altogether [105]. Since dc resistance does not affect the CC mode of operation, voltage 
transient carries trap information alone. 

2. 5. 2. 2 Leakage current 

Device prepared under poor vacuum condition and surface preparation, ion implanted 
devices often exhibit high reverse leakage current. If the capture of the leakage related 
carriers in the space charge layer is larger than the thermal emission rate from a particular 
gap state to the band, then the apparent emission rate is given by [106] 

Sapp = e„ + (2.16) 

where is true emission rate, J^ak is leakage current, cr„ is capture cross-section. Par- 
ticularly, for the case of broad density of states, leakage current has dramatic effect on 
capacitance transient [107]. A correction scheme for leakage current effect on the capac- 
itance transient has been proposed by Dmowski et al. [106]. 

2.6 Deep Level Transient Spectroscopy (DLTS) 

Though several methods are known for extracting time constant from the capacitance 
transient directly, the advantages of spectroscopy in transient analysis was first realized by 
D.V. Lang [71] who introduced the rate window concept in the form of DLTS to deep level 
characterization. DLTS converts the exponential transient signal into a spectrum which 
is amenable to standard spectroscopic interpretation. The signal S(T) is constructed from 
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the capacitance transient C(t,T), at a fixed rate window {tj, t-z) ■ 

SiT) = CituT)-C{t2,T) 


S(T) has a maximum at a temperature Tm- For an exponential 
e„ at temperature Tm is given by 


fin = 


h — h 


and peak value of the signal f(r) is given as 




transient, emission rate 


(2.18) 


(2.19) 


where 



Note that the signal height depends on strength of the transient Co and trap concentration 
nr- From different choice of rate windows, the dependence of emission rate on temjxTat.ure 
can be obtained in a single temperature scan. The ratio r in Eqn. 2.20 is chosen to 
optimize the signal/noise ratio and peak shape resolution. Along with high .sensitivity, 
DLTS provides a quick and simple way to identify all the traps in the depletion layer of 
a diode provided one is able to change trap occupation by external excitation such as 
voltage pulse. From the lineshape of the spectrum, extent of nonexponentiality can be 
judged. 

Nonexponentiality due to high trap density is avoided by implementing the DLTS 
measurement in constant capacitance mode (CC-DLTS) of operation where voltage tran- 
sient instead of capacitance transient is recorded as a function of temperature [108] and 
processed. However, we do not perform CC-DLTS measurement in this work due to large 
change of capacitance with temperature. 
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Temperature (K) 


Figure 2.4: Effect of series resistance on DLTS spectra for different Q values. Simulated peak 
corresponds to a trap with Et= 0.5 eV, cr = 1 x 10“^® cm^, C=40 pF; tl=7ms, t2=14ms. 

DLTS signal will assume complicated form in presence of very large trap concentration in 
the depletion layer [110]. 

2.6.3 Trap parameter extraction 

In case of exponential transient, DLTS signal peak temperatures corresponding to differ- 
ent rate windows will follow an Arrhenius relationship. Slope of ln{en/T‘^) versus l/T 
plot yields activation energy of the trap and intercept gives capture cross-section (a). 
The cross-section determined in this way may be affected by electric field, degeneracy of 
the level and temperature dependent capture rate. More direct estimation of cr is done 
by varying pulse width technique in which fraction of occupancy for a particular level is 
monitored as a function of fill pulse duration t/ [98, 111]. Trap concentration .s usuinated 
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from DLTS signal peak height (using Eqn. 2.22) for a chosen value of t i/ ^ i ratio. Howe\-er. 
this relation holds only for small trap concentration compared lo the l^aekgratmnd doping. 
For larger trap concentration, trap density can be estimated from (. C-DLIS method or 
from transient amplitude measured in constant capacitance mode of operation. ( onren- 
tration profiling of traps in DLTS measurement is done by selecting the region of interest 
through appropriate changes in magnitude of majority carrier pulse for difft'renf .scans in 
temperature. 

Broadening in energy for capture barrier or emission energ}' is determineti from the 
broadening in DLTS peak shape. A broad peak may results from contimiou.siy di.striltuted 
states in the band gap as in the case of amorphous Si. Majority and minority carrier trap 
can be distinguished from the sign of the peak. 

Many variants of the conventional DLTS have been proposed for characterization of 
semiconductors including frequency scanned DLTS, Laplace transform DLTS. Fourier 
transform DLTS, optical and scanning DLTS, interface trapped charge DLTS. lockin DLTS 
etc. [112]. In the next section, we will describe a variant of frequency scanned DLTS 
technique which we have used extensively in our work. 

2.7 Time Analyzed Transient Spectroscopy (TATS) 

This is an isothermal spectroscopic technique where time window is varied keeping the 
temperature constant unlike DLTS where temperature is varied keeping the time window 
fixed. Isothermal spectroscopic techniques have been used in the literature in various 
forms [113j 114, 115]. TATS is one such isothermal spectroscopic technique based on 
difference signal in time window. The TATS signal is given by 


S{t)=C{t,T)^C{t + jt,T) 


(2.23) 
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where C represent isothermal capacitance transient at temperature T and 7 is an ex- 
perimentally chosen constant. For exponential transient of time constant Tg, S(t) has 
a maximum when plotted against ln{t) and maximum occurs at time tm given by the 
relation 


T^m 

Zn(l -t- 7) 

The peak value of the TATS signal is given by 


(2.24) 


- C'o2^^2_^^(1 + 7) ^ (2.25) 

where Co is the prefactor of the transient C(t,T). 

There are several advantages of TATS analysis over DLTS. TATS is a spectroscopy 
in the time domain alone and therefore possible temperature dependence of the transient 
prefactor containing occupancy etc. does not occur. Additionally, in methods involving 
temperature scanning such as DLTS, the lineshape is dependent on trap parameters and 
the temperature regime of observation. 

In contrast, the width of the TATS peak depends only on the parameter 7 , which is 
chosen to optimize resolvability and signal to noise ratio. The full width half-maximum 
(FWHM) of the TATS lineshape is a constant given by [116] 


FWHM = ln{U) - ln{t-) = 2.496 


(2.26) 


for 7 = l,t+ and being the times corresponding to the half maximum of the TATS peak. 
This can be obtained in the same manner as Eqn. 2.21 was derived for DLTS lineshape. 
Because of this, it is very easy to detect and estimate the degree of nonexponentiality from 
TATS lineshape. A convenient way of characterizing nonexponentiality is to represent it 
by a stretched exponential transient of the form 


C ~ Coexp 



/ 3<1 


(2.27) 
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where (3 is the stretching factor. This form of transient, also kiifwn lUS illiain- Watts de- 
cay [117], has been invoked in many physical systems including defect anai>-sis ; 1 1 S. 1 igi. 
Constructing the TATS signal for transient of Eqn. 2.27, it can be shown Ity differentiation 
that signal goes through a peak at tm given b}^ 

/^n /3 _ /?In(l + 7 ) 

V/ “(H-7)^-l 

which reduces to the expression in Eqn. 2.24 for /3 = 1. Simple calculation shows that 
the peak position in time tm does not change appreciably on lowering J. Specific, ally, for 
7 = 1, the peak position (^) is 0.990 times the perfect exponential case even for a 3 as 
low as 0.5. Hence, estimation of the time constant is independent of the stretching factor 
P in most cases. The experimental inaccuracies are much larger than expc'cted due to 
nonexponentiality factor p. 

Nonexponentiality in transient normally shows up as broadening of the lineshapc*. If we 
use Eqn. 2.28 for peak position then the FWHM of the TATS signal is proportional to 
( 1 //?) for the case of 7 = 1 . Hence, a quantitative measure of nonexponentiality is obtained 
simply by inspection of TATS spectral lineshape. For more precise determination of the 
P factor, one can use precomputed values of FWHM as a function of p. For example, for 
the case of 7 = 1 , we have simulated the FWHM of the TATS lineshape for various values 
of P, which fits into the following straight line formula : 


1 _ FWHM - 0.078 
P ~ 2.417 


(2.29) 


for determination of p correct upto three decimal places. In many cases, it is more 
meaningful to give a physical interpretation to the broadened peak in DLTS or TATS 
spectra. A distribution in deep energy levels can give rise to nonexponential transient of 
the form [ 120 , 121 ] 


C'(i) = f “ 9iEai)Coex.p{[-en(EaPt]dE^ 

*3 0 


(2.30) 
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where the distribution function often assumed to be Gaussian around the mean energy 
value £’ao is 

Here, 25\/2 ln2 is the FWHM of the broadening of the energy level. This situation arise 
in case of alloy semiconductors and defect located in strain field as in the case of less 
perfect crystal. From TATS lineshape fitting, exact value of Ea^ and g{E — Cj) can be 
estimated. In case of DLTS lineshape fitting, other temperature dependent factors may 
complicate the analysis. 

For better resolution of emission time constants, TATS can be performed using higher 
order spectroscopies [122, 123, 124, 125]. For any gamma value, TATS signal of different 
orders are given as 

1st order : Si(t) = C{t) - C{{1 + 'y)t) (2.32) 

2nd order: 52(t) = C'(t) — 1.5C'((1 + 7 )t) + 0.5C'((1 + 7 )^t) (2.33) 

Zrd order : Sz{t) = C{t) — + 7)t) + g^((l + T)^^) “ (2.34) 

A hierarchy of higher order spectroscopies can be used, where suitability is dictated by 
signal to noise ratio. In case of large signal, higher order TATS is very useful in resolving 
close peaks as shown in Fig. 2.5. Here two peaks with T 2 /T 1 —A which are not so well 

resolved in first order spectra is resolved clearly in 3rd order spectra. We have taken 

advantages of such higher order analysis for our studies in this work. The formalism can 
be applied to the case of other spectroscopies namely DLTS, PATS [126] and TWTS 
[125]. Although we have mentioned many advantages of TATS over DLTS, it should be 
pointed out that for TATS analysis, data is required over very wide time scales for a 
given temperature. For example, peaks occurring in the temperature range 80K-250K in 
DLTS will have time constant variation in the range of 10® sec - 10"® sec for a particular 
temperature within the above mentioned should be noted 
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log(time/s) 

Figure 2.5: Comparison of different orders of TATS to resolve two closely peaks : (1) 

first order TATS, (2) second order TATS, and (3) third order TATS. 

that while DLTS is used for surveying defects, TATS can be used for detailed analysis 
of defect phenomena related to a particular defect or defects lying close in energy. In 
particular, TATS is more versatile in analyzing different kinds of nonexponentiality of 
transients encountered in this work. 

2.7.1 Analysis of nonexponentiality : A case study 

We present here a case study in analyzing nonexponentialit}' of photo-induced current 

transients in imdoped semi-insulating GaAs (SI-GaAs). The superiority of TATS over 

conventionally used photoinduced transient spectroscopy (PITS) [127], which is based on 

DLTS like analysis, is demonstrated in this section using our own experimental data 

^This case study involves current transients instead of capacitance transients. Only commonality of 
techniques with the rest of this work is emphasized. 
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In SI-GaAs, at the end of an intrinsic light pulse, normally one expects to observe 
decaying photocurrent transient leading to positive peak in PITS spectra. However, most 
often negative peak arising out of rising current transient have also been observed in 
these materials. In spite of inherent nonexponential nature of the photoinduced current 
transients, several workers have attempted to fit multiple exponentials [128, 129] for the 
observed transient. 

For acquiring experimental data, we have used commercial undoped SI-GaAs sample 
cut from 400-fim thick wafer. We used the planar structure with electrical contacts in the 
form of strips on front surface by evaporating Au/Ge alloy in vacuum. The metal contacts 
are annealed at 450°C in flowing forming gas. The samples were mounted on a liquid 
nitrogen cooled cold-finger type cryostat with temperature controlled to within ±0.1K in 
the temperature range 90K-350K. A copper-constantan thermocouple was mounted on 
the sample holder close to the de\dce. Light pulses were provided from the He-Ne laser 
source (A=6328 A) with the help of a mechanical chopper with cut off times less than 0.2 
ms. The current transients are digitized at each temperature with the help of Keithley 
Source Measure Unit (Model 236). The complete setup, shown in Fig. 2.6, is PC controlled 
except for temperature control. Each acquired transient is then used for further analysis 
using TATS. The time constants are calculated using Eqn. 2.24 and stretching factor or 
nonexponentiality is estimated using Eqn. 2.29. 

2. 7.1.1 Typical current transients and PITS 

Typical photocurrent decay curves are shown in Fig. 2.7 for two different temperatures. In 
SI-GaAs both these types of transients are commonly observed. In our discussion, we will 
refer to transients labelled as (a) in Fig. 2.7 as decaying transients and to curves of type 
(b) as rising transients. Deca3dng photocurrent transient in these materials is interpreted 
to be due to slow thermal release of carriers from traps filled during light pulse. Most of 
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Beam splitter 



Figure 2.6; Schematic of photocurrent transient measurement .setup. 


the excess carriers would have recombined at a much faster rate normally due to band 
to band recombination. Hence the slow decaying transient carries information regarding 
the emission rate of carriers from the defect center. However a ’clean' interpretation of 
rising transients such as curve (b) has not been possible and many models have been 

advanced mainly through photoinduced current transient spectroscopy (PITS) studies 
[ 127 , 128 , 130 ]. 

As far as transients of type (b) are concerned, note that in most earlier studies only 
the rising part of the transients have been analyzed. However our data clearly shows that 
all such transients also have a slow decaying component reaching a minimum and then 
rising to a saturation. We have taken care to confirm that the decaying component is 
not related to any artefact of experiment such as time constant of chopper or R-C circuit 
response. This is further borne out by the fact that the slow decaying component is highly 
temperature dependent and is always observed along with negative transient. Hence we 
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PITS analysis questionable. Hence, PITS is not well suited for detailed studies involving 
questions of phenomena, though it can be used as a general purpose survej^ technique. 
Specifically, it fails miserably in trying to understand negative peaks or rising transients 
in SI-GaAs. 

2. 7.1. 2 TATS of positive peak 

As has already been pointed out in section (2.7), there are several advantages associated 
with TATS analysis, it being a completely isothermal technique. To demonstrate its 
efficacy we choose at first to analyze the decaying transients which lead to positive peak 
in PITS spectra in the temperature range 165-200K. In Fig. 2.9(a), we show T.A.TS signals 
for these transients at several temperatures. The emission time constant is obtained 
from the peak position using Eqn. 2.24. The corresponding Arrnenius plot is shown in 
Fig. 2.9(b) which yields an activation energy of 0.41 eV with corresponding capture cross- 
section of 6.1 X cm^. This can be attributed to the level called EL6 reported in the 
literature [130]. 

As discussed in section (2.5), inspection of the TATS spectra directly show that tran- 
sients are nonexponential. This is clearly shown in Fig. 2.10 where experimental TATS 
points are shown by symbols and a perfect exponential of time cdhstant corresponding to 
the maximum (tm) is shown as a dashed curve. From experimental FWHM, the stretch- 
ing factor p is obtained using Eqn. 2.29 and the corresponding TATS spectra is shown 
by bold line in the same figure. This clearly shows the power of TATS in determining 
degree of nonexponentiality of transients in a straightforward way. Similarly, stretching 
factor determined at different temperatures range between 0;67-f].90, the lower value be- 
ing obtained at lower temperatures. The significance of this systematic deviation from 
exponentiality with temperature is not clear to us but it could be related to recapture 
of carriers. The height of TATS peak at different temperatures do not scale with (l/tm) 





Figure 2.9: (a) A series of TATS spectra in a SI~GaAs sample at different temperatures : (A) 
174. 3K, (B) 179. 7K, (C) 184.5K, and (D) 188. 8K. (b) Arrhenius plot for deep trap corresponding 
to positive TATS peak shown in (a). 
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Figure 2.10: Comparison of experimental spectra with calculated ones to demonstrate the 
degree of nonexponentiality. The experimental spectra (with symbols) are broader than the 
spectrum corresponding to perfect exponential (dotted) line. The bold line corresponds to a 
stretched exponential with /3=0.668 determined from the FWHM of the experimental curve. 

as current transient prefactor includes e„, emission rate of trap. This is mainly due to 
the change in the degree of occupancy achieved for a fixed exposure time at different 
temperatures. Normally this would contribute to lineshape distortion in PITS spectra. 

2.7.1.3 TATS of negative peak 

We use TATS to analyze the nature of nonexponentiality of rising transients (Fig. 2.7) 
observed at higher temperature. As in PITS spectra, rising transients lead to negative 
peaks in TATS spectra for several temperatures. Note that the sample shows two negative 
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peaks shown in Fig. 2.11(a) as opposed to what is usually observed in PITS spectra. This 
is possibly due to the fact that the smaller negative peak gets subtiieiged in aheaiiy highly 
distorted PITS hneshape. 

The dominant negative peak in TATS spectra is due to the rising |)art of tite current 
transient as shown in curve (b) of Fig. 2.7. Note that the lAFS .signal fails to zero sharply 
on the faster side of the peak. In fact it is initially positive corresponding ti) the decay- 
ing part of the transient. The TATS signal goes through zero at a time approximately 
corresponding to the minimum of the transient. Hence, it is clear that hy no stretch of 
imagination the transient corresponding to the negative peak can !>e con.sidered i*xponen- 
tial or mildly non-exponential by choosing to focus only on the later rising part of .such 
transients. 

However, for thermally stimulated relaxation processes, use of Arrhenius plots to esti- 
mate the approximate energy involved is a robust procedure not withstanding ntjnoxpo- 
nentiality. Hence, we treat the TATS peak shifts with temperature as indicator of shift 
in time constant and hence plot an Arrhenius plot given in Fig. 2.11(b). The activation 
energy so obtained is 0.79 eV with an estimate of capture cross section from the inter- 
cept as 2.3 X 10~“ cm^. The activation energy of negative peaks found in earlier works 
[130, 131] normally range between 0.65 to 1.2 eV. Therefore, it has always been attributed 
to midgap centers, though a definite identification has remained controversial. It has been 
suggested earlier that presence of large dark current plays a crucial role in giving rise to 
negative peak [128]. In contrast, note that we have obseix'ed negative peaks even at 
temperature well below room temperature where the dark currents are negligibly small. 

Figure 2.12 shows TATS of transients at a particular temperature (290K) for various 
values of filling time i.e. the time for which the laser light is on. Change in height of 
TATS peak of the main negative peak is an indicator of slow increase in the corresponding 
filled trap concentration. However, it is seen that the second negative peak height remains 
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Figure 2 . 11 : (a)TATS spectra with 7=1 in SI-GaAs samples at different temperatures : (A) 
293.9 K, (B) 297.7 K, (C) 302.5K, and (D) 307.4 K showing two negative peaks, (b) Arrhenius 
plot for the dominant negative peak of (a). 
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Figure 2.12: TATS spectra for SI-GaAs at temperature 290 K for different filling times : (A) 
20 ms, (B) 70 ms , (C) 100 ms, and (D) 200 ms. 

constant indicating that saturation of the peak occurs in time scale faster than the shortest 
filling time scale used in these experiments. The extent of trap filling is found to be 
highly temperature dependent and studying filling behavior is rewarding by an isothermal 
spectroscopy such as TATS. It also demonstrates that the height of conventional PITS 
peaks cannot be considered as healthy indicator of extent of filling. 

To explain such kind of nonexponentiality revealed from TATS, we have proposed a 
simple kinetic model [105] without attempting to pinpoint any physical mechanism. The 
most significant feature of our model is to introduce a new rate process which acts as a 
delay between capture of a carrier and its emission. For simplicity, we will assume that 
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onh one dominant electron trap and one band, say, conduction band is involved in the 
process. An electron gets captured at an empty trap taking the defect to a metastable 
state from which it must relax to the ground state. It can then finally emit the carriers 
to the band. We can write down the rate equations as follows : 

dA^ 

— —A^A^C + ( 2 . 35 ) 

dA* , , 

— = A+A+a - rA* (2.36) 

with the constraint that A° + A* + A+ = 1, where A+, A*, A° are the concentration of 
empty, occupied metastable state and occupied ground state respectively, r is the rate of 
conversion from A* to A° and C is the product of capture cross section (cr), average thermal 
velocity {vth) and the total trap concentration. It is also assumed that the concentration 
of electrons in the band is equal to the concentration of empty traps, an assumption 
mostly valid in these materials in the dark. The rate equations are written with all 
concentrations normalized with respect to the total trap concentration. The desired delay 
between capture and emission is introduced by the negative term in the Eqn. 2.36 where 
it is treated like an emission term in so far as it depends only on the concentration of the 
metastable state. The rate equations can conveniently be summarized in a state transition 
diagram as is given in inset of Fig. 2.13 [116]. 

These rate equations are solved numerically on a personal computer. The choice of pa- 
rameter fin is taken from the experimental Arrhenius plot, normalized capture coefficients 
and initial concentrations are chosen to correspond approximately to EL2 parameters in 
SI-GaAs. The rate r is chosen so as to produce qualitative features of our data. Fig- 
ure 2. 13 shows the TATS of a particular solution obtained numerically and compares it 
with a typical experimental TATS curve for negative peak. The ability of any model 
to mimic cusp like behavior of the transient or positive to negative crossover in TATS 
spectrum puts severe restriction on the class of models that need be considered. The abil- 
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Figure 2.13: TATS spectra for the experimental data (symbols) and calculated data (bold line) 
from the proposed model with EL2 parameters. The chosen parameters as defined in Eqs. 2.36 
and Eqs. 2.37 are : en=19.5 s”\ 0=1x10^ s“\ rc=l, A+(0)=0.011, A*(0)=0-96. 
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Wafer -> 

n'^n-Si 

p^p-Si 

Orientations : 

[lOOl 

[11 ij 

Epilayer thickness : 

22 /i m 

22 pm 

Doping concentration : 

1.2x10^® P /cc 1.5x10''’ B /a: 

Back contact (Ohmic) : 

A1 

Al 

Schottky contact : 

Au 

Al, Pd 

Implanted species : 

Ar'*', Au'^ 

Ar'*', Au'^ 

Respective energy (MeV) : 

1.45 , 4.6 

1.45, 4.6 

Ar'^ ion doses (cm'"^) : 

5-10x10'^ 

1-5x10'-* 

Au^ ion doses (cm”^) : 

5 X 10® 

5x10® 

Annealing (Post Imp.) : 

160-600°C 

160-600'’C 


Table 3.1: Details of sample preparation conditions. 


ion irradiation. 3MV Tandem Van de Graaff accelerator facility at Institute of Physics, 
Bhubaneswar, India, was used for Au"^ ion implantation. 

3.2.1 Schottky diode fabrication 

For the present study, epitaxially grown silicon wafers of both n-type and p-type were 
used. Details of the device preparation conditions are listed in Table 3.1. The wafers are 
cleaned in ultra pure trichloroethylene (TCE), acetone and methanol sequentially in an 
ultrasonic cleaner for few minutes at room temperature. Then the wafer is etched in buffer 
hydrofloric acid (HF) solution for few seconds to remove surface oxide and then cleaned 
in deionized water thoroughly and dried. The cleaned wafers are put immediately in a 
vacuum chamber for deposition of A1 (thickness ~ 1500A) on the back side of the wafer 
for making ohmic contact. The metal is deposited using a thermal evaporation system. 
In case of non-epitaxial wafers i.e. CZ wafers, back contact is annealed at 485°C for 30 
minutes in vacuum for ohmic contact formation. 

The wafers with ohmic contact on back side are then put into an electron beam 
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evaporation system for deposition of gold on n-type Si, and a.1 nmi ni iim /p a.]1 a.HiiiTn on 
p-type Si for Schottky contact formation. This evaporation system is equipped with a 
liquid nitrogen trap to have a clean vacuum in the deposition chamber. The evaporation 
is done under 10~® Torr pressure and about 2000A thickness of metal is deposited on 
the front side (mirror finished) of the wafer. The thickness is monitored using a quartz 
crystal controlled digital thickness monitor during deposition. High purity (99.99 %) gold 
wires and aluminium from ingot are used for deposition of Schottky contact. Metal masks 
each having many circular holes of 1 mm diameter were used during metal deposition for 
defining Schottky diode area. The wafers are scribed to obtain pieces of convenient size. 

Ion irradiation (described in the next section) are performed on finished devices or on 
bare wafprs from front side of the wafer. For n-type Si samples, irradiation is done prior to 
any contact formation and for p-type Si, irradiation is done on finished, devices and on bare 
wafers which are cut into pieces. In case of Schottky diodes prepared after ion irradiation, 
the wafer is carefully degreased and subsequently cleaned in buffer HF solution for oxide 
removal. The wafers are cut into small pieces having one/two Schottky diode(s) on each 
piece and mounted on standard T05 header for measurements. The contact from metal 
side of each Schottky diode is taken with fine gold wire (25 diameter) bonded with 
polymer based conductive silver paste (A and B compound, Elteck Qorp., India). Other 
terminal is pasted with the body of the header. The epoxy contact is cured in oven at 
70°C for 30 minutes for proper strength of the bonding. A cap is put onto the header 
for protection from mechanical damage. All the devices underwent this low temperature 
heat treatment for curing prior to any measurement and these devices will be referred to 
as as-implanted. 
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3.2.2 MeV ion implantation 

High. 6n6rgy Ar"^ and Au”^ ions are used for implantation on siiit on. High t’lH'ru} Ar ions 
are produced using Van de Graaff accelerator. It is a widely used and ronmu'.nl>' awiilable 
electrostatic accelerator. A schematic diagram of the accelerator sy.stem is .shown in 
Fig. 3.1 This accelerator (xModel AiV-2000) is a precision high intensity .source of 

positive ions. The accelerator produces a beam of ions which is inten,s<'. homogeneous, 
stable and controllable over a wide range of energies. It is provided with thre^e gas bottles 
and presently they are filled up with He, Ar and Kr gases. An r.f. ion source i.s used to 
produce positive ions. 

The mass and energy^ analysis of the accelerated ion beam is done by an analyzing 
magnet. The analyzing magnet system consists of a 16 KGauss electro-magnet cotipled 
with a d.c. regulated (0-30 V) power supply. During operation the magnet is cooled to 
280K by circulating chilled water through the coils. The stainless steel switching chamber 
has five ports at angles 0*’, ±20° and ±45°. The energy^ and the mass analyzed beam can 
be directed to the experimental chamber through any one of the five ports. 

The accelerated ions obtained from the accelerator reach the expcrinn'iit:'.! chamber 
through a clean and highly evacuated beam line. The ion beam after leaving the analyzing 
magnet passes through beam stopper followed by a beam profile monitor. The beam profile 
monitor is used to monitor the position, shape and uniformity of the beam. The beam 
then passes through a specially designed beam defining X-Y slit arrangement mounted on 
a double cross chamber which is used for selecting a uniform portion of beam of required 
size onto the target placed in the irradiation chamber. 

The heavy ion (Ar) irradiation chamber (shown in Fig. 3.2) is mounted on a turbo 
molecular pump to obtain a hydrocarbon free clean vacuum. In the center of the chamber 
a vertically movable sample holder is placed on which many samples can be mounted at a 
time. The holder is electrically isolated from the rest of the chamber. A current integrator 



Control Room 



Figure 3.1: Schematic diagram of the 2 MeV Van de GraafF accelerator system. 
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Figure 3.2; Schematic diagram of the heavy ion irradiation chamber. 
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(ORTEC Model 439) which measures the total charge accumulated during irradiation is 
connected to the sample holder. A cylinder made up of aluminium sheet surrounding the 
sample holder acts as secondary electron emission ( from sample) suppressor. 

For Au"*" ions irradiation, 3 MV tandem Van de Graaff accelerator facility at Institute 
of Physics, Bhubaneswar, India was used. The irradiation was performed on cleaned 
n-type and p-type Si wafers at room temperature at an energy of 4.6 MeV with beam 
current of 0.2 nA for 60 seconds. This resulted in a dose of 5x10® cm~^ in both n-Si and 
p-Si. With the help of a beam scanner, an area of 10mm x 5mm was exposed to Au^ ions. 
The irradiation was performed 2-4° oif surface normal on Si for avoiding channeling. 

1.45 MeV Ar'^ ions are irradiated at a dose rate 5 x 10^^ cm“^ s~^ for various duration 
resulting in total doses of 5 x 10^^ cm~^ & 1 x 10^^ cm~^ for n-type Si and 1 x 10^'^ cm~^ 
& 5 X lO'-^ cm“^ for p-type Si. The current was chosen such that negligible beam heating 
occurs during irradiation. The irradiation is performed at normal incidence to the surface 
of the wafer which have 4° misalignment cut to avoid channeling. An area of 3.5x3. 5 
mm^ uniform beam was falling on the sample at normal incidence. 

3.2.3 Thermal annealing 

For high temperature annealing of the irradiated samples, the irradiated wafers after ul- 
trasonic cleaning with TCE, acetone and methanol are kept in a quartz tube which is 
evacuated by diffusion pump equipped with a liquid nitrogen trap. The furnace used for 
annealing has a ceramic tube, wound with nichrome wire and insulated with asbestos pow- 
der. A uniform temperature with control of ±2K could be obtained in the central portion 
of the furnace. The furnace is kept in a horizontal movable platform (see Fig. 3.3). After 
annealing, the cooling of the sample was done under room temperature. As-implanted 
wafers are annealed in vacuum for | hour at temperatures 400°C' and 600°C. For low 
temperature annealing, the finished devices (ion irradiated) are kept in oven set at 160°C7 
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Figure 3.3; Schematic of the furnace used for vacuum annealing at iiigh teinperanin^s. 

for 2 hour, 1 hour and 2 hours in sequence. Few samples were stored at room temperature 
(20°C'-35°C') for more than one year and some of these samples were used for studying 
possible annealing due to room temperature storage. 

3.3 Experimental Setup and Measurements 

work involves capacitance measurements as a function of temperature, time and 

tage. The same setup is used for transient as well as steady state measurements. 

pt for temperature programming, the setup is made fully computer controlled for a 

ty of measurements. A block diagram of capacitance measurement setup is sho^vn in 

Fig. 3.4 In the following subsections various parts of measurement setups and procedure 
are described. 
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Figure 3.4; Block of the capacitance/voltage transient and DLTS/TATS measurement 

system. 
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the capacitance range is chosen as 1 pF or 3 pF, respectively or its multiple of 10. For 
a typical C-V measurement, at first temperature is monitored and voltmeter is triggered 
(internal) to measure the applied voltage and capacitance simultaneously in two channels 
of the high speed voltmeter. The control is done through IEEE-488 interface bus. The 
set up is fully computer controlled except temperature control. The data is transferred 
from internal buffers of the high speed voltmeter to the computer for further processing. 

3.3.3 Capacitance transient measurements 

The sample is kept under reverse bias, and a desired temperature is achieved through 
proper control of heater current in the cryostat and liquid nitrogen level. Bias voltage 
pulses and trigger pulses (for synchronous measurements) are generated from the PCL 
718 card along with the help of a multiplexer (Fig. 3.7). Since the analog outputs from 
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plug-in boards in the computer are very noisy, a simple trick is used to oht aiii clean pulses 
with programmable vmltage levels and duration. The noisy software progiiunnutiile analog 
outputs from the card are heavily filtered and fed to the variotis input.s ui a multiplexer 
(e.g 4051). Bias voltage pulse is generated by s%vitching analog channids. Tlie selection 
of multiplexer channel is controlled by setting digital outputs of the same Ixiard (DO. Dl. 
D2). Filling pulse widths down to few micro seconds could be easily achirwed. The trigger 
pulse for synchronization are also obtained from the digital outputs t)f the plug-in board. 
The capacitance meter (Boonton 72B) had a fast response time option of 50 /rs installed 
in it. For capacitance transients, the analog output of the rapacitam-e meter is digitized 
using high speed voltmeter Keithley 194A with 16-bit accuracy. 

For each transient, a maximum of 30,000 data points are stored at sampling rates 
ranging from 1 ^s to 1 s. An acquired transient with 30,000 data points at ecpial intervals 
covers about 4| orders of magnitude in time. For longer durations, two transit'nt.s are 
acquired with different sampling rates. The data is generally e()iiipre,'--.s(‘fi and stored in 
binary format. A typical transient measurement consists of the following .st'quence of 
operations : 

1. Device is kept under reverse bias to ensure emptying of traps above Fermi level (Ef) 
in the depletion region. 

2. Temperature is stabilized in the cryostat to the desired degree of accuracy. 

3. Filling pulse is applied leading to partial or complete occupancy of traps in the 
depletion region. 

4. The high speed voltmeter (Keithley 194A) is triggered for data acquisition at the 
end of the fill pulse. 

5. The acquired data is transferred from Keithley 194 A to the PC for further processing 
at the end of the emission transient. 



Figure 3.8: Schematic of the feedback circuit for constant capacitance control 

3.3.4 Voltage transient (constant capacitance) measurement 

In general, by applying a voltage pulse across the device, one measures capacitance tran- 
sient which carries trap information. However, to avoid nonexponential transient due to 
large trap concentration, capacitance is held constant during carrier emission measure- 
ments by dynamically varying the applied voltage during the transient response through 
a feedback circuit. Figure 3.8 shows the schematic for a constant capacitance feedback 
circuit. The analog output proportional to the capacitance (in the range 0-3 V) was 
used as feedback. The reference was generated from PCL 718 AD card. Control was 
obtained using a proportional integral (PI) controller. To improve stability, all redun- 
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dant/unintentional poles in the feedback loop were eliminated. Zero biasing of the sample 
during filling is done by collapsing the feedback using a JFET lor the d(!siied duration. 
The relevant pulse is applied at the gate of the JFET for the desiied duration. This 
scheme was essential for higher temperature transients as the C-V curve vtis nearly flat 
in the allowed voltage range. 

The sample capacitance is a nonlinear function of temperature and bias voltage. The 
response time of the controller w'as kept usually within 1 ms to ensure stability over a 
large temperature range and for low noise. Fast response times up to 50 /is could be 
obtained by careful tuning of the parameters of the PI controller irrespective of choice 
of bias, filling and temperature. Filling time experiments are carried out using such fine 
tuning at a particular temperature. 

The measurement steps are exactly identical to the steps of acquiring capacitance 
transient as described earlier. 

3.3.5 Deep level transient spectroscopy (DLTS) & time ana- 
lyzed transient spectroscopy (TATS) measurements 

For deep level transient spectroscopy (DLTS) measurements the same setup as shown 
in Fig. 3.4 is used. The capacitance waveform (C-t) is digitized at each temperature 
during a temperature scan. Only one temperature sweep of the sample is required since 
the entire C-t waveform is obtained at different temperatures. Typically, a filling pulse 
of duration 50 ms was used in most of the DLTS measurements. The device is cooled 
under reverse bias and during heating cycle, at each temperature typically 10 transients 
are taken and averaged. This improves the signal to noise ratio dramatically. For this 
reason, temperature sweep is done at slow rate. Typically, DLTS signal corresponding to 
7 different rate windows are stored in a single temperature scan for obtaining Arrhenius 
plot. On line plot of DLTS signal is done on the computer screen for a particular rate 
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window. One can store the entire waveform to check the exponentiality of the transient. 
Various signal processing functions can be performed on C-t data to improve signal/noise 
ratio, analyze nonexponentiality, resolve closely spaced peaks. Special care was taken to 
eliminate noise and electromagnetic interference as described in the next section. 

Due to large change in capacitance during temperature scanning, DLTS is not per- 
formed in constant capacitance mode of operation. This would require high reverse bias 
to be applied through the feedback circuit to keep the capacitance constant. 

For time analyzed transient spectroscopy (TATS) analysis, at a stabilized temperature 
the entire transient (C-t) is acquired with 30,000 data points and data in logarithmic time 
is extracted after proper filtering. Typically, only 600 data points are stored in logarithmic 
time scale. Signal is constructed in the way described in the previous Chapter. For 
obtaining an Arrhenius plot, transients at several temperatures are taken. Transients can 
be acquired for various filling times at one fixed temperature for capture cross-section 
measurement. 

When the capacitance is held constant with the help of feedback circuit and isothermal 
voltage transient is taken to construct spectroscopic signal, that measurement is termed 
as constant capacitance (CC) TATS. In this work most lineshape analysis of spectroscopic 
data is carried out with CC-TATS measurements as the nonexponentiality due to large 
trap concentration is absent in such transients. 

3. 3^6 Thermally stimulated capacitance (TSCAP) measurement 

In thermally stimulated capacitance (TSCAP) measurement, the device is cooled under 
reverse bias and at lowest temperature the traps are filled with majority carriers at zero 
bias or forward bias. Then the device is reverse biased, and heated at a constant rate 
to monitor the steady state capacitance as a function of temperature. Capacitance steps 
are observed as traps emit their earners. Similar scans without the filling pulse at low 
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temperature are referred to as unfilled XSCAP curves and used as referent <* for isulating 
capacitance changes due to trap emission. 

3.3.7 Current vs. voltage (I-V) measurement 

For I-V measurements, programmable voltage is provided by PCL 718 card along with 
extra circuitry for obtaining both positive and negative voltage. Current is measured by a 
GPIB controlled Keithley picoammeter (Model 495). The sample is mounted in the cryo- 
stat (as described earlier) to monitor and control temperature. At a stable f c'lnpi'rature, 
current and voltage data is taken for logarithmic /linear variation in voltage. The setup 
is computer controlled except for temperature control. 

3.3.8 Impedance measurement 

The implanted devices having large density of traps are tested for low frequency capaci- 
tance, impedance, loss, series resistance measurements. The measurement is done using 
a GPIB controlled Hewlett Packard impedance analyzer (Model 4192A) which can mea- 
sure in the frequency range 5Hz-13MHz with test signal level lOmV-lV. At first, zero 
adjustment is done in open and close circuit configuration separately. Impedence and loss 
measurement at different temperatures was performed using an ac signal level of 10 mV 
at a frequency of 1 MHz. This set up is also computer controlled. 

3.3.9 Noise elimination zmd filtration 

The sensitivity of the transient measurement setup largely depends upon the degree of 
noise control and filtration. Despite large signal due to high trap concentration, stringent 
noise control was essential for higher order spectroscopy and lineshape analysis. 

Two major sources of noise were electromagnetic interference leading to 50 Hz pickup, 
and high frequency noise from the computer. By appropriate shielding, grounding and 
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filtering, noise in the voltage bias/capacitance reference was reduced to less than 1 mV 
peak to peak. Heavy filtration dampens the response time and prevents fast switching of 
biiis levels necessary for filling pulse. To overcome this problem, two heavily filtered bias 
levels were independently provided from PCL card and a multiplexer was used to switch 
between tw'o levels as already described in section 3.3.2. The transient signal is processed 
using moving window averaging and further multi-point averaging whenever necessary. 
After filtration, the data set is reduced to 600 evenly spaced points in the logarithmic of 
time. This data is then subjected to the required spectroscopic operations to obtain the 
spectra. 
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Voltage (V) 

Figure 4.2: Comparison of C-V characteristics of unimplanted, Au'^ and Ar^ implanted n-Si. 

lower voltages. Since depletion width (W) is proportional to 1/C, the figure clearly shows 
that zero bias (or, even for any particular applied bias) depletion width systematically 
increases on increasing the dose of implantation. 

Similar qualitative features, though to a lesser extent, are observed for Au'^ implanted 
silicon even though the doses used are only 10^ cm~^. Figure 4.2 shows comparison of C-V 
characteristics at room temperature for unimplanted, Ar+ implanted and Au'*' implanted 
Si. As the range and dose of Au"'" ions are much smaller compared to Ar"'' ions in n-Si, 
onfe expects less effects of implantation on C-V characteristics for Au"'" ions. 

Figure 4.3 shows d(l/C'^)/dV vs. depletion width W, as is normally done to obtain 
conventional carrier concentration profiles, for unimpalnted and implanted samples ob- 
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Depth (yum) 


Figure 4.3: Apparent carrier concentration profile for unimplanted, Au"*" and Ar"^ implanted 
Si for different doses. 

tained by differentiating the curves of Fig. 4.2. For unimplanted samples, carrier profile is 
seen to be uniform and in accordance with the background doping of the epitaxial layers. 
For implanted samples, the apparent carrier profiles show a sharp rise for the lowest W 
corresponding to flat region in the respective C-V characteristics. The depths at which 
these sharp rise occur are very large i.e. at about twice the zero bias width of the control 
diodes in case of Au'*’ implanted Si and about 4 times in case of Ar"’’ implanted Si. The 
sharp rise in concentration profile is due to presence of flat region in C-V characteristics. 
This feature in carrier profile is an artefact due to large trap density in a region of deple- 
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Figure 4.4: Vacancy distribution for 1.45 MeV Ar"*" (5 x 10^^ cm and 4.6 MeV Au"*" (5 x 10® 
cm“^) ions as predicted from TRIM simulation. 

tion layer. The location of this trap dominated region lies far beyond the range predicted 
by TRIM [15] (Transport of Ions in Matter) which is only 1.23 fj,m. Large concentration 
of traps are expected to be present in as-implanted Si, but its location is unusual. 

The TRIM predicted vacancy concentration profile for Ar'*' and Au+ ion implantation 
are shown in Fig. 4.4. The peak concentration of the defect profile shows that for Au+ 
ion implantation (dose 5 x 10® cm~^) defect peak concentration is about 1500 times lesser 
than the Ar+ ion (5 x 10^^ cm~'^) case. Note that implanted dose of Au+ ion is 4 orders 
of magnitude less compared to Ar"^ ion case. Being heavy in mass Au"^ ion produces 
significant amount of defect even for such a low dose. The damage peak positions are 
slightly shifted due to diflierence in the implant energy. However, this predicted position 
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of the peak defect concentration is quite shallower than the mciisureci \<ilue <ts shown in 
Fig. 4.3. This is an important observation in the present work where defects are seen to 
be present at a place much larger than the expected range from TRIM simulation. 

The flatness in C-V curve has been reported earlier in Mg'^ and Ni^ implanted Si 
[132] and in He ion implanted Si [133, 134, 135] without a discussion of its origin. Singh et 
al. [136] have attributed occurrence of such flat regions to interface states in Xi/nCdF 2 
Schottky diode using ideas suggested by Fonash [137]. High series resistance in the 
device can give rise to such feature [88] as it was shown with simulation in Chapter 
2. Consequently, this yields wrong carrier concentration profile. However, independent 
measurement of dissipation factor (a;RsC where is series resistance and C is actual 
capacitance) at various temperature using an impedance analyzer showed no significant 
effect of series resistance in the capacitance measurement of irradiated device. The dis- 
sipation factor was found to be varying in the range of 0 to 0.17, larger value for higher 
temperature, in the temperature range of 100K-310K. The series resistance effect become 
significant when dissipation factor is comparable to or greater than 1. Moreover, jis has 
been shown in the last chapter, large series resistance give rise to a minima in the 1/C‘^ vs. 
V curve rather than flatness as observed. Hence, this effect is a significant characteristics 
of the damaged region in the depletion width. We go on to present more parameters 
dependence on C-V curves in order to isolate the origin of these features. 

4.1.2 Temperature Dependence 

C-V measurements were carried out at various temperatures in the range 100K-310K. 
Figure 4.5 shows 1/C^ vs. V plot for low dose implanted Ar"^ ions in Si measured at 
various temperatures in the range 100K-310K. In these measurements, the sample is 
initially kept at steady state reverse bias at a constant temperature, and then C-V data 
is taken by changing the bias from high to low reverse bias. The change in capacitance 
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Figure 4.5: C-V characteristics of as-implanted n-Si measured at temperatures (Tl) 303. 9K, 
(T2) 278.3K, (T3) 257.4K, (T4) 237.9K, (T5) 217.2, (T6) 186.6K, (T7) 155.5K and (T8) 103.5K. 

with temperature is seen to be large as shown in Fig. 4.5. 

Note that as temperature is decreased from Ti to Tg in the figure, the depletion layer 
edge for any particular bias (say, zero bias) lies deeper in the sample. Further, the extent 
in voltage of flatter region gets reduced and the curves at the lowest temperatures are 
fully linear indicating constancy of carrier concentration. These distinct features taken 
together points to the involvement of traps on controlling the nature of the curves. Large 
changes in capacitances, say at zero bias, is an indicator of changes in the charged trap 
occupancy within the depletion layer with temperature. Clearly trapped charge is such 
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that the depletion width is large at lower temperature. For higher temperat«re.s and 
lower voltages, the traps seem to affect the C-V curve, their orcupancy l)eing controlled 
by crossing of the Fermi level with trap level within the depletion layer. Tlu'se sy.stematics 
were observed in ail our damaged samples. The increase in depletion width in presence of 
damaged region seems to be a consequence of acceptor nature of the dominant traps in the 
damaged region. These interpretation further gets supported by hysteresis mtuistirements 
described in the next section. 

4.1.3 Hysteresis Effect in C-V 

The trap dominated nature of the damaged layer is confirmed by observation of hystere- 
sis behaviour in the C-V measurement. Figure 4.6 shows IfC^ vs. plots at a certain 
low temperature for different voltage sweep rates in both directions of voltage change. 
A typical scan begins with decreasing reverse voltage from a steady state value. Note 
that the C-V characteristics depends on both sign and rate of voltage sweep clearly show- 
ing hysteresis effect. The occurrence of hysteresis can be easily understood in terms of 
presence of large number of trapping centers in concentration comparable or larger than 
the background doping. During C-V measurement, if the voltage sweep rate is made 
faster/comparable to the trap emission rate, one observes a difference in capacitance be- 
tween decreasing and increasing voltage sweep. For the fastest sweep rate the difference 
between the two curves are maximum and for slow sweep rate upper curve approaches the 
lower curve as shown in Fig. 4.6. It can be noted that the C-V characteristics are identi- 
cal for different sweep rates when the voltage is decreased owing to fast capture process 
at the traps, while it is different on increasing the voltage since it is now controlled by 
rate of emission. Thus it is clear that flat region or nonlinearity in C-V is due to trap 
Ailing during decreasing reverse bias. From transient measurement data presented later, 
we And that the sweep rates used here are comparable to the emission rate of a single 
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Figure 4.6: 1/C^ vs. V plot showing hysteresis behaviour for decreasing and increasing reverse 
bias with different voltage sweep rates. Arrows indicate the direction of voltage sweep. 

major trap identified in this work. The zero bias depletion width of most devices having 
damaged layer are 2-3 /xm larger than their corresponding undamaged device. Also note 
that the depletion width is larger when the reverse bias is increased after allowing capture 
during C-V measurements. This implies that after capture of electrons, the defects get 
negatively charged in so large number that the depletion region is widened to uncover 
required amount of positive background charge to maintain charge balance condition. 
Hence these traps play the role of dominant compensating center. 

Similar hysteresis effect has been observed for Au+ implanted Si as well. In Fig. 4.7, 
we choose to display such a behaviour with one difference. In this case the degree of 
forward biasing controls the degree of filling of traps. The sweep rate is kept constant for 
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Figure 4.7: C-V characteristics with decreasing and increasing bias showing hysteresis effect in 
An'*' implanted n-Si. Slanted arrow indicates the direction of voltage sweep, and different set of 
curves were obtained for different degree of forward bias injection. The forward bias voltage for 
each is indicated by vertical arrows above the voltage axis. 

different cycles. During the decreasing voltage cycle we take the device to fonvard bias 
indicated as negative voltage in the figure. With increasing forward bias, we see larger 
hysteresis effect i.e. difference between curves corresponding to decreasing and increasing 
voltage cycle increases. The arrows in the figure indicate the voltage upto which the 
forward bias is applied before starting the increasing voltage cycle. It can be noted 
from the largest forward bias curve in Fig. 4.7 that initially the difference in capacitance 
between decreasing and increasing voltage cycle is negligible. This is partly due to sweep 
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Figure 4.9: Forward I-V characteristics of 160°C annealed sample at 302K showing ohmic and 
space charge limited (SCL) current region. Voltage corresponding to nearly vertical section of 
the characteristics is trap filled limit voltage (Vtfz.)- 

replaced by a metal and i-region is nearly intrinsic due to the above defect level. Transport 
behaviour on metal/implanted (amorphous) Si contacts have been earlier explained using 
a model of charge injection into a thin layer of trap rich amorphous Si [138]. Similar 
model has been proposed for proton bombarded GaAs [139]. 

4.2.2 Temperature Dependence 

Temperature dependence of I-V characteristics for 160°C annealed {Ar'^ implanted) n-Si 
are shown in Fig. 4.10. The two distinct regions of current conduction are seen clearly 
for three different temperatures. The notable feature in these curves is decrease of Vtfl 
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Figure 4.10: Forward I-V characteristics of Ar'^ irradiated (160"C annealed) Si measured at 
different temperatures showing shifts in Vtfl- 

with decreasing temperature. The observed shift of V^fl can be explained on the basis 
of trap filling with lowering temperature at a fixed bias. As the Fermi level moves closer 
to the conduction band in n-Si upon reduction in temperature, more traps get filled and 
thus number of unfilled fraction reduces with lowering temperature. Thus trap filled limit 
occurs for lower voltages at lower temperatures. 

We will demonstrate in the next chapter similar effect in TSCAP spectra of Ar"*" ion 
implanted n-Si where trap filling is observed even without applying filling pulse. Thus 
unambiguous determination of space charge limited current in irradiated device proves 
the presence of fully trap controlled region in the depletion region of the irradiated diode. 
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Space charge limited current measurements have been used in the literature to deter- 
mine energetic distribution and spatial profile of bulk traps in semiconductors [91, 140]. 
However, nonuniform spatial distribution of trap, voltage and temperature dependence of 
occupancy does not allow accurate determination of energetic distribution of traps from 
I-V characteristics. Nevertheless, for purposes of rough estimation, under the assumption 
of uniform trap distribution over the full depletion layer, we can use Eqn. 2.7 of Chapter 
2 in this case. For Vxfi,= 1-7V (Fig. 4.10), estimated trap density is Nt = 5.5 x 10^“* 
cm~^ for layer width of 2 /im. It is to be noted that the Vtfl measures that fraction of 
total trap that are empty in thermal equilibrium. Moreover, effective length over which 
trapping and detrapping is occurring during bias change is quite small compared to the 
total depletion width. Thus the trap concentration peak value is much larger than the 
value estimated above. For our purposes, it suffices to note that very large concentration 
of traps occur within the depletion layer controlling both its I-V and C-V characteris- 
tics. In the next section we present results on similar measurements for high temperature 
furnace-annealed samples. 

4.3 Effect of Annealing on C-V and I-V Character- 
istics 

We have annealed few irradiated samples at higher temperature prior to Schottky diode 
fabrication. The irradiated wafer pieces are annealed in vacuum (« 2 x lO'® Torr) for 
30 minutes at 400‘’C and 600"C respectively. The effect of annealing on C-V and I- 
V characteristics of the devices further help in relating trap activity to the degree of 
damage. 

4.3.1 C-V Characteristics 

Figure 4.11 shows l/C^ vs. V plot for Ar+ ion irradiated n-Si before and after annealing at 
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Figure 4.11: Comparison of C-V characteristics at 300K of unimplanted, as-implanted and 
annealed n-Si implanted with 5 x 10'-^ cm“^ Ar"^ ions. 

various temperatures. We have justified previously the cause of large depletion width in 
as-implanted device. It is clearly seen that high temperature annealing decreases the trap 
concentration. This leads to systematic changes in both the major features in 1 /C‘^ vs. V 
curve. W ith increase in annealing temperature (i) the zero bias depletion width recovers 
towards the control samples, and (ii) the flat region for low bias reduces progressively. 
The temperature range of annealing (400"-600°C) itself suggests that the dominant defects 
controlling electrical behaviour are point defects, since extended defects begin annealing 
out only at temperatures above lOOD^C. Comparison with unimplanted device shows that 
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Figure 4.12; Hysteresis in C-V measurement at 298.2K during decreasing and increasing voltage 
sweep for 400°C annealed sample. 

600°C is not sufficient to remove all the defects completely as the depletion width for 
the implanted samples is larger than the unimplanted case. The near parallel nature 
of the curves in Fig. 4.11 for high reverse bias region indicates absence of defects in the 
deeper region. For 400®C annealed samples, we could observe hysteresis effect in C-V 
measurements at room temperature as shown if Fig 4.12. It was obtained for fast sweep 
of voltage during decreasing and increasing bias voltage. This is also indicative of changed 
time constant of dominant defects controlling hysteresis. This indeed is further confirmed 
from transient measurements presented in the next chapter. 
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Figure 4.13: Forward and reverse I-V characteristics of Ar"^ implanted samples after 400°C 
and 600°C annealing. 

4.3.2 I-V Characteristics 

We show in Fig. 4.13, comparison of I-V characteristics at 300K for 400°C and GOO^C 
annealed samples. The reverse saturation current is similar for these two samples. But the 
forward characteristics of 600° C annealed samples is considerably improved as compared 
to the case of 400°C annealed samples. For the later case space charge limited current 
above trap filled limit voltage is clearly seen. Though Vtfl is not so large as it was for 
as-implanted sample, presence of still large trap concentration is clear. It is expected that 
these defects can be annealed out by annealing above 600°C. These defects are anticipated 
to be due to point defect complexes rather than extended defect which forms after high 
temperature annealing in some cases. 
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4.4 Discussion 

4.4.1 Summary of Major Experimental Features 

It would be appropriate here to reiterate the salient features in the results presented in 
the preceding sections. 

• High dose heavy ion irradiation causes damage which manifests in C-V and I-V 
characteristics of Schottky diodes made from these samples. Specifically, a very 
striking feature is that the capacitance does not change with bias for a voltage 
range in the low bias regime. 

• The damaged layer is obsen’-ed to be within the depletion layer with its width being 
much larger than control samples. The width is also larger than ion range or damage 
location predicted from TRIM simulations. 

• Hysteresis and temperature dependence of C-V characteristics, are clearly indicative 
of involvement of traps in damaged layer embedded within the depletion layer. The 
observation of space charge limited conduction in forward bias also indicated the 
presence of a dominant trap in large concentration. 

• The increase in irradiation dose increases the zero bias width. Annealing at tem- 
peratures in the range 400-600°C reverses the effect to different degrees indicating 
recovery and damage relaxation. 

4.4.2 Qualitative Explanation of Major Features 

Let us first consider flatness in low-bias regime of C-V curves or bias independence of 
capacitance and consequent nonlinearity of 1/C^ vs. V curves. It is well known that a 
trap can be occupied in the transition region of the depletion layer where the trap level 
is below the Fermi level. Figure 4.14 shows a band diagram for a Schottky diode with 
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Figure 4.14: Energy band diagram for a Schottky diode with one acceptor trap level at energy 
Er communicating with conduction band. Yr is applied reverse bias and Yd is built-in voltage. 

one uniformly distributed deep level at energy Er which gets negatively charged when 
occupied with electron. Changes in occupancy occur within depletion layer at the location 
where Er and Er cross each other due to band bending. In our case we are dealing with 
large concentration and nonuniform spatial distribution of traps. It is worth mentioning 
here that in standard cases of small trap concentration, band bending is mainly controlled 
by background shallow dopants. For this case, in contrast, trap concentration is so large 
that it controls band bending in the region of its presence and which in turn determines 
the crossing between Er and Ep. Therefore, the location of the crossing is determined by 
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4 . 4.3 IVIodsI S imu.l3.tjioiis of C“V Chopcictsristics 

The purpose of model simulations is to arrive at the essential requirements on charge 
profile that would reproduce the major features of C-V characteristics. We seek to obtain 
reasonable estimates of concentration of charges and the distances at which they are 
located without attempting to provide numerical fit to the experimental data. 

Several simplifying assumptions are required in order to keep the results of simulation 
transparent. To account for the presence of dominant trap indicated in the experiments, 
we assume that there exists only one midgap deep acceptor trap due to damage. The fact 
that its energ}' is in the middle of the bandgap would be shown to be fully justified on the 
basis of spectroscopic measurements in the next chapter. The region of depletion layer 
from the surface to the location of the damage induced trap is assumed to be uniformly 
compensated having low charge density. In reality, the degree of compensation would 
be vart'ing over distance with a profile more or less similar to that predicted for damage 
distribution. Compensation would mainly be due to deactivation of shallow dopants and 
presence of deep defects in the region. However, the edge of the depletion layer is unable 
to penetrate into this layer during our capacitance measurements preventing changes in 
charge distribution in that region. Therefore an assumed effective compensation would 
suffice to take into account the voltage drop across this region. We also assume that 
the metal-semiconductor barrier height is same as that for control samples (i.e. standard 
value of 0.8eV for Au on n-Si Schottky). This value only appears as a constant in the 
magnitude of voltage. 

One dimensional Poisson equation is solved numerically for a particular choice of trap 
profile to obtain the band bending in the semiconductor. The voltage dependence of 
capacitance is obtained by integrating over different widths of depletion layer. The trap 
occupation, band bending, and carrier concentrations are determined self-consistently 
during numerical integration. For the purposes of trap occupation, Fermi occupancy 
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Figure 4.15: Charge profile due to trap concentration Nx and background doping showing 
different parameters considered for simulations. 

function is replaced with usual step approximation i.e. all traps lying below the Fermi 
level are considered to be occupied giving rise to a step in the occupancy at the Er- 
Ep crossing point. Non-degenerate carrier statistics was only needed for all simulations. 
Depletion approximation was not explicitly assumed, though most of the results reported 
are well within the ambit of this approximation. 

The parameters of simulation are specified in terms of a charge profile as depicted in 
Fig. 4.15. The relevant parameters are 

1. the mean position (xo) of trap profile (assumed to be a Gaussian), 

2. the FWHM {2\/2ln2 a) of the Gaussian distribution 

3. the maximum concentration of trap Nx, 
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4. the width of the compensated layer (xc) and the level of compensation specified in 
terms of fraction (f) of background doping (N^). 


The Poisson equation can now be written as 

_ Pt + Pd 
dx^ e 


(4.1) 


where 


Pt =~qNTexp 


Pd = 


{x - Xoy 
2 ct 2 

f QNd/f if X < Xc 
I qNd X > Xc 

Typical results of a simulation are shown in Fig. 4.16 where energy band diagram, 
trap concentration profile and the resulting 1/C^ vs. V are depicted in sequence. The 
C-V characteristics shown reproduced the major features in the experimental cur^-es. The 
trap concentration and location are adjusted to obtain capacitance and voltage close to 
experiments. Recall that the flatness in the C-V curve is due to the fact that changes in 
trap occupation account for extra voltage drop without necessitating changes in depletion 
width. 

In order to evaluate how the choice of parameters aifect the calculated C-\* curves, 
we perform simulations varying each parameter at a time holding others constant. Fig- 
ure 4. 17(a)-(d) shows the sets of C-V characteristics which clearly display different depen- 
dences. They can be summarized as follows: 


• The sharpness of the transition from flat to linear region is controlled by FWHM of 
the trap distribution (Fig. 4.17(a)). 

• The zero bias width is relatively insensitive to changes in trap concentration (Fig. 4.17(b)) 
but very sensitive to the choice of its mean position (Fig. 4.17(a)). The mean posi- 
tion largely controls the level of capacitance of the flat region without affecting the 
linear portion. 
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Figure 4.16: Calculated energy band diagram, trap occupancy profile and corresponding simu- 
lated C-V characteristics for a Schottky diode with one acceptor trap level at energy Ex'. Two 
sets of energy band diagrams are shown corresponding to cases when : (1) all traps in the 
distribution are occupied, (2) only the shaded portion of the charge profile is occupied. 
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Figure 4.17: Different parameter dependence of simulated C-V characteristics for a Schottky 
diode with one acceptor trap level at energy : (a) FWHM of Gaussian distribution, (b) trap 
concentration (Nr), (c) mean of trap profile (x^), and (d) compensation (N^/f). 
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ature is decreased the bulk Fermi level of the semiconductor moves up in relation to the 
metal Fermi level causing the crossing point to move towards the surface. This leads to 
increased occupation of the traps; hence increased negative charge forcing the depletion 
layer in turn to move out. At low temperature the traps filled in the process stay occupied 
leading to both larger zero bias depletion width, and linear 1/C^ vs. V curves. 

Therefore, the information regarding electrically active defects would be from the 
Er/Eir crossing point. In principle, by manipulating this crossing point in space one 
can obtain charge profile information. However, the problem of identifying the absolute 
location in space where this crossing occurs is not straightforward. The measurement 
of capacitance only gives the location of the edge of depletion layer and the crossing 
point relative to this edge would depend on the details of the profile inside. We have 
shown that the change of slope from the flat region to the linear region in 1/C^ vs. V 
curves is rather sharp as represented in apparent concentration profile (Fig. 4.3). The 
trap concentration profile at the crossing point is very sharp and can be treated as an 
’’effective electrical interface” between trap-free and trap-dominated regions of the sample. 
The nature of this crossing point is distinctly different from that of a hetero junction or 
insulator-semiconductor interface capable of inversion. We have in course of our work 
ruled out these possibilities on the basis of C-V studies over extended region of bias. The 
doses used were below amorphization threshold and it was several orders of magnitude 
lower in case of Au'*' ion implantation. Yet the nature of characteristics are similar. 

The deactivation of shallow dopants due to damage and charges due to electrically 
active traps produce a compensated region with low charge density till the end-of-range 
of the implanted ions. The degree of compensation would be expected to have a profile 
more or less similar to damage distribution as obtained from TRIM simulations. 



Charge (cm* ) 


Nd=1.2x 10^^ cm'^ 


1 



1.e+15 

0 

-1.e+15 
-2.e+15 
-3.e+1 5 
-4.e+15 
-5.e+15 
-6.e+15 
-7.e+15 



Figure 4.18: 

(h\ TV, 
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Figure 4.19: Hysteresis effect shown with simulated C-V curves assuming a fixed rate of emission 
from one acceptor trap level at energy Er- 

of emission of carrier from the trap into the simulations. A typical set of calculated cur^'es 
demonstrating hysteresis in close agreement with experiment is shown in Fig. 4.19. 

4.4.4 Physical Processes : Defect Migration k Clustering 

Armed with insights from simulation, we are now in a position to discuss the most per- 
plexing aspect of the results in this chapter. In the last section we have already seen 
that except for the last fraction of a micron, the depletion width is mostly due to the 
distance of the sharp edge of the trap profile. Typical depletion width in samples contain- 
ing damage is between 2.5/im-4/im. In contrast, the prediction of TRIM simulations give 
an end-of-ion range of 1.23/im only. Hence the question arises as to why the trap profile 
extends so much beyond the location of the end-of-ion range damage. There is also the 
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the mechanism of recovery is unclear at this stage, it is most probably due to spatial 
variation of degree of clustering. It is known that binding energy for interstitial clusters 
increases with cluster size given by Euin) = 2.1 - 1.45(V^ - where n is the 

number of interstitials in the cluster [52]. Hence higher temperature is needed to repair 
regions containing larger clusters which are closer to the damage. Though more detailed 
studies are required to pin down the mechanisms, our work demonstrates that capacitance 
based studies can serve as an important tool in such an endeavour. 

4.5 Summary & Conclusions 

In this Chapter we have presented results on steady state capacitance-voltage (C-V) 
measurements on Schottky diodes of n-Si containing damaged buried layers. It has been 
shown that characteristic features such as occurrence of bias independent capacitance 
region, hysteresis, strong temperature dependence of capacitance and space charge limited 
conduction are controlled by large concentration of a dominant acceptor trap within the 
depletion layer. The effective electrical interface between trap free and trap dominated 
region is defined by crossing between trap level and Fermi level. Model simulations of 
C-V characteristics were carried out to explain major features and to obtain essential 
requirements of the charge profile. A sharp trap profile is seen much beyond the ion range 
and its location is found to be dependent on irradiation dose and annealing temperature. 
The detection of sharp profile of traps at large distances has been explained on the basis 
of interstitial migration and clustering. 


Chapter 5 

Deep Level Spectroscopy of Defects 
in Deep Buried Layer 

The essential physics of ail forms of thermally stimulated space-charge spectroscopy is 
based on the thermal emission transient at some fixed temperature. In Chapter 2 we 
ha\e discussed the formalism and analysis of transient signal and derivation of spectro- 
scopic signal. One of the principal advantages of spectroscopic analysis, such as DLTS, 
is the presence of a flat baseline i.e. there is high common-mode rejection of the steady 
state capacitance variations in spectroscopic signals. Therefore, it is possible to obtain 
extremely high sensitivity e.g. AC/C ~ 10’® for C~100 pF. Moreover, lineshape analy- 
sis in spectroscopic signal is extremely convenient for interpretation of related physics of 
defects in comparison to transient analysis. However, the task of extending use of these 
spectroscopic methods to the emerging class of defect dominated materials poses consid- 
erable challenge and MeV ion induced damage in buried layers is a case in point since 
validity of most standard analysis may be questionable. 

In the previous Chapter, we have seen manifestation of various unusual features in 
steady-state capacitance and current measurements due to large density of defects in the 
damaged layer. In this Chapter, we will see more pronounced effect of these defects on the 
transient measurements. We carefully monitor capacitance transients and use a variety 
of spectroscopic techniques, identifying in the process sources of difficulties in standard 
interpretation. We have used TSCAP, DLTS as survey techniques. Having realized the 
problems of conventional DLTS based on capacitance transient measurements, we have 
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utilized the constant capacitance (CC) voltage transient measnrenn’nts. and TATS has 
been used for spectroscopic analysis. Large change in capacitance with temperattire does 
not allow implementation of DLTS technique in constant capacitance modes. Moreover, 
other advantages of isothermal spectroscopic technique as mentioned in Chapter 2 com- 
pared to thermal scan technique dictates the extensive use of TATS in CC mode for this 
work. 


5.1 Capacitance Transients : Some Unusual Fea- 
tures 

The basic step in all our spectroscopic measurements involves acquisition of capacitance 
or voltage transients. Hence it is worthwhile to concentrate on the features of transients 
at the outset. The transients are monitored at the end of a filling pulse across the test 
diode during -which applied reverse bias is taken to zero to fill the traps with carriers. 
During the filling voltage pulse the time response of the capacitance is termed irs capture 
transient, and its time dependence at the end of the filling pulse is termed as emission 
transient. In general capture of carriers at traps is faster than carrier emission rate at 
a particular temperature, though capture process can be thermally activated to make it 
slow enough at lower temperatures. 

Figure 5.1 shows a set of typical capacitance response on application of a voltage pulse 
(from reverse bias -4V to OV, for a filling time of t/) to Schottky diodes which have 
damage embedded in the depletion layer. Following the application of zero bias, initially 
the capacitance increases sharply due to collapse of the depletion region but almost im- 
mediately it decreases to low capacitance value comparable to reverse bias capacitance 
thus showing a ’spike-like’ behaviour in capacitance response during the filling pulse. Flat 
region towards the end of the filling pulse indicates no further change in depletion width 
for longer filling duration. This has been checked for long durations in many cases to 
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Figure 5.1: (a) Voltage flUiug pulse of duration (t/) 20 ms on irradiated devices, (b) Cap^itance 
transient response at room temperature for 5 x 10“ cm-" Ar+ as-implanted, (c) 5 x 10 cm 
Ar+ and lOO'-C annealed and (d) 5 x itf cm'" Au+ as-implanted samples. 
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detect any possible presence of slow components in capture transient. On reapplication of 
reverse bias voltage (Vr), a very large magnitude of capacitnnre change is observed due 
to emission of trapped carriers. Note that the transient magnitude is large in spite of the 
fact that the net change in capacitance at the end of filling pulse is very small. 

In normal cases only an increase in capacitance corresponding to zero bi;rs is expected, 
sometimes followed by a small transient due to filling of traps within the depletion layer 
from carriers in the Debye tail at the depletion edge. In our case the increase in capacitance 
is seen to be much larger than steady state zero bias capacitance indicating that the 
depletion edge does penetrate for a short while the beyond location of the apparent barrier 
observed through C-V measurements. The return of the depletion width is then clearly 
due to filling of traps leading to large decrease in positive charges within that period. 
Hence ’spike-like’ response is direct consequence of encounter of the depletion edge with 
sharp trap profile with huge concentrations. As has been argued in the hist Chapter, the 
dominant traps must be assuming negatively charged on capturing electrons in order to 
explain such a large depletion width. The occurrence of large magnitude majority carrier 
emission transient further proves beyond doubt that it is trapping within depletion layer 
which controls the observed response. 

The magnitude of the capacitance spike is observed to be temperature dependent. 
.A.t lower temperature, the spike height reduces and capacitance response increasingly 
resembles the normal case. This is also consistent with linearity of the 1/C^ vs. V plot 
at lower temperature as discussed in the last Chapter. The change in Fermi level with 
temperature increases steady state trap occupation and thus reduces the spike height as 
well as transient magnitude at lower temperature. 

Form a closer look at the transient, it is apparent that the steady state of capacitance 
at reverse bias is achieved rather prematurely and suddenly before the characteristic time 
the transient ought to have taken to decay. In this sense, the transient behaviour is nonex- 
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ponential, and e\en taking into account nonexponentiality due to large trap concentration, 
such transient behaviour cannot be explained. A proper study of nonexponentiality re- 
quires use of spectroscopic techniques where lineshape analysis can reveal detail nature of 
the defect phenomena responsible for such unusual features. Since emission transient is 
of large magnitude, standard implementation of the DLTS analysis should be dealt with 
caution. In the next section we will present results from DLTS and thermally stimulated 
capacitance (TSCAP) measurements on these samples. 

5.2 Broad Survey Using DLTS 

The methodology' and setup used for DLTS measurements are given in earlier chapters. 
The control diode was tested for any trap signature using DLTS measurements. No 
traps could be detected upto detection limit of 10“^ of background doping. This ensures 
that any traps detected in irradiated sample would be irradiation induced. The samples 
used for this study were irradiated prior to Schottky diode preparation i.e. irradiation 
was performed on bare silicon wafer pieces with back contact made by thick aluminum 
deposition. Devices were made on these wafers after proper etching/cleaning of surface 
oxide or any other contaminants. 

A set of typical DLTS spectrum for samples irradiated with two different doses of 
Ar'*’ ions and low dose Au"*" ions are shown in Fig. 5.2. For the chosen rate window, three 
majority carrier related peaks labelled A, V 2 and D1 are observed in the temperature range 
90K-310K, irrespective of ion species and dose. The occurrence of a minority carrier peak 
HI has been sample dependent. Though we do not study its systematics, appearance of 
minority carrier peak has been discussed later. Note that the doses used are considerably 
high so as to produce large concentration of defects. The peak A occurs in very small 
concentrations and it is due to well known A-center (oxygen- vacancy pair) [20]. This 
defect has an activation energy of Ec-0.17 eV with a electron capture cross-section of 
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~ 2 X 10 cm . This center is detected in most irradiation works with bulk silicon and 
has been extensi\ely studied [83]. Since it is associated with oxygen, it is less likely in 
epitaxial wafers having low oxygen content. This justifies its very small peak height in 
our DLTS spectra. The peak labelled V 2 is a well known peak due to divacancy center 
(V 2 observed in all ion irradiation processes in n-Si [21, 24, 28, 142]. The peak labeled 
D1 is a newly observed feature and it is the most dominant among all other peaks shown 
in Fig. 5.2. 

The concentration of V 2 related level as calculated from peak height is ~ 5 x 10^^ 
cm This concentration {Nx/Nd ~ 0.05) is small compared to that expected due to 
high irradiation dose used in this work. In MeV Ge implanted Si with dose 1.2 x 10^ 
much larger concentration (Nx/Njj ~ 0.16) of divacancy centers are reported by Svensson 
et al. [18]. The observed small concentration of divacancy centers even for such a high 
dose may be due to several reasons. Firstly, it may be due to incomplete filling of the 
traps using zero bias filling. As the zero bias depletion width is very high as concluded 
from C- V measurements, all the trap levels are not filled during application of filling pulse. 
This results in a peak height smaller than expected. Secondly, due to use of high dose 
rate {~ 10^^ cm"^ s”’-) in our irradiation, divacancy production rate is smaller. Lower 
generation rate of point defects with increasing ion dose rate was reported by Svensson 
et al. [18]. Thirdly, evolution of defect complexes and clusters can be responsible for 
reducing divacancy concentration [19]. Comparatively, larger density of divacancy centers 
are observed in Au"^ implanted n-Si as can be seen in Fig. 5.2. The dose used for Au'*' 
implantation is 4 orders of magnitude lower than the dose used for Ar"*" implantation. 
However, due to heavier mass of the Au"*" ions, divacancy production rate is quite high. 
This is obviously smaller than the concentration expected from TRIM simulation as a 
large fraction of vacancies and interstitials recombine during the irradiation. In most of 

o / 

the low dose irradiation work, doubly negative charge state of the divacancy level (V 2 ) 



Deep Level Spectroscopy of Defects in Deep Buried Layer 


127 



Figure 5.4: DLTS spectra of Ar"^ implauted sample annealed at 160^'’C showing large D1 peak. 
This peak shape is skewed at the right side of the peak. 

centration in all our samples. This has not been observed so far mostly because, to our 
knowledge no other DLTS study uses unanneaied samples with as high doses as in our 
case. An unusual feature in the DLTS lineshape of D1 is the sharp rise in the curve at 
the right side of the peak corresponding to premature end as previously pointed out in 
transients shown in Fig. 5.1. In the rate window implementation of DLTS, if the two times 
ti and t 2 are chosen such that ti fall in the normal region and ta fails in the prematurely 
flat region of the transient, such a narrow peak (or sharp rise) would be expected in DLTS 
spectrum. This particular spectrum shown in Fig. 5.4 corresponds to the transient shown 
in Fig. 5.1(b). In such a case, the peak position shown in Fig. 5.4 may not be the true 
peak position and hence use of Arrhenius plot for estimation of trap parameters will not 
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3 deld correct value. However, if one ignores the effect of such nonexponentiality on the 
Airlieniiis plot, rough estimation of trap activation energy corresponds to a midgap level. 

Note that nonexponentiality due to high trap concentration may affect the DLTS spec- 
trum. Moreo\er, due to large change in steady state capacitance with temperature, the 
DLTS peak heights are strongly dependent on choice of rate window. In a temperature 
scanning method the spectral lineshape would be distorted in such cases. The temperature 
dependent peak height also seem to indicate different extent of filling at different tem- 
peratures. The superiority of a isothermal spectroscopic method has been demonstrated 
in Chapter 2 with a case study. We will see in the subsequent section the usefulness of 
such technique to isolate possible causes of nonexponentiality reinforcing its usefulness to 
defect phenomena. However, it is to be noted that the use of DLTS has enabled us to 
make a broad survey of defects present in such materials with ease, and it further dictates 
the necessity of using different spectroscopies for a particular defect of interest. 

5.3 Detection of Defect Using TSCAP 

TSCAP measurement is one of the simplest methods of detecting high density of traps 
in semiconductors. We have carried out TSCAP measurements for detection of dominant 
trap and for estimation of their concentrations. The measurement is done in constant 
voltage mode as large change in capacitance does not allow the measurement to be made in 
constant capacitance mode. Though this method has limitations as discussed in Chapter 
2, we performed these measurements to get a better estimate of the density of traps. In 
these experiments reverse bias is applied only after cooling to low temperature and is hold 
at zero bias considerable amount of time for the filled case. Hence, filling of traps to a 
greater degree as compared to DLTS is ensured. 

A typical TSCAP spectra for as-implanted (Ar"*" ions) Si is shown in Fig. 5.5 for zero 
bias filling and no filling case. The heating rate (^) was chosen to be 2K/minute and 
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Figure 5.5: TSCAP signal for as-implanted Si with and without zero bias filling. Applied 
reverse bias Vr.=4.5V. Arrows indicate capacitance step positions corresponding to V 2 and D1 
levels. 

both the curves are taken in heating cycle. The filled curve shows one extra capacitance 
step at lower temperature (165K) apart from the increased step height for the other one 
at 210K indicated by arrows. Thus, TSCAP spectra is dominated by two defect levels. 
The low temperature step is due to divacancy related level ¥3 (Ec-0.42) eV which we have 
discussed earlier. Note that zero bias filling is essential for this trap to be seen in TSCAP 
spectra. It can also be noted that the change in capacitance corresponding to this step 
height is considerably large as compared to that found in DLTS measurement. The large 
step at 210K is due to damage related trap (Dl) identified in this work. This step height 
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is a better indicator of the large magnitude of concentrations present in the as-implanted 
samples. Fioni the step height in TSCAP, one can estimate the trap density using Nd-dlL 
where dll’ is change in depletion width due to capacitance step. For Fig. 5.5 this amount 
is 1 X 10" cm-'A Similar values were also estimated from C-V measurement at room 
temperature. W'ith known heating rates and reasonable values of capture cross-section as 
found from transient measurement, we estimated the trap energy which confirmed that it 
is the same major trap as identified in DLTS measurement. 

The large step height in unfilled curve of Fig. 5.5 indicates that traps are filled just by 
lowering the temperature and with zero bias filling these traps are additionally filled in 
larger amount. This can be explained on the basis of shift in Fermi level with temperature 
lowering which in turn causes occupation of large amount of additional trap at the Ef/Et 
crossing. Due to large occupation of the traps, the band becomes almost fiat over the 
distance such defects are present. This trap filling without application of filling pulse 
indicates that trap occupancy is controlled by quasi Fermi level which in turn is controlled 
by trap. 

Figure 5.6 shows TSCAP spectra without zero bias filling for three different reverse 
biases. It is evident that for all these reverse biases, room temperature capacitance is 
almost identical indicating that C-V characteristics are flat around this temperature. This 
trap filling even at such high reverse biases indicates that trap occupancy is controlled by 
E/? movement with temperature. Note that such a behaviour is also in accordance with 
the temperature dependence of ^^tfl ^ shown in I-V characteristics of ion damaged 
Schottky devices. 

Thus the TSCAP has been an useful tool to study such defects with large concen- 
tration. These unusual features including the dominance of D1 defect in damaged layer 
motivates us to study of this defect in detail using TATS. 
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Figure 5.6: TSCAP signal for as-implanted Si taken under different reverse bias without zero 
bias filling. 

5.4 TATS Measurements 

In as-implanted (e.g. Ar'*' ) sample, because of large density of traps, nonexponentiality 
in transient is only to be expected. By performing the transient measurement in constant 
capacitance (CC) mode, we can avoid the nonexponentiality due to high trap concen- 
tration. Moreover, effect of series resistance on the capacitance transient can be severe 
causing ambiguity in conclusion regarding majority and minority carrier peaks in spec- 
troscopic signal. Though independent measurement of dissipation factor {u>RsC) leads 
to the conclusion of negligible effect of series resistance on capacitance measurement, we 
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altogether overcome this effect by performing measurement in CC mode. In this mode, 
a feedbac k tiu uit controls the capacitance of the device to a predefined set value and 
coriesponding voltage change is monitored as a function of time. Thus in CC mode, we 
measure the voltage transient instead of conventional capacitance transient. 

As discussed in Chapter 2, TATS signal is constructed from the capacitance/voltage 
transient by varying the time window keeping temperature constant. For a chosen cor- 
relation, peak position in signal gives trap emission rate at a particular temperature. As 
the emission rate varies with temperature, TATS signal peak position also moves. It is 
also important to know the sign of the signal in TATS to differentiate between majority 
and minority carrier peaks. In CC-TATS method, a majority carrier related emission 
results in a positive peak whereas in conventional TATS (CV-TATS) it gives rise to a 
negativ'e peak. Note that in conventional DLTS majority carrier related peak is negative. 
Henceforth, the terms positive peak and negative peak in TATS is to be understood with 
respect to particular methodologv' of analysis. Selection of suitable range of temperature 
for studying details of a particular defect is essential in TATS analysis. The detailed 
study of V 2 related peak will be discussed in the next Chapter where interesting features 
of this peak are revealed by performing transient experiments with varying filling times. 
In the next section, we will concentrate on the D1 level for improved determination of 
defect parameters. 

5.4.1 TATS of dominant peak D1 

In Fig. 5.7, CC-TATS signal is shown for transients taken at different temperatures for 
Ar"^ irradiated at two different doses. All the measurements are performed by keeping the 
depletion width (capacitance) constant corresponding to quiescent reverse bias. During 
filling, the device is taken to zero bias by special circuitry as described in Chapter 2. As 
the temperature is decreased, peak heights are reduced becoming undetectable at very 
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Figure 5.7: CC-TATS signal of defect D1 in as-implanted n-Si with (a) 5 x cm -2 (b) 
1 X 10^^ cm~^ At"'’ at various temperatures for fixed depletion width. 
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low temperature. An inspection of lineshape in CC-TATS spectra shows that the peaks 
are distoited on the right, especially peaks corresponding to high temperatures. The 
reduction in peak height can be attributed to the fraction of trap filled or emptied being 
lesser at lower temperature. Similar features were also observed in I-V, C-V, TSCAP and 
DLTS characteristics. A temperature dependence of transient amplitude can give rise to 
lineshape distortion in a temperature scanning method such as DLTS. This in turn does 
not allow proper estimation of trap emission rate. In such cases superiority of isothermal 
spectroscopic methods such as TATS are clearly established. In both the samples with 
different doses, this temperature effect is pronounced. 

As the defect is present in large concentration, this is expected to be responsible for 
carrier compensation in the damaged region which lies within the depletion region. As 
the measured depletion width is very large compaxed to the unirradiated sample, increase 
in depletion width is due to presence of high density of negatively charged traps (when 
occupied). Though these traps are located within the depletion layer, they are occupied 
due to unusual band bending at the defect location. 

5.4.2 Nature of nonexponentiality 

As mentioned in the previous section, lineshape of TATS peaks reveals that the peaks are 
not of proper symmetry as expected due to an exponential transient. For convenience of 
presentation from this point onwards, we would like to refer to peaks due to exponential 
transients as exponential peaks and that due to nonexponential transients as nonexpo- 
nential peaks. As discussed in Chapter 2, inspection of the lineshape reveals that the 
peaks shown in Fig. 5.7 are narrower than exponential peaks. Narrow peaks can arise 
from transients whose emission rate is increasing with time as opposed to the case of an 
stretched exponential (stretching factor 0 <. 1 ) transient where effectively emission rate is 
reducing with time. In no physical situation such narrowness is known to occur. Physical 
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peak is more distorted causing larger narrowness. 

However, questions may be raised regarding the uniqueness of such a fit where both 
peak height and time constant are treated as free parameters in the fitting process. In 
order to test for the exponential character of the portion of the signal prior to apparent 
peak, we resort to method of point by point analysis of transient to evaluate time constant 
as described by Mangelsdorf [144]. We indeed find that within experimental errors a single 
time constant is involved till the time of the sudden cessation of the transient. Once the 
time constant is fixed, the peak height is adjusted to match the lineshape of the left 
side of the peak which would be unique. The nature and quality of fitting to more than 
thousand experimental points for each transient is convincing enough to consider that pair 
of parameters obtained are stable points in the solution space of fitting. The parameters 
so obtained shows that the time constant obtained from apparent peak position and peak 
amplitudes are both underestimates. 

Figure 5.9 shows CC-TATS spectra with fitted exponential curves for two different 
samples implanted with different doses. Left peak marked Pi is due to divacancy related 
peak and it is seen to be broader than exponential peak. We will discuss more details 
about this peak later. The major peak (Dl) height is seen to be smaller in case of samples 
implanted with high dose and lineshape is less skewed. Figure 5.10 shows Arrhenius plot 
for the major peak (Dl) in four different samples prepared with different dose and 
annealing conditions. The estimated energy level and capture cross-section for the Dl 
defect are presented in Table 5.1. The energy level of the defect Dl is seen to vary between 
0.49 eV to 0.56 eV depending on processing conditions. The magnitude of capture cross- 
section is typical of a deep neutral center [145] accepting an electron to become negatively 
charged. The variation in emission energy can be attributed to the defects’ sensitivity 
to the degree of disorder in its environment or the degree of relaxation in the dcimaged 
layer. The relaxation of the damaged layer with relatively low temperature annealing will 
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Figure 5.9: CC-TATS spectra of D1 with fitted exponential peaks for two different samples 
implanted with different doses. Left peak marked Va is due to divacancy related peak and it is 
not fitted. 


Table 5.1: Trap parameters for the major peak (Dl) for different set of samples. 
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Figure 5.10: Arrhenius plot for the major trap (Dl) in Ar"*" implanted four different samples 
: low dose (1) as-implanted, (2) 160°C annealed and in high dose (3)as-implanted, (4) 160°C 
annealed n-Si. 

be discussed in a separate section. However, it is clear from the TATS spectra that the 
defect Dl stabilizes to deeper energy state for relaxed material in the case of low dose 
implantation or low temperature annealing. Prom fitting it has become clear that the 
peak height seen in DLTS or CC-TATS spectra is an underestimate. The actual peak 
height is shown with the fitted broken line in the CC-TATS spectra. For the spectra of 
Fig. 5.8 for 160°C annealed samples, the estimated trap concentration is 5.5 x 10^'* cm“^. 
This concentration is sufficiently large to cause nonexponential transient in conventional 
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meAsurenient. For high dose or as-implanted sample this concentration is smaller. This 
may be due to smaller extent of trap filling in those cases. Due to near intrinsic nature 
of the damaged laver, available carrier density for trap filling may be considerably low 
and thus a fraction of the traps show up in emission transient. Thus the observation of 
smaller peak in high dose implanted sample is an artefact of lesser extent of trap filling. 
Note that this concentration estimation is similar to that estimated from C-V hysteresis 
measurement or from I-V measurement. 

The occurrence of sudden cessation of the emission transient leading to premature 
termination of peak in DLTS and CC-TATS spectra can be understood as follows. On 
restoration of reverse bias immediately after the filling pulse, the depletion width increases 
considerably owing to large increase in negative charge in the damaged layer, thus pushing 
the Fermi level down below the trap level in the damaged region. This leads to exponential 
emission of carrier from extra trapped charge till the quasi-Fermi level hits the trap level 
quickly bringing emission transient to an end. Note that the phenomena is crucially 
dependent on the fact that the emitting center itself controls the quasi-Fermi level in the 
damaged region and the overall band bending through its occupancy. Hence, in general 
premature termination is to be expected in a p-i-n like structure where the intrinsic i- 
region is compensated due to a deep level. 

Since this phenomena occurs in time domain spectroscopies even in constant capaci- 
tance mode, it is a clear pointer to non-constancy of quasi-Fermi level in the region of large 
trap concentration in the damaged layer. It is the crossing of the quasi-Fermi level with 
Et that decides trap occupancy. As emission proceeds the split in the quasi-Fermi level 
also keeps changing till it meets the trap level in the damaged region restoring quiescent 
reverse bias occupancy. 

Observation of similar sudden termination of transients have been reported earlier in 
electron irradiated material and ascribed to reasons such as entropy driven metastability 
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Figure 5.11: CC-TATS spectra (solid line) of D1 measured at temperatures (Tl) 212.7K, (T2) 
217. 7K, and (T3) 228.6K. The dotted lines corresponds to peak fitted with exponential transients 
showing broad, exponential and narrow peaks respectively. 

[146] or recombination currents [147]. In our opinion, mechanisms responsible for anoma- 
lous phenomena such as these are best understood when studied in the time domain using 
isothermal transient spectroscopies. 

For this major peak, it is also observed that narrowness changes with temperature 
of measurement. Figure 5.11 shows the temperature dependence of the narrowness of 
this peak for Ar"^ implanted sample. With lowering temperature, peaks becomes less 
narrow and below certain temperature, it becomes broader than peak expected from 
exponential transient. We have already explained the cause of narrowness in peak shape. 
However, peaks broader than exponential is expected in such materials where defect are 
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Figure 5.12: CC-TATS spectra of low dose and high dose implanted Si. Peak D1 is fitted 
(dotted line) with Gaussian broadened energy of different FWHM. 

formed in a strained environment due to lattice damage. The degree of broadening in 
lineshape is a measure of degree of disorder in the immediate vicinity of these defects. A 
convenient measure of degree of disorder is in terms of spectroscopic lineshape broadening 
using Gaussian broadened activation energy of the trap level. CC-TATS measurements 
at lower temperature reveals a Gaussian broadening of trap energy with FWHM (full 
width at half maximum) of 25 meV for high dose (1 x 10^^ cm“^) .Ar'*’ implanted sample. 
This broadening is very small (6 meV) for low dose implanted sample. This is shown in 
Fig. 5.12 for as-implanted sample where experimental peak (solid line) due to Dl is fitted 
with a dotted line assuming Gaussian distribution in energy. Its sensitivity to strain and 
lattice disorder shows that it is possible to use it as a probe in the study of degree of 
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disorder for high dose implantations. 

5.5 Annealing of Ion Damaged Layers 

In order to understand the stability of the defects discussed above we tinncalcd at 400°C 
and 600°C the irradiated samples prior to device fabrication. In some samples, devices 
fabricated from as-implanted wafers were annealed at relatively low temperature of 160°C 
for different durations. We have already discussed the effect of furnace annealing on the 
C-V characteristics in Chapter 4. From the depletion width measurement it w;is found 
that defects in the damaged layer get annealed under these heat treatments resulting 
in reduced zero bias depletion width. Note that the zero bias width was ver\- large 
before annealing due to the presence of large number of negatively charged defects in the 
depletion layer of the junction. In this section we will examine the defect stability with 
respect to room temperature storage, relatively low temperature oven annealing and high 
temperature furnace annealing from spectroscopic measurements. 

5.5.1 Room temperature usage and low temperature oven an- 
nealing 

During the course of this work it was observed that for some high dose samples, there are 
minor, though noticeable, change in DLTS spectra with room temperature storage and 
usage. Figure 5.13 shows a set of DLTS spectra from high dose Ar'^ implanted n-Si. Curve 
1 and 2 in the figure show the spectra taken under the similar conditions on the same 
sample in an interval of two weeks. Though the V 2 peak seen to be exactly matching for 
the two spectra, the damage related D1 peak is shifted towards higher temperature for 
the same rate window. The room temperature was in the range of 20°-35‘’C for this period 
and intermittently transient and steady-state capacitance measurements were performed 
on these sample. The increase in peak height and shift in peak position indicates that 
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Figure 5.13: A set of DLTS spectra for (1) as-implanted, (2) after room temperature storage 
and usage of two weeks, (3) IGO^C, 30 minutes annealed, and (4) 160°C, 3 hours annealed 
samples. 
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Figure 5.15: Comparison of TSCAP signal from as-implanted, 400®C and 600°C annealed 
samples taken under similar conditions. The dotted curves refers to signal without applying zero 
bias filling pulse. 
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for proper estimation of trap activation energy, nevertheless an approximate estimation 
refers to a midgap trap level. Hence, the major defect created in as-implanted samples 
are not annealed out by 400°C annealing, rather they form more stable defects complexes 
upon annealing. By 600°C annealing these defects are almost annealed out as shown 
in Fig. 5.14. Though the spectra corresponding to as-implanted and 400°C annealed 
sample showed apparent increase in peak height of the major peak along with growth of a 
minority carrier related peak (H2), C-V measurements on these samples indicated gradual 
decrease of defect density with annealing at 400°C and 600°C. With higher temperature 
annealing, the C-V behaviour approaches to that of unimplanted sample. It can be noted 
that though DLTS spectra shows negligible majority carrier related peak height in 600°C 
annealed sample, C-V characteristic showed a substantial difference between unimplanted 
and 600”C annealed samples. This indicates that all the defects beyond the effective 
electrical interface are not completely annealed out by 600°C annealing for 30 minutes. 
It has been reported earlier that damage related defects begin to anneal by about SOO^C 
[149]. Our results are in agreement with this. The presence of positive peaks in DLTS 
spectra indicates that it may be due to minority carrier related trap levels. Note that the 
peak (H2) height is quite large in the spectrum. Though the definite origin of these are 
not clear at this moment, the appearance of minority carrier due to mild inversion in the 
damaged layer may cause such peak in DLTS spectra. 

We have also studied the defect levels present in annealed samples using TSCAP 
measurements. A comparison of TSCAP signal from as-implanted, 400°C and GOO^C 
annealed samples are shown in Fig. 5.15. The capacitance is monitored as a function 
of temperature in the heating cycle. The dotted curve corresponds to TSCAP signal 
without zero bias filling of the trap and the solid line corresponds to signal after zero 
bias cooling and subsequent heating at reverse bias. In as-implanted sample, difference 
between dotted curve and solid curve corresponds to trap filling. It is also seen that even 
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without application of zero bias, traps can be filled by simply lowering t einpruat ure due to 
shift in Fermi level. In 400‘’C annealed sample, a downward or neuaf ive st<‘p can be seen 
after zero bias filling of traps. This is again indicative of pre.semm of minority carrier trap, 
corresponding to the positive peak in DLTS. The change in capacitance in as-implanted 


sample is (juite large compared to other annealed samples. In GOO C iiimt aled sample, 
only a smaller negative step is seen without any majority carrier lelateii I positive) step. 


5.6 Possible Origin of Defect D1 

Though proper identification of a defect is rarely done on the basis of defect activation 
energy and capture cross-section, an attempt can be made on the basis of significant hints 
obtained otherwise from this study. Recall that this defect is present in a region much 
beyond the ion range predicted by TRIM calculation. According to TRIM simulation, the 
defect profile (mainly vacancies) lies close to the projected ion range. In our case, the ion 
range was expected to be about 1.23 iJ,m with a range straggling of 0.18 /<m. However, 
we have observed this defect at a depth of ~ 3 fim from C-V studie.s. lids is somewhat 
surprising that the defect profile can be so different from the TRIM jmcdicted profile. 
Though very few studies are reported in the literature regarding the defect profile in high 
dose unannealed ion irradiated Si, preliminary studies by others with low dose MeV ions 
have indicated that defect profile is broad in the deeper region whereas its concentration 
is reduced near the surface. In most cases, it has been attributed to vacancy diffusion 
or field enhanced migration of vacancy type defects from the irradiated zone toward the 
deeper region [23]. 

In earlier works of low dose MeV heavy ion implanted Si, this defect is not found 
[20, 21]. However, in GeV Ar"*" irradiated and annealed Si many defects levels apart from 
divacancy level are found to be formed [25]. Most of those defects were believed to be 
formed during post-implantation annealing. As our sample did not undergo any high 
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temperature post-implantation annealing, the dominant defect level observed by us has a 
different origin. 

The trap energy corresponds to a midgap level with capture cross-section tjTjical of a 
neutral center. Occurrence of midgap levels has sporadically been noted in the literature 
in ion irradiated samples without assigning any particular origin [29, 150, 151]. In some 
cases such defects have been attributed to dislocation related levels in samples annealed 
after high energy- implantation [25]. In some cases it has been referred to as damage 
related level without any further detailed studies [151]. In low dose implanted n-Si, 
electrical studies have normally revealed only vacancy-related defects namely oxygen- 
vacancy, phosphorus-vacancy, divacancy etc. These point defects are well characterized 
in the literature mainly in the study of electron irradiated silicon. Occurrence of other 
vacancy complexes are reported by Lee et al. through EPR studies [39]. In case of 
heavy ion implantation, point defect complexes are also likely apart from the above well 
known defects. Various defect reactions involving vacancies, interstitials and impurities 
lead eventually to more stable defects. As the vacancy and interstitials can migrate even 
at room temperature, survival of defect complexes are dependent on various factors. It is 
also known that many of the thermally stable defects may not be electrically active. 

EPR study on ion implanted Si reveals presence of centers related to dangling bond, 
di-interstitials, planar tetravacancies etc. [34]. Dangling bond centers are believed to 
be most dominant in high dose implanted amorphous sUicon [152, 153, 154). Note that 
in our sample one major defect level (Dl) is dominating the electrical signature. As we 
have worked with preamorphized silicon, it is tempting to attribute defect Dl to either 
Henglin g bond (D-center) or higher order complexes of intrinsic defects [156]. It has been 
shown recently that structural relaxation of amorphous silicon (a-Si) produced by MeV 
ion implantation is controUed primardy by point defect processes [69] and hence it would 
be interesting to investigate the role of accunrdation of Dl defects for higher doses in 
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creation, relaxation and eventual regrowth of such lasers. 

More recently it has been shown from molecular dynamif’ siirmLu i.nis and high res- 
olution electron microscopy observations that accumulation of i. ...• . i . .... (D-D) pair 
occurs in heavy ion implanted layers leading to homogeneous ainorphizatiftn bt'vond a 
threshold dose and only a defected material below that dose [uoj. Di-interstitials has also 
been studied by Lee et el. [40] and by Tan [56] in ion irradiattxl Si. Further. Privitera 
et al. [51] have recently demonstrated long range migration of Si self interstitials. More 
recently, molecular dynamic simulation in conjunction with .Monte Carlo simulation for 
defect process have predicted long range migration of defects and cluster formation even 
at room temperature [52]. In our case, since Ar+ damage is a source of seif-iiuerstitials 
[45], the involvement of interstitial complexes such as di-interstitials are strongly indi- 
cated. The migration energy of the di-interstitial has been estimated to be 0.6 eV [40], 
and it is therefore immobile species at room temperature. Recent cxju'riiut'nui! study 
on electrical signature of point defect to extended defect evolution upon liigh dose im- 
plantation has indicated also interstitial cluster related defects [19j. Collating all these 
information together, there is a strong indication that the dominant defect D1 is related to 
di-interstitial pairs or similar interstitial clusters. Since we observe prcvUmunantly point 
defect signature, the defect responsible is probably a point defect participating in a clus- 
ter. Though di-interstitial itself is only observed by EPR in the -f ve charge state, it may 
be driven to modified states, where it can bind electrons due to clustered environment. 
The fact that it’s trapping parameters are sensitive to damage environment in term of 
dose used and annealing also suggested that it is a deep level whose relaxation state is 
strongly determined by the state of clustering. The energy level brodening is only mod- 
erate showing that it is the long range component of the defect potential that samples 
only an average disorder in the environment. Also recall that we are observing deep levels 
in a thin slice of damaged layer at the ’eflFective electrical interface’ in our samples. The 
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sate of clustering ma.v be fairly homogeneous within this slice, though it may be varying 
with distance from the damaged layer. The changes with annealing would then be due to 
ripening of cluster containing point defects. 

The more definiti\e recognition with EPR would be helpful in such heaw dose im- 
planted samples. Ho'we\er, several conditions must be satisfied for this complimentary 
suppoits. The damaged layer must be having sufficient number of defects for detection 
and the clustering should not broaden EPR signal beyond recognition. Electron occupa- 
tion may have to be induced optically. Transient EPR studies analogous to DLTS would 
then be most convincing experiments provided one does detect EPR signals. 

Photocarrier lifetime measurement in ion-implanted a-Si have shown the presence of 
midgap states with capture cross-section an ~ 6 x cm^ [155] which is similar to 
the observed values shown in Table 5.1. The significance of defect D1 also lies in the 
observation that it is sensitive to strain and disorder in the environment. The observed 
disorder measured by Gaussian broadening parameter shows that we are still in the weak 
disorder limit. 

5.7 Summary and Conclusions 

In this chapter, we have studied the defects in the deep buried layers using a variety of 
junction spectroscopic methods. DLTS and TSCAP has been used as a general purpose 
survey technique to detect all possible defect levels in the band gap of Si in the tem- 
perature range 90K-320K. TATS measurements are made in constant capacitance mode 
to overcome problems due to high defect density, series resistance and strong tempera- 
ture dependence of capacitance. Results and conclusions based on these spectroscopic 
measurements on as-implanted and annealed samples can be summarized as follows . 

1. The DLTS spectra of the ion irradiated Schottky diodes with buried layers showed 
presence of peaks due to oxygen- vacancy ( A-center) , divacancy ( V 2 ) and a damage 
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related new defect (Dl) levels. 


2 . 


The DLTS and CC-TATS spectra is dominated !)>■ 
large concentration compared to V 2 related detes t. 


prcst'p.re nf It] in ven- 

Tlst* encri’v levels of the* defect 


Dl corresponds to a midgap level with capture cro.ss-.sectien tyiiieal uf a nentral 
deep level. The energy level position, varied from t).49 e\' to fl.r)2 e\'. is .sensitive to 
the degree of disorder and relaxation of the damaged layer. 


3. The dominant defect related transient specially at higli temperature shows prema- 
ture termination of the transient leading to skewed peak shape in DI.T.S or CC-TATS 
spectra. The unusual capacitance spike during trap filling in the damaged layer is 
also associated with such observation. 


4. Nonexponentiality of the transient was analyzed using CC'-'rATS spect ra and showed 
that the measured peak height and time constant both are underesiimaied in the 
measured spectra. The process of sudden termination in transient is explained 
on the basis of quasi-Fermi level movement due to trap occupation change during 
emission. 


5. Spectroscopic lineshape analysis showed energy brr)adeniiig of about 25 me\' for the 
major peak due to defect Dl for high dose implaihaMi samples. Negligible broadening 
of lineshape was present for the low dose implanted sample. 

6. Low temperature heat treatment of ion damaged sample showed relaxation of the 
damaged layer causing defect energy to go deeper in the band gap and more defect 
was accessible to the transient measurements. 

7. Relatively higher temperature furnace annealing showed gradual change of the defect 
configuration indicating formation of more stable configuration of the same defect. 
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Most of these defects could be annealed out by vacuum annealing at 600°C for 30 
minutes. 

8. No electrical signature of extended defect formation is seen from arm paling studies 
as doses used were below amorphization threshold and annealing temperatures were 
lowc 

9. The possible origin of the major defects (Dl) are discussed on the basis of self- 
interstitial migration and formation of complexes like di-interstitials or their clusters. 
Possibility of dangling bond related center for this defect is not completely ruled out 
as specific orientation the defect complexes may give rise to unsaturated dangling 
bonds. 

10. Minority carrier related peak is detected on annealed samples from DLTS and 
TSCAP studies. Occurrence of such a peak may be due to possible inversion in 
the damaged region causing minority carrier available for trap filling. 

11. The detected majority carrier defects are seen to be present much beyond the pre- 
dicted damage profile from TRIM calculation. 
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the Fermi level governs the energies at which formation of D-defects occur (157, 158]. 
Alternative approaches have been to explain these differing energies by defect relaxation 
mechanism or even simply by occupation statistics of pre-existing distribution. Recently 
junction capacitance transient measurements or in other words, observations of charge 
relaxation mechanism have been used to seek the nature of defect relaxation in a-Si. A 
study of capture of carriers at gap states led Cohen and coworkers [159] to propose a 
novel hieraichical defect relaxation mechanism in which polarization around negatively 
charged detect cause relaxation driving its emission energy deeper into the gap. However, 
an alternate explanation of their transient data was proposed by Farmer and Su [160], 
who invoked a parallel mechanism of charge redistribution from shallower to deeper levels 
through multiple trapping at exponentially distributed gap states. In either case compar- 
ison with experimental transients required convolution with broad density of states and 
hence model distinction remain equivocal. The class of mechanisms in which a hierarchi- 
cal relaxation is invoked can be viewed as serial mechanisms, whereas in contrast parallel 
mechanisms invoke participation of independent defects coupled through some common 
constraint. This difference in viewpoint has been central to many controversies in the 
defect literature, specially in the case of disordered semiconductors and compound semi- 
conductor alloys. Controversies regarding distinction between hierarchical and parallel 
mechanism is common to many disorder systems and a succinct critique in this regard 
has been developed by Palmer et al. [161]. 

This brief introduction to the wider debate in disordered systems has been necessary 
since vve present in this chapter unambiguous evidence of occurrence of both parallel 
and serial mechanisms of charge relaxation in the damaged layers in different regimes of 
sample conditions. We show that the study of defects in preamorphized silicon provides 
important conclusions about charge relaxation mechanism in amorphous silicon as these 
materials can be treated as precursor to the amorphous phase. It is possible to study 
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a variety of defect related phenomena more convenientH in i tjntrolli.ci disordered 

materials without having to take recourse to complex dccomolution prut udurc to compare 
theoretical predictions with observed experimental features. Apart from this, these studies 
are, of course, intrinsically important in understanding the nature of disorder in ion 
implanted material. 

6.3 Design and Principle of Experiments 

Experimentally our motive behind the measurements that we describe in this chapter was 
to monitor relative trap occupancies of multiple trap with different filling time duration, in 
order to investigate the nature of coupling between the trap levels observed in irradiated 
Si. In conventional cases where more than one deep level is present in the material under 
study, the occupancies of the deep levels (i.e as monitored by the magnitude of the TATS 
peak) increase monotonically with filling time and saturate for large filling times indicating 
its concentration. The sample is taken to zero-bias for a predetermined certain duration 
ranging from 100 fj,s to few seconds, and the resulting transient is monitored. One set 
of measurements constitute use of many different filling times at a fixed temperature. 
The temperature is chosen such that both peaks (peak V 2 and Dl) occur within chosen 
time scale of measurements. In order to unravel any possible coupling kinetics of capture 
and to achieve sufficient resolvability, we resort to TATS. The implementation of higher 
order spectroscopy in transient analysis allows higher resolvability between closely spaced 
peaks. The choice of order of spectroscopy is guided by the required compromise between 
resolvability and signal to noise ratio. 
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Figure 6.1: Second order TATS spectra of peak PI for filling times varying from 100 yus to 15 
s. Arrows indicate the apparent peak position for two different centers PlA, PIB. 

6.4 Trapping Kinetics At Low Temperatures : Evi- 
dence For Charge Redistribution 

6.4.1 Occupancy change with filling time 

Since the divacancy peak (V 2 in the literature has been associated with two different 
defects, we use TATS in order to resolve possible multiplicity in its lineshape. Moreover, 
for verv small filling time and large filling time, we found that the divacancy related peak 
(V 2 ) is apparently shifted. To resolve any closely spaced peaks which may cause such 
phenomena, we perform higher order TATS analysis on the emission transients due to 
different filling times. Figure 6.1 shows second order CC-TATS spectra of defect corre- 
sponding to peak V 2 for several durations of filling time. For convenience of nomenclature, 
we denote this peak as PI in this chapter. This clearly shows that the lineshape is consti- 
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tuted of contributions from two distinct emitting centers for low filling times, And becomes 
progressively dominated by the slower peak for large filling times. This also shoves that 
changes in occupancy of the constituent peaks with filling time is nonnionotonic. An 
unambiguous demonstration of the two peak nature of this spectra is clearly central to 
the understanding of carrier capture kinetics. 

To quantitatively evaluate occupancy changes of all the three peaks PlA, PlB and Dl, 
we use first order TATS spectra for filling times varying over five orders in magnitude. 
A typical set of such spectra is shown in Fig. 6.2 (frames (a) to (e)) for progressively 
increasing filling times. Note that peak Dl initially increases primarily at the expense 
of peak PlA and then saturates. The relative concentration of the constituents have 
been evaluated by fitting contributions from three different defects. The emission time 
constants of PlB and Dl have been obtained from their respective peak positions in the 
spectra in the high filling limit. The emission time constant (tc) of PlA is then estimated 
from experimental spectra. All time constants (rg) are kept fixed for varying filling times, 
while the relative concentrations of the peaks are treated as free parameters in the fitting. 
In the frame e of Fig. 6.2, a fitting to the peak Dl is also shown by dotted line. This 
required a Gaussian broadening with FWHM (full width at half maximum) of 25 meV in 
energy. However the amount of broadening needed for the same peak for the lower dose 
(5 X 10^^ cm~^) sample was 6 meV. Peaks PlA and PlB show similar small broadening 
in the case of high dose samples. Similar feature of variations in trap occupancies with 
filling is also shown in Fig. 6.3 for low dose implanted sample. Larger peak heights are 
observed for this sample. 

The relative occupancies obtained from the analysis have been plotted as a function 
of logarithm of filling time in Fig. 6.4(a) and (b) for samples irradiated with two doses. 
The occupancies of peak PlA initially increase going through a maximum, significantly 
enough, at filling times approximately equal to its emission time constant. The increase 
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Figure 6.3: First order TATS spectra at 211.6K for diflferent filling times {tf) showing pro- 
gressive changes in occupancy for low dose (5 x 10^^ cm”^) implanted sample. Peak PI is 
fitted (dotted line) to two different centers PIA and PIB. Peak D1 can be fitted to a Gaussian 
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Figure 6.4: Occupancy of individual peak as a lunciion m o 

^ \ C , ini3 nm-2 and Cbl 1 X 10^^ cm“^. Arrows mdicate the 

different doses of Ar"*" ions : (a) 5 x 10 c I / 

measured emission time constants for the corresponding cent 



Trapping kinetics of defect states in n-Si 


164 


of peak D1 occurs when PlA starts decreasing. Peak PlB also goe.s through a nijiximum 
and its contribution is small to the growth of peak Dl. These broad feattires are common 

to both the samples irradiated to different doses. 

Clearly this is a case of coupled carrier kinetics with multiple trap.s. There are a prion 
many possible mechanisms which can lead to such a coupling viz. nniltiple charge states 
of the same defect, internal conversion from a metastable to stable state or even some 
more elaborate defect relaxation process. The class of mechanisms in wltith hieiarchical 
relaxation [161] is invoked can be viewed as serial mechanisms, whereas in contrast 
parallel mechanisms invoke participation of independent defects coupled tlirough some 
common constraint. This difference in viewpoint has been central to many controversies 
in the defect literature, specifically in the case of disordered semiconductors and compound 
semiconductor alloys [162]. 


6.4.2 Simulating experimental features 


In the present case, the key observation is that the shallower state.s lo.se ehargc.s to the 
deeper ones, and that too for filling times larger than their characteristic emission time. 
This in itself pins down the mechanism to redistribution of charge.s througli multiple 
trapping. For filling times (tf) shorter than time constant of emission (r^,) of a particular 
level, the occupancy increases in proportion to the product of its capture constant and 
number density of empty states. For tf > r^, emitted electrons from a shallow trap is 
recaptured by deeper states with increase in filling time. The situation can be modeled 
simply by following rate equations for independent traps, 


driTi 

dt 


= -eiriTi + Cin[NTi — nn), i = l to m 


( 6 . 1 ) 


with the crucial constraint on number of electrons is, 

m 

" = - nr.), 


( 6 . 2 ) 



167 


Trapping kinetics of defect states in n-Si 

defcTt: for Inghor filling time D1 dominates. Note that in the e.xperimental data, though 
the int<.r„u.diat,. level PIB participates in the process of charge redistribution, its peak 
occiipimcv occurs for filling time larger than its emission time constant. This may be due 
to fieki enhancement of emission rate during the transient. 

6.4.3 Discussion 

The factors that made possible such an unambiguous observation of charge redistribution 
is worth listing and each of them has an important bearing on both phenomena and design 
of experiments. 

• Firstly, a necessary condition for observation of the parallel mechanism of redistri- 
bution of electrons through conduction band is that the number of trapping sites 
available is much larger than the number of electrons available for capture. This 
condition is clearly satisfied in our experiment in the region of large damage. We 
not only have large number of traps, but also the region is converted to highly com- 
pensated region. Application of bias across the sample lets us control the number of 
carriers in the disordered region conveniently by raising or lowering the Fermi level 
in the crystalline region. 

• The study of an embedded disordered region is also significant in another respect. 
The change in charges in the compensated region reflect in the width of the depletion 
layer, whose edge lies in the undamaged region. This allows convenient use of 
standard high frequency capacitance measurements to monitor charge relaxation in 
the damaged region. 

• Also note that the time scale of dynamics varies from milliseconds to several seconds. 
For divacancy level whose capture cross-section is 3 x 10"^® cm^ the filling time m 
normal case with carrier concentration of 1 x 10^^ cm-^ would be of the order of 
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few tens of nanoseconds. In the damaged region however, the carrier concentration 
is down to near-intrinsic levels and hence the dynamics is considerably slowed down 
enabling clear observation of charge redistribution effect in a convenient time scale. 

Charge redistribution among multiple traps have been observed in case of silicon re- 
lated DX centers in AlGaAs [163, 125]. However, the dynamics is more complicated 
[125] in that case due to specific model of local alloy disorder, in which the same donor 
atom relaxes in a variety of local environments with distinct energ}'’ levels. In case of 
hydrogenated a-Si, model distinction was difiicult since both the serial model of Cohen 
et al. [159] and the parallel model of Farmer and Su [160], used convolution with a 
broad density of states. Our results unambiguously confirm that it is not necessary here 
to invoke hierarchical relaxation to explain progressive deepening of energy levels leading 
to slowing down of emission transients with filling time. In our study we have been able to 
simulate conditions of those experiments without going over to the amorphization regime 
and thus avoiding the complexity of having to deal with any particular distribution of 
continuum of levels. 

Note that the fully developed D1 peak was non-exponential to the extent that its 
density of states required a Gaussian broadening with FWHM of 25 meV as pointed out 
in Fig. 6.3e. On the other hand no significant broadening was observed for the same level in 
the sample irradiated with the lower dose. This broadening therefore provides a measure of 
inhomogeneities or disorder that the defect samples in its environment. It would require 
a more detailed study to track this level of broadening along with possible increase in 
multiplicity of levels to the full amorphization limit. It is important to point out the 
central role that a careful choice of spectroscopy plays in distinguishing between these two 
different effects of disorder i.e. increase in multiple levels and degree of energy broadening 
of these levels. The mean emission energy is the result of high degree of localization in 
real space provided by the short range part of the defect potential, whereas its broadening 
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can be viewed as resulting frcm the long, possibly Coulombic tail of the potential which 
sample s a much larger volume of its environment. The existence of this broadening in our 
spectra foi high dose sample can be considered as a post facto confirmation that we are 
indeed in the high disorder limit. 

The defect level corresponding to peak PIB is located at — 0.42 eV with capture 
cross-section of ~ 3 x 10 cm^ as obtsdned by CC-TATS analysis at different temper- 
atures. It is commonly attributed to divacancy V 2 (-/O). The peaks PlB is due to the 
divacancy related level. It is interesting to note that effects such as motional averaging of 
two configurations and sensitivity to damaged related strain have been invoked to explain 
DLTS lineshapes of (-/O) related peak [20]. The broadening of the individual peaks 
TATS spectra can be attributed to lattice strain in the damaged layer. Indeed it was 
found from DLTS study of hea\w ion irradiated n-Si that the divacancy related defects 
with energy Ec-0.42 eV was formed in a region of lattice accommodating more strain. 

The origin of the Dl defect is discussed in the previous chapter. It is confirmed from 
this study that these are responsible for carrier compensation in the damaged layer and 
less number of free carrier makes the carrier dynamics considerably slow to monitor such 
phenomena. 

In summary, we have studied charge relaxation in severely damaged regions in silicon 
by embedding them in the depletion layer of otherwise crystalline material. Multiple trap 
levels, which have moderate broadening of energy for heavy damage, have been resolved. 
We have provided unambiguous transient based spectroscopic evidence of the occurrence 
of charge redistribution among multiple levels in the compensated regions. 
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6.5 Trapping Kinetics In Relaxed Materials : Evi- 
dence For Novel Defect Relaxation 

During the course of filling time experiments, a novel phenomena was consistently ob- 
served in samples with either low dose implantation or those partially relaxed after an- 
nealing. It was most strikingly observed in (i) low dose Au"*" (5 x 10® cm“^) implanted 
n-Si and (ii) Ar"*" implanted samples in which partial relaxation was induced by annealing 
at 400‘'C for 30 minutes. The presence of a preamorphized region with comparatively 
lesser degree of disorder is common to these set of samples. 

Figure 6.6 shows CC-TATS spectra at room temperature for a series of filling times 
(t/) varying from 20 ms to 1 minute. The feature of decreasing signal at the lowest 
emission of times (ms range) is due to the falling portion of the peak corresponding to the 
major trap studied in the previous chapter. Apart from that, the well formed peak shown 
in the figure is due to a new distinct peak whose time constant is strongly dependent 
on the filling time, becoming slower for longer duration of filling pulses. Note that the 
time constant of emission changes 4-5 orders of magnitude as does the filling time. Also 
note that the height of the new peak is approximately similar for all different filling times 
showing that the corresponding defect is located in the same region of peak and has nearly 
the same degree of occupancy irrespective of duration of pulse. Though the figure shows 
only a selection of curves, the time constant was observed to be continuously changing 
even when monitoring with small changes in filling time. In other words, the emission time 
constant of carriers from the defect progressively increased with filling time and hence this 
is not a case of a defect relaxing through discrete states in energy. At the same time, the 
lineshape analysis shows that the peaks are approximately exponential for lower filling 
times becoming slightly broad for largest filling time. This shows that for any particular 
filling time, we are observing emission from a discrete defect level. Had it one been from 
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Figure 6.6: CC-TATS spectra for different filling time (t/) at a fixed temperature showing 
evolution of a new peak whose time constant increases with increasing filli ng tunes. 
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Figure 6.7: Change of emission time constant (r) with filling time (tf) for the evolved peak of 
Fig. 6.6. This shows stabilization of the relaxed state energy (Ejpc ln(r)) for longer filling times. 

a part of broad density of states in energy, we would still have observed broadening for 
all the peaks. The peak evolves from a time constant of 6.5ms to 14s for filling times 
in the range 20ms to 59s at a particular temperature. This three orders of change in 
time constant would translate to about 180 meV change in terms of activation energy. 
Figure 6.7 shows log of time constant as a function of filling time. This is indicative of 
occurrence of hierarchical relaxation through continuous metastable states, stabilizing in 
a state with the longest time constant. This gradual slowing down (with r and tf obeying 
a power law) of emission or deepening of energy is typical of hierarchically constrained 
dynamics [161]. The emission time constant is dependent on resident time of electron in 
the defect site. This dependence on history and power law dependence are signature of 
non-Markovian process possessing a memory of the duration the electron spends at the 
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Figure 6.8: CC-TATS spectra at different temperature for a fixed filling time (t/=ls) showing 
temperature independence of the evolved peah. 

occupied defect [159]. 

Another intriguing feature of this phenomena is that the time constant is nearly tem- 
perature independent for any particular filling time, t^^. This is illustrated in Fig 6.8 where 
TATS spectra for a filling time of Is are shown for different temperatures varying over 
30K. The time constant of major peak on the left of the figure does indeed change as 
would be expected of a normal thermally activated process. However, the peaks that 
evolves with t/, does not show shift in time constant with temperature. This shows that 
the process associated with the relaxation is nearly athermal. Prom the athermal nature, 
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Figure 6.9: Power law dependence of r and tf for spectra shown in Fig. 6.6. Solid line refers 
to a best fit to relation : r ~ (l+ty/te//)“ . 

we conclude that the part of the free energy changing with defect relaxation can be con- 
sidered to be proportional to temperature T. Incorporating the power law dependence of 
filling time (tf) and emission time constant (r), the effective energy associated vith the 
relaxation can be written as 

Ueff ~ KbTln[{l -f- tf/teff)°‘] (6.7) 

where it is characterized by an effective time constant of relaxation t^// and exponent a. 
In Fig. 6.9, we show such a dependence of r and tf where the linear portion of the plot can 
be fitted with te//=0.15 ms and 0 = 1 . 2 . The nature of this energy term also suggests that 
it can be viewed as an entropic term contribution to the measured enthalpy [146, 164]. 
Physically it would mean that traps undergoing relaxation during electron occupation, 
have a distinct entropy most probably owing to changes in degeneracy factor. 
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Trapping kinetics of defect states in n-Si 
be dependent on the nature of damage which may allow multiple final configurations for 



log (time/s) 


Figure 6.10: CC-TATS spectra with varying filling time at a fixed temperature for Au+ irra- 
diated (5 X 10® cm~^) n-Si. Multiple peak structure in evolved peak for longer filling times can 
be noted. 

long term filling of trap. 


In summary, in ion damaged deep buried layers evidence for charge redistribution 
among multiple traps have been shown to occur through trapping kinetics studies over 
5-6 orders of magnitude in time. At relatively higher temperature, for partially relaxed 
samples, a novel defect relaxation process is observed where defect energy progressively 
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deepens with trap filling time and this process is nearly athermal. It is important to reit- 
erate here that isothermal spectroscopy such as CC-TATS is absolutely essential in study 
of phenomena such as described in this chapter. In a temperature scanning spectroscopy, 
such as DLTS, sensitivity to filling time and its temperature dependence on one hand, and 
possible occurrence of athermal process on the other can completely distort manifestation 
of such phenomena beyond recognition. 



Chapter 7 

Defects In Ion Irradiated p-Si 

The p^e^ ious three Chapters were devoted exclusively to n-type silicon. This chapter, in 
a nutshell, is devoted to similar studies on p-type silicon. The scope is limited to explore 
major differences and similarities with n-type case hoping to identify commonality of origin 
of traps and related phenomena. As pointed out in Chapter 1, a perusal of literature shows 
that irradiation induced defects has been studied to a lesser extent in p-type Si, in general. 
It is specially true for MeV heavy ions, where most of the studies have concentrated on 
evolution of defects observable by cross-sectional TEM with annealing. We focus on 
electrical characterization using C-V and deep level spectroscopies in samples irradiated 
with MeV Ar"*" and Au'*’ ions. During the course of this work we have used both epitaxial 
and CZ p-type wafers. 

We will keep this chapter brief by avoiding repetitions of explanations of similar effects 
already discussed in n-type case, and by focusing only on features specific to p-type 
samples. 

7.1 C-V Characteristics 

A number of Schottky diodes were made on unimplanted and implanted wafers. The 
control diodes (made on unimplanted wafer) are tested by C-V measurements at various 
temperatures. The plot of vs. V was seen to be linear indicating uniformity of 
shallow doping (boron) with concentration 1.5 x 10 ctti in the high voltage region. 
Figure 7.1 shows l/C^ vs. V plots for Ar+ ion irradiated p-Si at several temperatures 
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Figure 7.1: Temperature dependence of C-V characteristics of Ar"^ irradiated p-Si Schottky 
diode. 

showing qualitatively similar features as in n-type case. Notable features that are common 
to in 1/C^ vs. V plots are : 

• occurrence of flat region indicating large trap concentration at the crossing point 
between Fermi level and dominant trap level, 

• transition to linear region at higher bias showing carrier concentration at the edge 
of the depletion layer. 

• increase in depletion width with lowering temperature due to increased trap Ailing, 


• increase in depletion width for higher doses 
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• o( c.nri( nee of depletion layer edge much beyond the range of the ions and hence the 
damage location. 

All thesf' fiMtures indicate, as explained for n-type analysis, the occurrence of defect 
migration and clustering. The differences are only quantitative in nature. These are 
listed biilow on the basis of studies on several batches of samples : 

• zero bicis depletion widths are about 60% of that of n-type case though the doping 
concentration are similar. 

• the transition from flat to linear region in 1/C^-V cun'^e is consistently less sharp 
than n-type case or with non-linearity extending over larger bias range. 

Figure 7.2 shows phenomena of hysteresis in C-V in close parallel to n-type case. This 
again indicates involvement of a dominant trap whose occupancy mainly controls C-V. The 
combinations of temperature and sweep rate needed to observe hysteresis are distinctly 
different from n-type case. 

7.2 Deep Level Spectroscopies of Defects 

7.2.1 Thermally Stimulated Capacitance (TSCAP) Measure- 
ments 

As described in Chapter 2, TSCAP is an useful spectroscopic technique for detection of 
deep trap with high concentration. To And dominant traps in ion damaged buried layers 
in p-Si, we measured the capacitance as a function of increasing temperature after the 
traps are filled with zero bias pulse at low temperature. Figure 7.3 shows a differentiated 
TSCAP signal for Ar^ ion implanted p-Si. Corresponding to the steps in capacitance 
values, two major peaks are detected in the temperature range 120K-230K. Though the 
detection sensitivity of TSCAP method is not so high as other transient spectroscopic 
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Figure 7.2: Hysteresis in C-V characteristics for decreasing and increasing reverse bias. Arrow 
indicates direction of voltage sweep. 

methods such as DLTS, it was used as a survey technique to find the dominant traps. 
The trap parameter estimation was done from the DLTS and TATS methods. These 
parameters were used to simulate TSCAP in close agreement with experimental spectra 
using Eqn. 2.9 of Chapter 2. 

7.2.2 Deep Level Transient Spectroscopy (DLTS) Studies 

No traps could be detected using DLTS in control samples upto detection limit of 10“^ 
of background doping. This ensures that any traps detected in irradiated sample would 
be irradiation induced. The irradiated samples used for this study was irradiated prior 
to Schottky diode preparation i.e. irradiation was performed on bare silicon wafer pieces 
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Figure 7.3: Differentiated TSCAP signal after zero bias filling of traps in Ar+ ion irradiated 
p-Si. Heating rate was 2.0K/min. 


%vith back ohmic contact made by thick aluminum deposition. Devices were made on 
these wafers after proper etchiug/cleaning of surface oxide or any other contaminants. 
A typical DLTS spectrum tor sample irradiated with 1x10'* Ar+ ions/cm" is shown m 
Fig. 7.4. For the chosen rate window ((2 ms)"'), three majority carrier related peaks 
(PI. P2. P3) and one minority carrier related peak (P4) are observed in the temperature 
range 120K-310K. The smaU peak PI has an activation energy of 0.37 eV and a capture 
cros^section of 4 x IQ-'* cm= as estimated from the Arrhenius plot shown in Fig. 7.5. 


r,., , J + c Trnlp tran commonly attributed to CiOi complex [171] in CZ 

This peak corresponds to a hole trap commo y 

, ■ 1 • nj 17 riei+ iTTinlflnted o-Si as well. This trap is observed 
substrates. We found this level in MeV Au unplan p 

in small concentration in epitacdal layers since oxygen content in these laymrs are small. 
The line shape of the DLTS peaks are Sued with peaks due to three exponent, al tran- 
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Figure 7.4: Typical DLTS spectrum of Ar'^ ion irradiated p-Si. Dotted line is a theoretical 
fit with three exponential transients corresponding to the three peaks Pi, P2 and P3. Time 
constant corresponding to the chosen rate window is 2 ms. 

sients as shown with the dotted line in Fig. 7.4. The trap parameters were chosen from 
Arrhenius plot obtained for each peak on the basis of DLTS results. Due to inherent 
nonexponentiality, Arrhenius plot may not yield accurate estimation of trap activation 
energies (Ex) and capture cross-sections (<7^) and therefore, ap value was adjusted so as 
to get the peak position at the same temperature as seen in the experimental spectra. It 
is clear that major peak (P2) line shape is broader than the standard line shape expected 
from exponential transients. Such broadened line shape would correspond to a Gaussian 
broadened activation energy for the corresponding trap level. The line shape distortion 
may be due to several reasons. Since the concentration of traps is expected to be high 
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Figure 7.6: DLTS spectra of Ar'^ ion irradiated p-Si plotted for diJfferent rate windows (Sn): (a) 
(2.2ms)~^, (b) (5.1ms)~^, (c) (15.2ms)~.\ (d) (30.1ms)~\ (e) (75.0ms)~^ and (f) {200.0ms)“^ 

corresponding to peak P3 are not possible since the peak is not so well resolved. 

Figure 7.6 shows a set of DLTS spectra for different rate windows. The reduction in 
peak height for slower rate windows for the major peak (P2) is substantial and presence 
of feature P4 distorts the peak P2. It is known that a strong temperature dependence in 
capture rate and series resistance [101], can lead to such features in DLTS spectra. In 
our case strong temperature dependence of capacitance is primarily responsible for peak 
height reduction as well as nonexponentiality in DLTS peaks. Due to the above reasons 
DLTS results have not been used for trap parameters estimation, rather they have been 
used as a survey technique for detection of different traps. Proper evaluation of the trap 
parameters have been made using time analyzed transient spectroscopy which is presented 
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Figure 7.7: CC-TATS spectra at various temperatures showing three peaks {P2, P3 and P4) 
for Ar"*" ion irradiated p-Si in a convenient range of time scale. 

as the problems due to high trap density and series resistance are sources of distortion in 
conventional TATS spectra, such conclusions may be in error. By using the technique in 
constant capacitance mode, it is possible to avoid the expected nonexponentialities due 
to large defect density in constant capacitance voltage spectroscopies [174]. Moreover, 
constant capacitance transients are not distorted due to any dc resistance present in the 
device. 

Figure 7.7 shows a set of CC-TATS spectra at several temperatures for Ar'^ implanted 
p-Si. In the range of time scale shown in the figure, we see two majority carrier peaks (P2, 
P3) and the minority carrier related peak P4.^ The minority carrier peak appears as a 

^Note that a majority carrier peak in DLTS appears as a negative peak in conventional spectra whereas 
it appears as positive in CC-TATS spectra. 
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fulh (if lopt.d feature in contrast to DLTS spectra. The peak corresponding to peak PI 
in DLTS .spi'ctra is too fast to appear for these temperatures. The lineshape analysis of the 
niajoi pe.ik (I 2) in CC-TATS spectra for low temperature transients show that the peak 
is bi Odder than « oi responding exponential transient. Thus lineshape broadening is indeed 
due to inhomogeneous environment. However, CC-TATS peaks for higher temperature 
transients are seen to be narrower due to premature termination of transient as it was 
in the (;ase of n-Si. But, parameter extraction at this temperature through fitting is not 
possible due to presence of another close peak slower in time constant. The reason for 
premature termination is same as was discussed in case of n-type Si. 

Figure 7.8 shows Arrhenius plot for the major peak whose emission rates have been 
obtained from CC-TATS spectra using low temperature transients. For comparison, time 
constants obtained from conventional DLTS spectra which results in non-linear Arrhenius 
behaviour, is also plotted in the same figure. Data from CC-TATS is most reliable and 
gives an activation energj' of -1-0.52 eV with an unusually high capture cross-section of 
9.9x10“^'* (rrn^. In Au^ implanted p-Si, we also observe the same damage related trap. 
Good match in Arrhenius plot for this major trap (P2) in Ar'^ and Au^ implanted p-Si 
has been obtained and shown in Fig. 7.9. For Au^ implanted sample, activation energy 
and capture cross-section was found to be Ev+0.52 eV and 7.9x10 cm^, respectively. 
A defect having similar parameters have been observed in heavy ion damaged Si and has 
been attributed to damage [21, 25]. However, such a level is not detected in damage 
created with light ions. Such high value of capture cross-section for this defect may 
be due to coulomb assisted trapping at the defect site. Similar high values of capture 
cross-section have been reported for grain boundary defects in Si [175]. 

As in the case of n-type Si, it is not possible to obtain the concentration profile of 
this defect since trap information is obtained only at the location of Er - -Ef crossing. 
However, it is quite clear that the defect concentration is as high as in the case of n-type 
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Figure 7.8: Arrhenius plot for the major peak (P2) obtained from CC-TATS and DLTS mea- 
surements. Solid line is a least square fit to the data from CC-TATS measurements. 

at the effective electrical interface and a conservative estimate^ would put it as 3x10^® 
cm~®. 


7.3 Capture Kinetics Of The Major Trap 

7.3.1 DLTS and TATS of capture transient 

The capacitance transients during capture are observed to be surprising slow. Capture 
transients are recorded at various temperatures and analysis shows that they are expo- 
nential in nature. As suggested by Ghosh et al. [176], for a straightforward measurement 

^For estimation of trap concentration, total amoimt of defect charge corresponding to amplitude of 
transient is calcnlated cind then divided by the approximate volume of the buried layer estimated from 
Gaussian damage profile. 
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Figure 7.9: Arrhenius plot for major peak (P2) in Au'^ and implanted p-Si employing 
CC-TATS measurements. 

of capture barrier, DLTS scans using capture transients were taken in the temperature 
interval of 175-300K. Typical capture DLTS peaks for different rate windows are shown 
in Fig. 7.10(a). Capture time constants have been obtained by both DLTS and TATS 
processing of capture transients and have been plotted in an Arrhenius plot in Fig. 7.10(b) 
yielding an activation energy of the 0.66 eV. Note that this is an unusually high activa- 
tion energy for capture processes and this value is larger than the barrier for emission. 
Therefore, this barrier cannot be associated with the capture cross-section of the defect 


alone. 
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Figure 7.10: (a) DLTS spectra for capture transients with rate windows : (1) 500 (2) 

125 5“^ and (3) 50 s~^ (b) Arrhenius plot obtained by DLTS and TATS analysis of capture 
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Figure 7.11: Filling time dependence of the occupancy of the major trap in Ar'^ implanted 
p-Si using varying pulse width technique in conjunction with CC-TATS 

7.3.2 Variable filling pulse- width technique 

For a more accurate measurement of the capture barrier, we performed the transient 
(capacitance/voltage) measurement using variable filling pulse width technique [177]. 
Figure 7. 11 shows the filling time dependence of occupation of the major trap (peak P2) 
at 216K determined using constant capacitance varying pulse width technique. Capture 
is seen to be exponential from the point defect nature of the traps. However, for the 
same sample using constant voltage transients, trap filling during capture is observed to 
be slow and approximately logarithmic in time. Hence, we conclude that such slow filling 
is an artefact of nonexponentiality due to large trap concentration. Similar logarithmic 
time dependence of capture has been reported from conventional DLTS measurements 
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for plastically deformed silicon with dislocations [1/3|. Therefore, a (,on\entional DLTS 
study can mislead one to conclude that defects are extended defects with time dependent 
capture barrier. A time dependent capture barrier will give rise to logarithmic filling time 
dependence as is the case for extended defects. In a recent study, Benton et al [19] 
have reported in p-type ion irradiated sample a level at Ei,+0.48 eV using conv'entional 
DLTS. They identify it to be due to extended defects {311} principally on the basis of 
logarithmic filling time dependence of occupancy. However, a more careful study of this 
level is required to resolve whether this is an extended or point defect. 

In our samples annealed at GOO^C, the major trap was found to have an activation 
energy of E^+O.Sl eV with capture cross-section of ~ 10"^^ cm^. The peak related to 
this level in DLTS was found to occur at a lower temperature than the above level (for 
the same rate window) reported by Benton et al.. 

Constant capacitance transient measurements using the variable filling pulse width 
technique at different temperatures reveal strong temperature dependence of capture ki- 
netics of the major trap. Figure 7.12 shows an Arrhenius plot of capture time constant 
yielding an activation energy of 0.66 eV, which is unusually high. Note that this value 
is identical to that found from direct measurement of the capture DLTS/TATS. Since 
the activation energy for the emission is smaller than the observed barrier for capture, it 
cannot be associated with the capture cross section of the defect. It must be partially 
due to a macroscopic barrier for the carrier to be captured. In fact, since the buried 
layer always remains within the depletion layer, part of the observed barrier may be due 
to temperature dependence of limited supply of holes at the buried layer during filling 
pulse. In addition, this trap most probably has considerable lattice relaxation associated 
with it. Separating out the contribution to activation energy from lattice relaxation and 
macroscopic barriers is not simple in the present configuration of the samples. 
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L n<icr stich circumstances, defects in p-t}T)e Si seemed to have stabilized quicker than in 


n-type Si. 

Figure 7.13 shows a set of TSCAP spectra for as-implanted and annealed samples. In 
each case, the dotted line curve is obtained for TSCAP spectra without zero bias filling of 
the traps and solid lines are after trap filling with zero bias. In both the cases, capacitance 
is recorded only during the heating cycle with controlled heating rate. For as-implanted 
and 400°C annealed samples, presence of multiples steps in TSCAP signal indicates the 
presence of several discrete traps. However, for 600°C annealed sample, only one large 
step is present indicating presence of one type of defect with large concentration. This 
study shows that in p-type Si, all the defects are not annealed out at 600°C in contrast 
to n-type Si where most of the damage related defects are annealed out 

DLTS measurements are also done on these samples to survey residual/secondary 


defects after annealing. Figure 7.14 shows DLTS spectra for as-implanted, and for samples 


annealed at different temperatures (160°, 400°, and 600°C, 30 minutes). Clearly, all the 
spectra are different from each other for same rate window. This indicates formation 
of new defects/complexes even for low temperature annealing. The defect present after 
600°C annealing occurs in very large concentration and thus causes large step m TSCAP 

spectra. 

In this annealed sample, the major trap was found to have an activation energj' of 
E.+0.51 eV with capture cross-section of ~ 10- cm^ The peak related to this level in 
DLTS was found to occur at a lower temperature than a simUar level reported by Benton 


et al. who studied defect evolution from point defect to extended defect with high dose 
implantation and subsequent annealing [19]. It is worth mentioning that in the report 


of Benton ei al. for high dose implanted annealed sample DLTS peak occurring at lower 
temperature is not analysed though it was dWnctly clear. We believe that these defects 

are of same origin i.e. interstitial clusters. 
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Figure 7.13: TSCAP spectra for (a) as-implanted, (b) 400°C annealed and (c) 600°C annealed 
samples with zero bias filling of traps (solid line). Dotted line represents spectra without zero 
bias filling. 
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Figure 7 14- Comparison of DLTS spectra for (a) as-implanted, (b) 160”C annealed (c) dOO'-C 
and (d) 600»C annealed Ai+ irradiated p-Si. The annealing time is 30 minutes in eadr 
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In the literature, formation of secondary defects upon annealing has been reported 
[179] for implanted p-Si. Such continuous change in DLTS spectra for different annealed 
samples is indicative of formation of defect clusters of varying size. The evolution of point 
defect cluster is only beginning to be studied at present and Benton et al. [19] term 
it as the study of ’’missing links” between intrinsic point defects and extended defects 
that they finally give rise to on annealing. In a study of low temperature annealing of 
defects in electron irradiated p-Si, formation of vacancy-oxygen complexes was reported. 
We have not seen any divacancy related level in our heavy ion irradiated p-Si. Moreover, 
temporal changes in DLTS spectra with minority carrier injection related effects has also 
been reported [180]. We have not observed any such effect. 

7.6 Discussion 

The electrical defect signatures described in preceding sections are different from those 
point defects normally associated with complexes of vacancies and interstitials with impu- 
ritites observed for low dose ion or electron irradiation. The difference could be ascribed 
to possible agglomeration of defects leading to distinct electrical signatures in case of high 
fluence irradiation. There have been extensive TEM studies in the fluence regime ilO^^ 
cm“^ of ion irradiation in Si, after high temperature annealing [56, 58]. These studies 
reveal formation of interstitial chains and related extended defects commonly referred to 
as {311} defects. At very high doses, the formation of amorphous region and the kinetics 
of amorphous/crystalline interface has been explained on the basis of accumulation of 
di-interstitial pairs [55]. Though the involvement of interstitial agglomeration has been 
beyond doubt from these studies, the electrical signatures of defects arising from such 
clustering have not been studied except for a recent report by Benton et al. [19]. In 
our case, the fact that these defects are detected much beyond the ion range is in itself a 
strong indication that interstitial clusters are responsible for the observed defect spectra. 
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StiKiuvs on clust(>r formation using molecular dynamics predict formation of clusters of 
varu'd siz(\s which m principle should reflect in broadening of electrical activation ener- 
gies. Our I ( Suits show that most spectra is dominated by point-defect-like trap signatures 
with ielati\ely small amount of broadening in energy, which could be either due to nar- 
row distribution of stable cluster sizes or inhomogeneous environment in the vicinity of 
the defect. The defect spectra goes through complex changes with annealing, and finally 
yielding a stable defect with large concentration after 600°C annealing as described in the 
previous section. The process of agglomeration seems to be more complicated in this case. 
In contrast, for the n-type case the major defect only displayed progressive deepening of 
activation energ}' with annealing. We do not' observe signatures of extended defects ei- 
ther in terms of broad structures in the defect spectra or logarithmic filling of occupancy. 
The annealing temperature of 600°C may not be sufficient for the formation of extended 
defects. 

The major defect is found to be responsible for compensation in the damaged region. 
Origin of the minority carrier trap peak can be explained on the basis of mild inversion 
occurring within depletion region. From analysis of minority carrier trap related peak, 
it was found that the time constant of emission of this peak and of the major trap in 
n-type are identical at same temperature. Thus it can be concluded that the dominant 
majority carrier trap in ion damaged n-Si and the minority carrier trap in as-implanted 
p-Si are of common origin and it acts as an efficient recombination center. The dominant 
majority carrier trap in p-typs Si is in all likelyhood related to di-interstitial or higher 
order interstitial clusters. 

7.7 Summary and Conclusions 

In this chapter, we have studied defects in high dose ion implanted p-type Si using a 
variety of steady state and transient spectroscopic techniques. As compared to n-type Si 
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case, the ion induced defects in p-type Si manifest slightly different features in spite of 
many commonalities. Results and conclusions based on C-V, DLTS, TSCAP, CC-TATS 
measurements on as-implanted and annealed samples can be summarized as follows : 

• C-V studies show presence of damage induced defects at a depth larger than that 
predicted by TRIM calculations. However, this depth is not as large as in n-type 
case. Migration and subsequent formation of defect cluster at a larger depth can 
explain such observation. 

• The DLTS spectra of as-implanted sample shows three distinct majority carrier trap 
related peaks and one minority carrier trap related peak. Well known CjOi center 
occurs in small concentration. The new major damage related trap has small energy 
broadening which has been ascribed to narrow distribution of stable cluster sizes or 
inhomogeneous environment in the vicinity of the defect. Minority carrier related 
trap peak is explained on the basis of availability of electrons in the depletion layer 
due to inversion or near intrinsic character of damaged layer. 

• By using CC-TATS, the activation energy of the dominant trap has been found to 

be -1-0.52 eV with a large capture cross-section {a ~ cm^). 

• Capture transient was found to be extremely slow and thermally activated. Using 
both capture DLTS and varying filling-pulse-width techniques a macroscopic capture 
barrier of 0.66 eV was found. This is partly due to large lattice relaxation associated 
with the state and partly due to macroscopic barrier limiting hole supply. 

• The defects in p-Si are less stable thermally. They form new secondary complexes 
upon annealing. No signature of extended defect formation is observed upon an- 
nealing upto 600°C. 
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• Estimation of time constant for the negative peak in CC-TATS spectra reveals that 
this peak is same as found in n-type Si. Thus the damage related major defect is 
an (dficient recombination center. 


Chapter 8 

Summary and Conclusions 

8.1 Summary and Conclusions 

Traditionally electrically active defects induced by irradiation in semiconductors have 
been studied for light particles and, to a lesser extent, for ions with (i) energies in the keV 
range, (ii) low dose, or (iii) after removal of damage through annealing. However, there 
has been a recent thrust, motivated by technological applications, in studies of defects 
due to high dose irradiations of MeV heavy ions. This thesis deals with several aspects 
of electrically active defects in deep buried layers in Si produced by MeV heavy ions 
using mainly capacitance based techniques. The work has involved critical evaluation of 
existing spectroscopic techniques and their suitable extension needed for studying heavy 
ion damaged semiconductors. These studies reveal several significant aspects concerning 
defect migration and clustering, electrical identification, effect of degree of disorder, and 
defect relaxation mechanisms. 

In order to separate out effects due to heavy ion damage alone firom electrical activity 
of implanted species, we have used noble gas Av^ ions, with energies around 1.45 MeV 
and doses cm“^ to produce buried layers in both n and p-type epitaxial silicon. 

Samples implanted with heavier Au'*' ions at 4.5 MeV and dose 5x10 cm , but not 
annealed to avoid electrical activation, have also been studied for comparison. Irradiation 
doses are chosen to be below amorphization threshold but comparatively higher than used 
by most other workers. Conventional Schottky barrier diodes are fabricated for charac- 
terization using capacitance based techniques. Both as-implanted and low temperature 
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carriers, large magnitude emission transient and its premature termination. These features 
are shown to be consequences of trap filling giving rise to large negative charge density 
within the depletion layer and are discussed in the light of crossing between quasi-Fermi 
level and the major trap controlling these phenomena. DLTS, in spite of its inadequacy 
in this situation to provide reliable defect parameters through standard analysis, has been 
useful as a survey technique, and does provide helpful guide to further experimentation. 
In case of n-Si, both DLTS and thermally stimulated capacitance (TSCAP) spectra in the 
range 90K-320K reveal occurrence of mainly two discrete trap levels of which one is the 
well known divacancy center (VJ Ec-0.42 eV) and the other is a newly found midgap 
trap. This new trap level is the dominant trap occurring in very large concentrations in 
all our lyleV irradiated samples and is responsible for controlling hysteresis in C-V and 
space charge limited conduction in forward I-V characteristics. 

(ii) Nonexponentiality : Isothermal Spectroscopic Study 

Transient measurements for time domain analysis (as in TATS) performed in constant 
capacitance mode were found to be most useful since they overcome known sources of 
nonexponentiality due to large trap concentration and series resistance. TATS lineshape 
corresponding to the major trap is seen to be broad for low temperatures and, narrow and 
skewed at high temperatures. The peak shape narrowness is linked to occurrence of pre- 
mature termination of transients which in turn is explained on the basis of time evolution 
of quasi-Fermi level and its crossing with the trap level within the depletion layer. It is 
possible to estimate the true time constant and peak height associated with these skewed 
peaks from a careful fitting of the useful part of the lineshape. The signature of the dom- 
inant trap is obtained from such analysis for various samples. The emission energy of the 
dominant trap is observed to vary from Ec-0.49 eV to Ec-0.56 eV depending on processing 
conditions such as low temperature oven annealing and irradiation dose. This sensitivity 
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of the defect parameters seem to be characteristic of preamorphized regions and is dne to 
varying degree of relaxation of disorder in its environment. At low temperatures, in the 
absence of skewness in lineshape due to premature termination, the damage related peak 
for high dose samples shows a Gaussian broadening in energy with FWHM of 25 meV. 
This !>roadening is disorder induced since it is observed to be only 6 meV for low dose 
samples. The capture cross-section values are typical of neutral trapping center. The 
possible origin of this defect h^ been discussed, most likely candidates being dangling 
bonds, higher order vacancy complexes or interstitial clusters. However, consistent with 
our observation of involvement of migrating species and in accordance with mounting ev- 
idence from recent studies on involvement of interstitials, di-interstitials or their clusters 
seem to be responsible for the dominant damage related trap. High temperature furnace 
annealing studies show near complete removal of the dominant trap on heat treatment at 
600 C for 30 minutes. In annealed samples minority carrier related traps could also be 
observed possibly due to mild inversion in the damaged layer. 


D. Charge Redistribution 

Trapping kinetics of multiple traps has been studied by monitoring relative occupan- 
cies of multiple traps after allowing filling of traps for different duration over 5-6 orders 
of magnitude in time. Surprisingly we observe two distinct mechanisms which could be 
isolated at different temperature regimes in different group of samples. At lower tem- 
peratures, filling time dependence of TATS peaks corresponding to multiple traps show 
nonmonotonic change in trap occupancy. With filling time evolution, it was possible 
to recognize that the peak attributed to divacancy alone, in reality, constitutes contribu- 
tions from two distinct emitting centers. The resolution of peak structure has been clearly 
demonstrated using higher order implementations of TATS. The progressive evolution of 
occupancies over many orders of magnitude in time reveals coupling in capture kinetics. 
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We show that this coupled kinetics is due to charge redistribution among multiple traps, 
in which deeper states gain at the cost of shallower states with time. Broad experimen- 
tal features in occupancy as a function of filling time could be reproduced by solving 
coupled rate equations for three independent traps with charge redistribution constraint. 
However, closer match with experiment is obtained when in addition a configurational 
interdependence of two traps is assumed. The mechanism of charge redistribution had 
been invoked earlier to explain charge relaxation in amorphous silicon, though it was 
considered controversial in that context due to lack of unambiguous experimental proof. 
This mechanism is expected to be operative whenever number of trapping sites are larger 
than availability of carriers for capture. 

E. Novel Defect Relaxation 

Evidence of a novel defect relaxation mechanism was found by studying filling time 
dependence of occupancies for low dose implanted samples (Au'*' at dose 5 x 10^ cm"^) and 
for partially relaxed (400'’C annealed) Ar"*" irradiated samples - presence of a preamor- 
phized region with comparatively lesser degree of disorder being common to these set of 
samples. The observed phenomenon could most conveniently be studied for temperatures 
above 270K. When filling time is increased at a particular temperature, TATS emission 
spectra shows a distinct new peak whose time constant is strongly dependent on the du- 
ration of filling. The time constant of emission shows a power law dependence on filling 
time. This dependence of emission characteristics on history, in terms of waiting time, 
is a strong indicator of a non-Markovian process of relaxation occurring during capture. 
The height of the new evolving peak remains almost constant for all filling times. Further, 
changes in ionization energy seem to involve an entropy-like term since the emission time 
constant is nearly temperature independent for a particular filling time. All these features 
suggest that the relaxation mechanism follows a hierarchically constrained dynamics sim- 
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ilar to those proposed for glassy systems. A closely similar defect relaxation mechanism 
has been invoked to explain capacitance transient data in hydrogenated amorphous silicon 
though it is more involved in that case due to presence of broad density of states. 

F . Comparison With p-type Silicon 

Studies on p-type Si with high dose irradiation showed qualitatively similar features 
as that of n-t}pe Si. From spectroscopic analysis, in addition to well known traps a 
new major trap and a minority carrier related trap were found. Major peak in this case 
is due to damage related trap and the DLTS lineshape is broader than peak expected 
from exponential transients due the strain in lattice. Proper evaluation of trap emission 
time constants, capture cross-section, trap concentration and spectroscopic lineshape were 
carried out using constant capacitance TATS analysis. The major trap is found to be 
located at Eu+0.52 eV having a large capture cross-section of ~ 10“^^ cm^. 

Carrier capture rate at this major trap is seen to be surprisingly slow and thermally 
activated. From DLTS-like processing of capture transients as well as from variable filling 
pulse width TATS measurements, capture barrier was estimated to be 0.66 eV which is 
quite large compared to the emission barrier mentioned above. A part of capture barrier 
may be due to a macroscopic barrier for holes to reach the trap levels at the damaged 
layer. With a high macroscopic barrier, a mild inversion or injection from near-intrinsic 
region are probable sources of filling of minority carrier traps which show up as negative 
peak in CC-TATS spectra. The time constant of the negative peak matches with that of 
damage related majority carrier peak in n-Si. Thus the dominant majority carrier defect 
in n-type and the minority carrier trap in p-type Si seem to be of common origin acting as 
an efficient midgap recombination center in heavy ion damaged Si. The majority carrier 
peaks in p-type also seem to have origin in interstitials and their clusters. 

In contrast to the case of n-Si, high temperature furnace annealing of damaged p-Si 
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upto 600^0 showed considerable change in defect spectra with appearance of new traps. 
The evolution of defects would require more systematic studies than envisaged in this 
work. 


8.2 Recommendations For Further Work 

Our work may be viewed as an attempt to lay a foundation for studies in defect dominated 
materials, such as heavy ion damaged Si, by overcoming several hurdles in the use of 
suitable techniques and interpretational dilEculties. There are several obvious extensions 
of our work some of which are listed below. 

1. Many of the interesting questions of damage and defect mechanisms require com- 
plimentary characterization techniques. Particularly, use of structural tools such as 
transmission electron microscopy (TEM), Rutherford backscattering spectrometry 
(RBS) channeling and positron annihilation spectroscopy (PAS) would give more 
coherent understanding of the damage related defect. Similarly, magnetic resonance 
studies such as EPR, ENDOR along with electrical studies would give details of the 
defect structure and symmetry for proper identification. 

2. Optical characterization techniques such as photocapacitance based spectroscopy 
would simultaneously be useful for such compensated material to develop better 
understanding of the defect. 

3. More systematic study of the defect in buried layers using our electrical charac- 
terization techniques should be done by varying different parameters (ion dose, 
implantation, temperature etc.) for ion implantation. Also, a more systematic and 
thorough study of evolution of defects with annealing is required to understand re- 
lationship between clustering of defects and their corresponding electrical activity. 
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During the course of this work, it was found that forward biasing of diodes some- 
times led to changes in defect spectra. This aspect needs further systematic and 
controlled studies. This may reveal that relaxation of the damaged layer could be 
induced by recombination events. 

4. During implantation, changes in defect spectra can be monitored in-situ to study 
defect kinetics at shorter time scales following irradiations into finished devices. 

We believe that extension of this work along these directions would open up many 
more possibilities in this regime of defect clustering in ion-implanted semiconductors and 
similar defect dominated materials. 
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heavy ions. There axe several simulation packages such as TRIM (Transport of Ion in 
Matter) [15], MARLOWE (based on Binary collision approximation programs) [16] 
which work remarkably well for low energies in predicting range and profile of ions and 
displaced atoms. Such computer programs normally use ballistic regime approximation 
and do not take into account defect processes. 

1.4 MeV Implant Induced Damage and Defects : A 
Brief Review 

The energy lost by the ions can be sufficient in displacing host atoms, which in turn 
can displace others resulting in collision cascades. Normally only a nuclear interaction 
result in displacement of host atoms. A host of different types of defects can results 
starting from simple point defects to their clusters, complexes and amorphizcd zones due 
to agglomeration. Radiation damage is produced only if the energy loss is larger than 
the displacement energ}’ of the lattice atom (typically ~ 15 eV for silicon). Therefore, 
the damage distribution can be expected to lie closer to the surface than that of the ion 
profile. An example of the difference is shown in Fig. 1.2, where vacancy and ion profiles 
are obtained using TRIM simulations. For low energy implants, the peak of damage 
profile generally occurs at ~ 0.7-0.8Rp with standard deviation of ~ 0.75ARp, where Rp 
and ARj, are ion range and straggling respectively. It has been observed that a simulation 
package such as TRIM normally yields ion profiles in close agreement with experiments, 
but the damage profile can be very different than predicted for heavy ion implants in MeV 
range. The final damage profiles in these cases is dependent on details of defect processes 
such as diffusion and agglomeration, in addition to cascade geometries. This is the reason 
why the results of the simulational packages cannot be taken at their face value at present 
for heavy ion implants in the MeV range. 

As described by Seidel [17], an implanted ion comes to rest in ~ 10~^^ s, which is 
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Figure 1.4: Schematic of the recovery of several interstitial related defects as studied by EPR 
(after ref. [44] and references therein). The relative intensity (number) is irrelevant between 
defects as annealing data are from different sources. 

in this thesis, we note that silicon interstitial generated by Ar"*" ion implantation has been 
shown to cause an order of magnitude increase in coefficient of phosphorus diffusivity at 
temperature less than 800°C . Argon ion implantation were shown to be source of long 
lasting silicon interstitials [45]. The current views on the life of an interstitial from the 
time of its generation has been succinctly summarized in a recent paper [46]. It has been 
shown that for sufficiently high damage levels, the resulting excess number of interstitial 
atoms are approximately equal to the imbalance in the interstitial vacancy concentrations 
resulting from an extra, implanted silicon atom. This is commonly referred to as ’+1’ 
model [47] which suggest that recombination of interstitials and vacancies is the dominant 
process operative at the initial stages of their generation. 
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mechanisms in presence of disorder. By studying such phorinmcna in a regime of 
disorder below amorphization, one hopes to avoid difficulties encoimterc’d in rhe 
more complex situation of amorphous silicon. 

Our aim in this work has not been to carry out comprehensive comparative stud- 
ies through variation in implantation parameters, rather the attempt has been to lay 
down a foundation for such studies using capacitance based techniques by exidoring ba- 
sic phenomena and possible interpretational hurdles. In other words our effort luis been 
to hone the capacitance based electrical characterization techniques to study electrical 
manifestation of defect phenomena due to MeV heatw ion implantation notwithstanding 
complexities due to heavy damage. More generally, these studies would be helpful in the 
study of complex inter-relationships between defects, damage and disorder. M(>\’ heavy 
ion damaged silicon provides a challenging setting for such studies. 
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scheme extensively in transient measurements to keep the sampling volunie same and to 
avoid non-exponentialities in transients due to large concentration of traps. X('(’(’s.sar}' 
instrumentation for implementation of this circuit is discussed in the next chapter. 

In practice, series resistance can be temperature dependent also. In a subsequent 
section, we will discuss way of overcoming such problems for analyzing capacitaiu'e tran- 
sients. 

2. 2. 3. 2 Excess leakage current 

Schottky barrier devices occasionally show large leakage current. After ion implantation, 
the device may become more leaky due to presence of recombination centers. The fis- 
sumption in the conventional profile equation is that the voltage is metisured across the 
space charge region only. If, however, excess leakage current flow an appreciable voltage 
can be developed across the quasi-neutral regions. This voltage is automatically included 
in the recorded voltage introducing errors in measured C-V characteristics [89|. 

2. 2. 3. 3 Minority carrier and interface trap 

When a Schottky barrier diode is reverse biased, the thermally generated electron-hole 
pairs are swept out of the space charge region and leave through the ohmic contact of 
the device. In case of MIS devices, thermally generated minority carriers drift to the 
insulator/ semiconductor interface to form inversion layer. The device is unable to remain 
in deep depletion and doping profile measurement will be in error. When the voltage is 
varied rapidly, inversion charge do not get time to appear. Cooling the device to low 
temperature works well to reduce the effect of minority carrier generation. Effect of 
interface traps is also of major concern in case of MIS structures. Its contribution is also 
reduced when the device is cooled. 
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2.6.2 Effect of series resistance and leakage current 

We have already discussed the effect of series resistance on the cnr^ariran T' transients. 
This can be extended to DLTS as well. Using Eqn. 2.12, measured DLTS .signal AC,,, can 
be writen as 

12.22) 

where Q = atr^C. This simple form of equation is valid for small change of Q with time. 

According to Eqn. 2.22, DLTS signal would vanish for Q = 1 and its sign is reversed 
for (3 > 1. Broniatowski tt al. [101] have proposed a method of detecting tlu' presence 
of series resistance in DLTS and a correction scheme, based on varying circnii impedance 
by means of additional series resistance. In practice, this is cumbersome to implement 
due to variation of and C with temperature. The distinction between negative peak 
due to minority carrier trap and peak due to series resistance can be made by noting the 
change in crossover point in temperature due to extra series resistance added [109. 102]. 
If no change is observed, the sign reversed DLTS peak can be attribtited to the presence 
of a minority carrier trap. 

The lineshape distortion in DLTS due to series resistance gives modified trap param- 
eters. Lineshape analysis of DLTS for different Q values show that larger the Q value. 
FWHM and peak height of DLTS will go down and peak will be shifted to higher tem- 
perature side for chosen window. Figure 2.4 shows a set of simulated DLTS .spectra for 
different Q values. For Q>1, peaks of opposite signs can be noted. This may be mistaken 
as minority carrier trap related peak. For temperature dependent series resistance, peak 
height will be a complicated function of resistance. 

Reduction in DLTS peak height with smaller rate window can be due to leakage current 
effect [107] as well as series resistance effect. As discussed in section 2.3.2, time constant 
will be changed and extracted trap parameter can be very different from actual parameter. 
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Figure 2.8; Typical PITS spectrum of SI-GaAs sample with ti =0.141 s, t2=0.279 .s. .showing 
negative peak along with broad positive peak. 


a broad positive peak prior to the negative peak and has been seen in most samples. The 
shape of the negative peak in PITS spectra is highly distorted owing to the presence of 
this broad positive peak. However, significance of this feature in relation to the negative 
peak has not been discussed in the literature. We also see a well formed positive peak at 
lower temperature in the PITS spectra which is discussed later in this section. 

In our experiments, as do other reports, we have seen that the strength of the transients 
is highly dependent on temperature for a constant filling time leading to distortion of 
PITS spectra. Moreover, most transients are non-exponentials to varying degrees making 




Chapter 3 

Experimental Details : Sample 
Preparation and Measurements 

3.1 Introduction 

In this chapter, details of sample preparation, diagnostic techniques and experimental 
set-up made in house are given. The diagnostic studies involve capacitance and current 
based steady-state or transient measurements in the temperature range of 77K-330K. 
Most experimental setups and their automations for various measurements were devel- 
oped in house using a few commercial instruments as basic modules. Flexibilities of 
in-house developments in the hardware and software were important for this work, since 
many designer experiments were needed to be carried out to decide suitable mode of 
characterization and for thorough appraisal of cross-sensitivities. 

3.2 Sample Preparation 

A typical sample preparation process involves following steps : (i) surface cleaning of 
Si wafers, (ii) irradiation with MeV ions (Ar, Au) on Si wafers and finished Schottky 
devices, (iii) vacuum annealing of implanted wafers or oven heating on finished devices 
for low temperature annealing, (iv) Schottky diode formation with back ohmic contact, 
(v) mounting the device on a header for taking epoxy contacts. Most of the studies in this 
work were done on as-implanted sample where step (iii) mentioned above is omitted. The 
2 MeV Van de Graaff accelerator facility here at IIT, Kanpur, India, was used for Ar+ 
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Figure 3.6: Biasing circuit for application of both forward bifis and reverse l)i;is. 


3.3.2 Capacitance vs. voltage (C-V) measurement 

A Boonton capacitance meter ( Model 72B with 50 fis response time] operated at 1 
MHz with a 15 mV test level was used for all capacitance measurements. Programmable 
voltage is generated with the help of a plug-in AD/DA card (Model PCL 718). It is a high 
performance, high speed, multi-function data acquisition card for IBM PC/^^T/AT. It has 
two analog outputs (range 0 to +5 V, with a resolution of 12 bit and 5 iis set tling time). 
For changing the polarity of the voltage, an opamp inverting circuit is implemented. A 
simple adder circuit whose schematic is shown in Fig. 3.6, is used to get voltage in the range 
of -5V to -I-5V. In most of the cases, voltage is changed from high value to low value while 
monitoring capacitance at each voltage. This prevents effect due to dynamical change of 
capacitance due to trap de-occupation at a temperature where trap emission time constant 
is comparable to the voltage sweep rate. For all C-V hysteresis measurements, voltage 
is swept at a predetermined rate from high voltage to low voltage and then from low 
voltage to high voltage. The diode is kept under reverse biased condition to start with. 
The applied voltage and capacitance meter analog output are digitized through two input 
channels of Keithley high speed voltmeter (K194A) which can store 64 KB data. Different 
ranges of capacitance meter gives different analog output which is IV or 3V maximum if 
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Voltage (V) 


Figure 4.1: Typical C-V characteristics at T=300K for n-Si Schottky diode irradiated with 
various doses of Ar''" ions. 

have been irradiated at room temperature with various implant doses in the range of 5- 
10x10^^ cm“^ of Ar+ ions and 5x10® cm“^ of Au"*" ions. Typically in a C-V measurement, 
bias voltage is changed from high reverse bias to zero bias so that time dependence of trap 
occupancy, if any, does not affect the steady state capacitance measurements assuming 
that carrier capture rate is very fast compared to rate of change of voltage. 

Figure 4.1 shows a typical 1/C^ vs. V plot at room temperature for unimplanted and 
Ar"*" ion implanted samples with two different doses. It can be noted that for unimplanted 
samples, the curve is linear in the measured voltage range indicating uniformity of shallow 
doping concentration. However, for implanted samples there are distinctly two regions of 
which the linear region (A) occurs for higher voltages and the flatter region (B) occurs for 
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Voltage (V) 


Figure 4.8: Forward and reverse I-V characteristics of Ar'^ implanted Si at room temperature. 

is low enough for unproblematic capacitance measurements as shown in the last section. 
However, note the sudden jump of forward current at a particular voltage. 

In Figure 4.9, we have plotted the log(I) vs. log(V) for Ar'*’ ion implanted sample 
annealed at 160°C for 30 minutes. This plot shows a near ideal case of crossover from 
ohmic conduction (I oc V) to square law space charge limited conduction (I oc V^) regime 
at the voltage corresponding to trap filled limit (Vtfl)- The nearly vertical rise in current 
accompanying the trap filling is an indicator that the associated defect is a discrete level 
[91]. In this case, the same defect controls hysteresis in C-V measurements due to its large 
concentration and this is the dominant level in ion damaged buried layer which acts as 
a near intrinsic region. Thus the structure under study is akin to p-i-n diode where p is 
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dynamic simulational studies of combined processes of defect generation, migration and 
clustering during and after implantation. Their atomistic simulation show migration of 
free interstitials and clustering at room temperature even several hours after initial fast 
recombination processes have died out. They also show that the depth distribution of 
remnant interstitials peak at approximately twice the implant range. Defect clustering is 
known to be more favourable in heavy ion irradiation due to occurrence of high density 
of collisional cascades in comparison to light ions. Similar efllects are also expected for 
high dose implants [19] . Therefore there is a growing body of evidence showing that 
the ion induced damage acts as a source of migrating interstitials which can combine to 
form clusters at distances beyond the damage location. Though the electrical activity 
of point defects and extended defects in implanted material have been studied in detail, 
there is a lack of similar studies on defect clusters. Its only recently, Benton et al [19] 
have reported electrically active traps associated with second and higher-order interstitial 
clusters in residual damage of Si implanted with high fluences. It appears reasonable 
to conclude that the dominant electrically active traps observed in our experiments also 
originate from interstitial clusters. We postpone a discussion on origin of the major defect 
till the next chapter where we report spectroscopic characterization of electrically active 
defects. The sharp trap profile at distances much larger than ion range would be fully 
consistent with the mechanism of migration and eventual cluster formation of interstitials. 
In fact, our results constitute a direct experimental evidence of such a mechanism. 

The fact that the effective electrical interface moves with changes in process variables 
such as dose and annealing also supports occurrence of a mechanism involving migration 
and clustering. With increase in irradiation dose, the interface between trap dominated 
and trap-free regions is expected to move deeper into the sample as there would be more 
abundant supply of interstitials from the damaged region. Similarly, the effective interface 
is observed to move towards the surface with increase in annealing temperature. Though 
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Though large concentration of defects are produced in ion-atom collisions, most of the 
primary defects, namely vacancies and interstitials, are unstable at room temperature and 
the} foini di\acanc} and other higher order complexes. Divacancy is known to have three 
charge state configurations which introduce different energy levels (Ec-0.42 eV, Ec-0.23 eV 
and E„-f-0.25 e\ ) in the bandgap of semiconductor. However, we do not observe the peak 
corresponding to Ec-0.23 eV related charge state configuration in the DLTS spectrum. 
This is mainly due to low production rate of - 0.23 level as compared to the - 0.42 
level in heavy ion irradiated n-Si. Note that the presence of a broad peak HI may cause 
suppression of other majority carrier related peaks which are of smaller concentration. In 
a detailed study of divacancy acceptor level profile in MeV ion irradiated n-Si, Svensson 
et al. [20, 143] reported that concentration of E'c-0.23 eV level is significantly lower 
than that of £'c-0.42 eV for high dose and heavier mass implants. They observed that 
concentration of the two levels are identical in tail region towards the surface while it is 
very different in the damage peak region. They argued that for heavy ions the rate of 
energy deposition is higher and A'a centers occur in lattice regions which accommodate 
lattice strain. This lattice strain causes broadening in the DLTS spectra. Inspection of the 
lineshape of the peak V 2 reveals that it is broader than the peak expected from exponential 
transient. Detailed analysis of the lineshape in the context of capture kinetics is carried 
out using isothermal spectroscopic technique and is presented in the next Chapter. It is 
also known that a defect level with similar trap parameters i.e. E-center (P-V complex) 
may contribute to the observed broadening. We have examined possible contribution of 
this level from lineshape analysis of as-implanted and low temperature (160°C) annealed 
samples. In principle, for silicon with large phosphorus doping, contribution of this defect 
may become significant. Normally it is recognized by studying its annealing behaviour 
since the defect anneals out at temperature of ~ 300°C. 

The defect responsible for peak D1 shown in Fig. 5.4 occurs consistently in large con- 
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log(time/s) 


Figure 5.8: (a) Isothermal voltage transient measured in CC mode for as-implanted and 160"C 
annealed sample, (b) Corresponding CC-TATS spectra (solid line), with exponential peak fitting 
(dotted line) assuming apparent peak as the true peak. The broken line is fitted with an 
exponential transient considering only left side of the peak. 
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Figure 5.14: Comparison of DLTS spectra on as-implanted, 400°C and 600°C annealed samples 
for same rate window. 

low temperature annealing is most probably due to changing sizes of clusters. 

5.5.2 Furnace annealing 

Effect of vacuum annealing of the irradiated Si at higher temperatures are shown in 
Fig, 5.14. In contrast to as-implanted sample which showed V 2 related peak and D1 peak 
(solid line), 400^0 annealed sample shows a different DLTS spectrum consisting of two 
majority carrier peaks and one minority carrier related peak. The absence of \ 2 related 
peak is expected since they anneal out at ^ 350^0. The peak position of the major defect 
Dl is shifted towards lower temperature side. Though DLTS analysis can not be relied 
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laboratory for dealing with similar questions asked in the context of amorphous silicon. 
This is mainly because the defect density of states remain discrete in our samples while 
other conditions prevailing in amorphous silicon such as large concent rat i.m of defects and 
presence of disordered environment are still satisfied. 

Before we present our results, we provide in the next sec tion a icrief ; .t; Kui-'iiud to the 
above mentioned debate and the perspective in which our results must !»> ^ieweci. Wo go 
on to provide experimental evidence for two distinct mechanisms of charge n>laxation in 
varying sample conditions and different temperature regimes of meiisimnumts. 

6.2 Models of Charge Relaxation in Disordered Sil- 
icon 

The presence of continuously distributed states over wide energies is one of the charac- 
teristic feature of electronic properties of disordered semiconductor .There is a v;ist and 
rich literature in relation to the determination of density states, origin, tdft'ctive corre- 
lation energies, and mechanisms of charge relaxation in these materials. These studies 
have been most intensive in hydrogenated amorphous silicon (a-Si:H) where the energy 
distribution is attributed to the dominant deep defect ’D’ associated with dangling bond 
[156]. It is also appropriate here to mention that , in contrast to a-Si:H, there has been 
very little studies on DOS (density of states) in the gap of pure amorphous silicon created 
by MeV ion implantation [155]. Without going into the details of these studies, which 
have been reviewed from time to time, it suffices to state here that a principal source of 
most controversies have been the necessity of using complex deconvolution procedures to 
extract information from measured quantities. 

There have been several approaches to resolve inconsistencies in connection with dif- 
fering density of states in a variety of samples both doped and intrinsic arSi;H- There have 
been attempts to explain such a variation in DOS using the ’defect pool Tuodel’ in which 
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c2=400, e2= 
c3=80, e3= 



log(tf) 


Figure 6.5: Simulated occupancy behaviour of three trap levels by solving rate equations under 
the condition Nj-> n, (a) assuming three independent trap levels, (b) assuming that fastest and 

slowest t.wn pmrficmrjif.innc nf fVia a m'lron in 
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During the capture pulse, successive emission and capture events can drive the defect to 
relax continuously to more stable states if the polarization cloud is unable to follow this 
process rapidly. In other words a relaxation process much slower than emission or capture 
processes would have the memory of net duration of occupation. 

The driving force, in fact, can be due to entropic variations in configurational space. 
Such entropy driven metastabilities have been invoked in the past in the context of sev- 
eral defects in semiconductors known to have metastable configurations which can be 
accessed by charge state change, optical excitation and rapid heating and cooling cycles 
[146, 166, 167]. Hamilton et al. [146] have used such an argument to explain the existence 
of a critical temperature at which certain configuration becomes unstable giving rise to a 
discontinuous DLTS spectrum. Though that case is distinct from the ones described here, 
the significance of entropic terms in determining relaxation is to be noted. Further, it is 
also significant that most of such observations have been reported in electron or ion irra- 
diated semiconductors. This is to be attributed to occurrence of multiple or continuously 
varying configuration for an electrically active point defect in the environment of irradia- 
tion induced damage. The temperature dependence of electronic ionization energ}- has in 
the past been conveniently represented by associating an ionization entropy resulting from 
electron-phonon interaction [168]. The estimates of entropy change on ionization from a 
deep defect varies between 3 — 4/s in the temperature range of our observations [169, 170]. 
Hence entropic contribution to ionization can be well in the range of 100 meV at room 
temperature. The net free energy change of 180 meV in our case can indeed be due to 
changes in ionization entropy. It is to be noted that samples used in this work provide 
a convenient laboratory for further studies on details of this relaxation mechanism. In a 
group of samples, studied in this context, we observe qualitatively similar behaviour, as 
described here, but with the difference that multiple traps evolve with filling time. An 
example of Au"*” irradiated sample is shown in Fig. 6.10. These show that relaxations may 
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forward bias, only a small percentage of the current is due to injected electrons, so minority 
carriers are not usually seen. However, in the present study, there is a large damage due to 
heavy ion implantation. As the major trap (P2) causes carrier compensation and changes 
the band bending depending upon the trap occupancy, we ascribe this negative peak to 
changes in occupancy of minority carrier traps in presence of mild inversion caused by 
band bending in the damage region. A simple calculation of band bending at zero bias 
shows that the interface at the damaged region would be mildly inverted. This seems to 
indicate that a minority carrier trap can indeed be filled within depletion layer. 

Another possible mechanism of small minority carrier generation may be due to the 
nearly intrinsic nature of the damaged region where a large concentration of traps can 
effectively pin the Fermi level to a position near the middle of the energy gap. Such 
an effect has been reported previously [28] to occur in ion damaged Si. The equilibrium 
electron concentration in the damaged region acts as a source of electron to fill the minority 
carrier trap. The slow rate of carrier capture in the damaged region is due to the presence 
of small majority carrier concentration. Minority carrier trap related peaks were observed 
by Jackson and Sah [178] in boron implanted n-Si after annealing at ~ 400°C and 
observation was justified by possible intrinsic hole concentration in the damaged region. 

7.5 Thermal Stability Of The Defects 

In order to carry out a preliminary study of thermal stability of these damage related 
defects, we have annealed the irradiated wafers at 400°C and 600°C for 30 minutes in 
vacuum, and subsequently fabricated Schottky diodes on them. Some finished devices 
(irradiated) are annealed at IGO^C for 30 minutes in oven to study low temperature an- 
nealing effects. It is also important to mention that no room temperature transformation 
of the defects has been seen in p-Si, in contrast to n-Si case. Note that, our samples were 
stored at room temperature for few months before carrying out electrical measurements. 
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B. Evidence for Defect Migration and Clustering 

Unusual features observed in steady-state C-V studies on Schottky diodes of n-Si, 
containing damaged buried layers, strongly indicated trap controlled nature of the dam- 
aged layer. The damaged layer is found to be extended to several microns beyond the 
range of the ions or damage profile predicted by Monte Carlo simulations using TRIM. 
The experimental features in C-V curves have been simulated using model charge profiles 
taking into account crossing of Fermi level with a midgap trap within the depletion layer. 
The simulations suggest that the irradiation damage seem to produce two regions . (i) 
a compensated region with most dopants being deactivated, and (ii) a sharp negatively 
charged defect profile with large concentration close to the depletion edge, with the deple- 
tion layer being much larger than expected range of ions. Our experimental observations 
lend support to recent results of molecular dynamic simulations according to which pri- 
mary defects such as vacancies and interstitials migrate and form clusters and complexes. 
Occurrence of a single dominating trap in the damaged region is established from C-V 
and I- V studies. Dopant activation and annealing of the major defect was found to occur 
in the damaged region by 600^0 vacuum annealing. The nature of variations in charge 
profile with annealing shows that the effective electrical interface moves progressively to- 
wards the surface. 


C. Spectroscopies of Defects 

(i) Conventional DLTS : Survey and Limitations 

lyansieat response of diodes — g damaged layer stored interesting and nnnsual 
features wMch include occurrence of a ’spite’ in capacitance respond during capture of 
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